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NOTICES 


An annual sum of £250 is available for premium awards for papers, including 


Technical Notes, published in the Journal. 


Members and non-members of the 


Society are invited to submit papers on any aspect of aeronautics 


NEWS OF MEMBERS 


Dr. J. H. ArGyris (Fellow), formerly Reader in the 
Theory of Aeronautical Structures, has been appointed to 
the newly created Chair of Aeronautical Structures at the 
University of London, tenable at Imperial College. 

Group CAPTAIN A. F. BANDIDT (Associate) has been 
appointed the Handley Page Representative in Australia, 
New Zealand and the Far East. 

B. D. BLACKWELL (Associate Fellow) has been appointed 
Engineer-in-Charge of a new department for the develop- 
ment of turbine blades of the Engine Division of the Bristol 
Aeroplane Co. 

P. Bretr (Associate Fellow), formerly Chief Aero- 
dynamicist of de Havilland Propellers Ltd., has been 
appointed Chief Project Engineer of the Company. 

R. E. CaTTERWELL (Graduate) has been awarded the 
James Clayton Fellowship by the Institution of Mechanical 
Engineers. The Fellowship is tenable at the College of 
Aeronautics, Cranfield. 

J. R. Ewans (Associate Fellow), formerly Deputy Chief 
Designer, has been appointed Chief Designer to A. V. Roe 
& Co. in succession to Mr. S. D. Davies who recently 
resigned from that position. 

B. S. Massey (Associate Fellow) has been appointed 
Assistant Chief Designer of the Engine Division of the 
Bristol Aeroplane Co. 

SQUADRON LEADER F. Murpuy (Associate), formerly the 
Chief Production Pilot of Hawker Aircraft Ltd., has been 
appointed R.A.F. Liaison Officer (Flying) to Hawkers. 

J. W. Piccotr (Associate Fellow) has been appointed 
Apprentice Supervisor at the Dunlop Rim & Wheel Co. Ltd. 

R. A. READ (Associate), formerly of Cardiff Airport, has 
been appointed Aerodrome Commandant at Bahrain, 
Persian Gulf. 

Dr. A. E. RUSSELL (Fellow) has relinguished the title of 
Chief Designer of the Bristol Aeroplane Co. and now holds 
the post of Chief Engineer (Aircraft Division). 

B. SARAVANOS (Associate Fellow), previously with 
Handley Page (Reading), has taken up a post as Lecturer 
in Mechanical Engineering at Coventry Technical College. 

E. R. SISSON (Associate) has been appointed Technical 
Service Manager to Smiths Aircraft Instruments Ltd. 

S. S. TRESILIAN (Associate Fellow) has been appointed 
Assistant Chief Engineer of the Engine Division of the 
Bristol Aeroplane Co. 

G. W. UNbERWoopD (Associate), formerly Structural 
Test Engineer at Gloster Aircraft Co. Ltd., has been 
appointed Deputy Chief of the Company’s Structural 
Research Department. 

JEFFREY WATKINS (Associate Fellow) has joined Solar 
Aircraft Co. as Senior Aerodynamicist. He was formerly 
with Westinghouse Electric Corporation, Philadelphia. 


BRANCH AREA MEETING—BIRMINGHAM BRANCH 


An Area Meeting is being held at the Birmingham Branch 
on 23rd September 1955. Mr. G. C. Best, Head Perform- 
ance Engineer, Orenda Engines Ltd., Canada, will lecture 
on the “ Application of a Digital Computor to Gas Turbine 
Engineering.” The meeting will be held at the Chamber 
of Commerce, New Street, Birmingham, at 7.15 p.m. 


ASSOCIATE FELLOWSHIP EXAMINATION 


Entry forms for Part If of the Old Syllabus of the 
Associate Fellowship Examination may be sent in up to 
31st August 1955. Entries for the Revised Syllabus closed 
on 30th June. 


GRADUATES’ AND STUDENTS’ SECTION—VISITS 


A visit has been arranged to the Royal Air Force, West 
Malling, for Saturday 24th September. This is an Oper- 
ational Fighter Station, within easy reach of London, and 
applications for the visit should be made immediately to 
the Hon. Visits Secretary, Mr. P. D. Stewart, B.O.A.C. 
Operations Development Unit, London Airport, Hounslow, 
Middlesex. 

A visit has also been arranged to the National Aero- 
nautical Establishment at Bedford for Wednesday 12th 
October. It is hoped to arrange for a coach party to 
travel from London for this visit, and in order to complete 
these arrangements, early application should be made to 
the Hon. Visits Secretary as above. 


HIGH SPEED PHOTOGRAPHY CONGRESS 


The Third International Congress on High Speed Photo- 
graphy will be sponsored by the Department of Scientific 
and Industrial Research and will be held in London from 
10th to 15th September 1956 at Government Offices, Horse 
Guards Avenue, London, S.W.1. 

Those wishing to participate in the Congress, either by 
contributing papers or by supplying exhibits for the 
exhibition, should apply to:—The Congress Secretariat, 
Dept. of Scientific and Industrial Research, Charles House, 
5-11 Regent Street, London, S.W.1. 


MASSACHUSETTS’ INSTITUTE OF ‘TECHNOLOGY FOREIGN 
STUDENT SUMMER PROJECT 1956 


The M.LT. is offering a small number of Scholarships 
to university graduates in Science and Technology to enable 
them to attend the Summer Session of the Institute from 
4th June to 25th September 1956. No applications in 
Architecture and Town Planning can be considered. 

Both men and women are eligible for these swards. 
Candidates must be British citizens under 28 years of age, 
in good health, who hold recognised university degrees in 
technical fields and have had at least two years’ work in 
either industrial or academic research since graduating. 

Application forms may be obtained after Monday 22nd 
August 1955 from The Cultural Affairs Officer, American 
Embassy, Room 302, 41 Grosvenor Square, W. No appli- 
cations can be considered that are received at the Embassy 
after Friday 14th October 1955. 


RoyaL AERONAUTICAL SOCIETY STAFF REQUIREMENT 
The Society urgently requires an experienced shorthand- 
typist. Applications should be made to the Secretary giving 
full particulars of experience. 


THE LIBRARY 


The Library has been presented with “ Automobilia and 
Flight,” January-July 1909. If any members know of 
previous or subsequent numbers the librarian will be glad 
to hear of them. He would also like to have information 
on “Aeronautical Engineering” which succeeded the 
American “ Aeronautics” circa 1915. 
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Definitions of the Thrust of a Jet Engine 
and of the Internal Drag of a Ducted Body 


SUMMARY: The problems which occur in defining the thrust of a jet engine and the 
internal drag of a ducted body are considered, and formal definitions and names are 
given for the concepts considered to be of importance. 


(Although already published as A.R.C. Current paper No. 190, this paper is considered to be 
of such general interest that it is published here, at the suggestion of the Aeronautical Research 
Council, so that it may receive a more permanent form of publication and a wider circulation —Ead). 


1. Introduction 

Recent measurements, in flight, of the thrusts of 
turbo-jet engines"~’ have emphasised the need for a 
precise definition of jet-engine thrust, and also of the 
corresponding drag of the body housing the engine. 
Such definitions would be of interest to the 
thermodynamicist (who is concerned with the internal 
performance of the engine) and to the aerodynamicist 
(who deals with the air flow external to the engine, in 
connection with either aircraft performance estimation 
or thrust measurement in flight). Each of these 
considers certain forces to be important, and it is clearly 
necessary that the various terms employed by either 
should have a precise, unambiguous and mutually- 
agreed meaning. 

The problem of defining the thrust of a jet engine 
and the drag of the body housing is, of course, only 
a special case of the general problem of the flow through 
and round ducted bodies; it is fundamentally the same 
therefore as the problem of defining the internal and 
external drag, say, of a piston-engine cooling system, or 
a ducted heat exchanger. In all such problems it is 
usually convenient to consider separately the internal 
and external flows, and the thrusts and drags associated 
with them. The definitions of these forces must be con- 
sistent, however; that is, the definitions must be such 
that the overall longitudinal force on the body is given 
by the difference between the total thrust and the total 
drag. As the manner in which the longitudinal force 
can be divided into thrust and drag is to some extent 
arbitrary, it is not surprising that different definitions 
have been used in different contexts. 

Attempts to resolve this problem were made inde- 
pendently by Jakobsson“? and in an unpublished paper 
by Warren, who considered what concepts were 
important and suggested a number of names and precise 
definitions for the various terms involved. Warren’s 
paper was discussed at a meeting of the Performance 
Sub-Committee of the Aeronautical Research Council 
in January 1953, when it was agreed that the subject 
required further consideration. At the Sub-Committee’s 
request, the Council allowed a temporary Panel to be set 
up to consider the problem of defining thrust and drag, 
and also the related problem of the definitions to be 
used in the description and analysis of the drag of an 
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aircraft. The Panel, which met first on Ist May 1953, 
was composed of the following members, who were 
nominated either by the Performance or by the Engine 
Aerodynamics Sub-Committees. 


NOMINATED BY THE PERFORMANCE SUB-COMMITTEE 
Professor W. A. Mair (Cambridge University), Chairman 
Mr. D. W. Bottle (A. & A.E.E.) 

Dr. R. C. Pankhurst (N.P.L.) 

Mr. T. V. Somerville (R.A.E.) 

Mr. C. O. Vernon (S.B.A.C.) 

Mr. C. H. E. Warren (R.A.E.) 

Mr. E. W. E. Rogers (N.P.L.), Secretary. 


NOMINATED BY THE ENGINE AERODYNAMICS SUB-COMMITTEE 
Mr. S. Gray (Eng., R(a)M.O.S.) 
Mr. H. Pearson (S.B.A.C.) 
Dr. J. Seddon (R.A.E.) 
Mr. A. R. Howell (N.G.T.E.) attended some of the 
meetings. 


The present Report deals with the Panel’s investiga- 
tion of the problem of defining jet engine thrusts and 
concepts related to this; a Report dealing with the 
analysis of drag will be issued separately. The present 
investigation is limited to the case where the thrust acts 
in the direction of the undisturbed stream (i.e. along the 
line of flight), the jet exit plane being normal to this 
direction. The more complicated cases of an inclined 
thrust axis, an oblique jet exit plane*, or a jet exit plane 
upstream of the engine cowling exit with internal mixing 
between cooling and jet flows, have been borne in mind 
throughout the Panel’s discussions, but it was felt that 
it would be better to defer consideration of these more 
complicated problems until agreement had been reached 
for the simpler case. 

During the investigation it became clear that it 
might be convenient to introduce names for certain 
quantities having the dimensions of pressure. The 
Panel found, when reviewing these quantities, that some 
of the definitions relating to pressure, as laid down in 
the British Standard Glossary of Aeronautical Terms“, 


*If the nozzle is not square-cut and there is still some con- 
vergence in the final part of the nozzle, it is clear that some 
deflection of the issuing jet from the free-stream direction 
will occur, entailing certain fundamental difficulties. Even 
if there is no convergence, some deflection will probably 
occur when the jet is choked. 
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were wholly unacceptable for compressible flow. The 
Panel would have liked to have suggested more suitable 
definitions, but was unable to agree on the names for 
certain quantities, or even on the need for names in 
certain casesf. 


Notation 

The symbols in this list have been used in writing 
this report, but it is not recommended that they should 
all necessarily be adopted for general use. 


A station area (normal to the direction of the 
undisturbed stream) 
A; effective nozzle area quoted by engine manu- 
facturers 
Ds internal drag according to Betz’s formula 
(see equation (25) ) 
D; internal drag according to Jones’ formula 
(see equation (24)) 
M Mach number, given by V/ /(yp/p) 
m_ rate of mass-flow 
P static pressure 
P, local static pressure of the external fluid sur- 
rounding the jet at exit 
P,, pressure of the undisturbed stream 
V_ resultant velocity 
y ratio of specific heats of fluid 
@ thrust, given by equation (1) 
Qin, intrinsic thrust (see equation (4) ) 
gross Jones thrust 
Jones thrust } (see equation (8) ) 
Op, gross Pearson thrust ) (see equations (10(a) 
@,, net Pearson thrust } and (b) ) 
O,.s. post-exit thrust (see equation (5) ) 
Oy. pre-entry thrust (see equation (3) ) 
0,. gross standard thrust ) (see equations (6) 
©,., net standard thrust and (7) ) 
p  fiuid density 
@ angle between local direction of flow and 
direction of undisturbed stream (positive 
when streamtube is diverging) 


e denotes the exit station (see Fig. 1) 

f denotes the station, corresponding to the 
“effective area” of the jet, where the flow is 
parallel to the undisturbed stream (see 
Appendix III) 

g denotes the station at which the static 

pressure of the internal fluid attains the 

local static pressure of the external fluid 

denotes the intake station (see Fig. 1) 

w denotes a station at infinity downstream, at 
which the internal static pressure is uniform 
and equal to the static pressure of the 
undisturbed stream (see Fig. 1) 

c© denotes a station at infinity upstream (see 
Fig. 1) 


+In any case, a revision of the names for various concepts 
relating to pressure is outside the terms of reference of the 
Panel, and is doubtless a matter for a more widely repre- 
sentative body than one set up primarily to consider the 
definitions of thrust and drag. 


2. The Problem 


The resultant force on a stationary ducted body in 
a steady stream of infinite extent can be obtained by 
considering the change of momentum of the fluid cop. 
tained within an arbitrary surface completely surround. 
ing the body. It is desirable, however, to divide the 
resultant force into internal “thrust” and external 
“drag” forces, the general case being considered as one 
in which the internal flow contributes a thrust rather 
than a drag, in view of the importance of jet engines, 
The results obtained will be immediately applicable to 


the related problem of internal drag, simply by inter. 


preting thrust as negative drag. The internal thrust and 
external drag could be defined as the sums of the com. 
ponents in the appropriate directions of the aerodynamic 
forces on the internal and external surfaces of the body, 
assuming that these surfaces could be precisely defined, 
Such an analysis, although quite valid, would be of 
interest mainly to the structural worker, who is con- 


cerned with the forces on the various parts of the body } 


itself. It would be of less value, as will be seen later, to 
the aerodynamicist and thermodynamicist, who are 
concerned mainly with the fluid and the changes in its 


state, and who require definitions of the forces in terms } 


of the rate of change of momentum of the fluid that 


flows through the body (internal flow) and of the fluid | 


that flows round it (external flow). 


In considering changes of momentum in the fluid } 


flowing past a ducted body, it has sometimes not been 


clear whether the relevant static pressures should be} 


measured absolutely or relative to some datum pressure, } 
For a closed (unducted) body this difficulty does not | 
arise, since the subtraction of a datum pressure is } 


equivalent to integrating a constant pressure over a 
closed surface and does not contribute to the overall 
force. But in considering parts of the surface and of the 
flow, as is done with ducted bodies, the derived forces 


depend on the datum pressure used. For a body at rest | 


relative to the surrounding fluid the undisturbed fluid 
pressure acts on all surfaces, and the resultant force is 


zero. When there is relative motion between the fluid 


and the body this equilibrium is disturbed and the} 


surface forces which are of significance are those due to { 
the changes in surface pressure arising from the motion. | 


Thus in the derivation of these forces all pressures must 
be referred to the undisturbed static pressure as 4 


datum. This practice is adopted in the present Report. | 
The question of the division of the resultant force} 


into internal “thrust” and external “drag ” components, ) 


resolves into a choice of the upstream and downstream } 


stations between which the momentum relations are 
applied. 


If the reference stations are taken at the entry and| 


exit, the calculated values of these components will inf 
general vary with any change in entry or exit area,) 
although there would not necessarily be any change in} 


the overall flow pattern. Thus if the thrust were defined | 


in this way, the engine manufacturer would not be able 
to express the performance of his engine by a compact 
series of curves, applicable to all installations. More- 
over, he would not be able to measure the thrust with- 
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Boundaries of pre-entry 
stream Cube 


FiGurE 1. 


out traversing the flow at both the entry and exit 
stations. There would also be difficulty in defining the 
location of the entry station precisely. The use of 
the entry and exit stations is therefore particularly 
unsuitable for the definition of thrust, and some more 
generally applicable definition is required. 

If a new reference station is considered, upstream of 
the entry station, then the calculated force associated 
with the internal flow will differ from that actually 
exerted on the internal surfaces of the body. It is shown 
in Appendix II that this difference is equal to the longi- 
tudinal component of the pressure forces acting on the 
boundaries of the “ pre-entry stream tube” (Fig. 1). A 
corresponding treatment of the external flow would 
yield a force differing from that actually exerted on the 
external surfaces by an equal and opposite amount. 
Physically it can be said that the flow in the pre-entry 
stream tube, induced by the engine, influences the 
external flow and thereby the pressure distribution on 
the external surfaces. Thus a force, equal in magnitude 
to the longitudinal component of the pressure forces on 
the boundaries of the pre-entry stream tube, is rightly 
attributed to the engine thrust, though usually acting on 
the external surfaces. This is in fact the basis on which 
engine thrusts have usually been derived in the past. 
Treating the external forces as positive in the drag sense, 
then a corresponding external drag is obtained which is 
equal to the force actually exerted on the external 
surfaces (as would be derived from the actual pressure 
distribution) increased by a quantity which has some- 
times been referred to as the “ pre-entry drag.” 

A similar argument can be applied to the flow down- 
stream of the exit station by introducing the concept of 
a “post-exit stream tube” separating the internal and 
external flows. It is concluded that the fundamental 
definition of the thrust of a jet engine should be framed 
in terms of the stations at infinity upstream and infinity 
downstream. 

Consideration of the pre-entry flow presents no 
difficulties, but there is in fact no physical post-exit 
stream tube, the flow downstream of the exit being com- 
plicated by the transfer of momentum and energy 


Boundaries of equivalent 
post-exit streamtube. 


Diagrammatic representation of the flow through a ducted body. 


between the internal and external flows in the wake of 
the body. It has been customary in the past to assume 
that the local static pressure in the external flow in the 
region of the jet exit is equal to the undisturbed static 
pressure. In general, these pressures are not equal and 
it is necessary to take account of this in framing 
definitions of thrust. Moreover, the pressure of the 
internal fluid at the exit station may differ from that in 
the adjacent external fluid, and shock waves may be 
present in either the internal or external flows, or both. 
A major problem, therefore, is to devise practical 
definitions of thrust which take account of the complex 
flow conditions downstream of the exit, and which 
enable the thrust to be determined from practicable 
measurements. This is considered in the following 
Sections. 

Some cautionary remarks about central bodies and 
side intakes are given in Appendix IV. 


3. The Thrust and Related Quantities 


The thrust obtained as a result of choosing reference 
stations at infinity both upstream and downstream of 
the body, and by assuming the existence of an effective 
post-exit stream tube, is given by the equation (derived 
in Appendix IT) 


where A is the area of the station (normal to the free- 
stream direction), p is the density of the fluid, and V is 
the local velocity, parallel to the free stream at both 
stations. Suffix “w” refers to the station at infinity 
downstream (where the static pressure is uniform and 
equal to that of the undisturbed stream), suffix “oo” 
refers to conditions at infinity upstream, and the inte- 
gration is made over the mass of fluid flowing through 
the body. 

The thrust given by equation (1) may be divided into 
three components: 


where ®,,. and ©... are the pre-entry and post-exit 
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thrusts, associated with the pressure exterted by the 
internal flow upon the external flow ahead of and behind 
the body, and 0), is the intrinsic thrust*, associated with 
the force on the internal surfaces of the body. 

The components of 0 may be written as 


= (pV? cos? P; dA Po V.°A (3) 


= Pi dAi, . (4) 


(5) 


where suffices “i” and “e” denote the intake and exit 
stations respectively, and ¢ is the local inclination of a 
streamline to the free-stream direction. 


The term 
A 


in equation (3) is the “ram drag”; by continuity this 
is equal to mV, where m, is the internal rate of mass- 
flow at the duct entry. 

The pre-entry thrust depends on the conditions at 
the entry and in the undisturbed stream, which are 
assumed to be known. The intrinsic thrust depends in 
addition on the conditions at exit, which in principle are 
assumed to be measurable or calculable. It follows 
from equation (2) that the use of © depends on the 
successful determination of the post-exit thrust, and 
thus in any practical definitions of thrust, an attempt 
must be made to relate the conditions at infinity down- 
stream to those at the jet exit; this approach is con- 
sidered in Section 4. 

The presence of a post-exit thrust depends on the 
existence, in the region surrounding the jet exit, of 
pressures differing from that of the undisturbed stream, 
and the value of the post-exit thrust depends on the 
magnitudes of these pressures and the subsequent 
behaviour of the combined internal and external flows. 
Hence any assessment of the post-exit thrust must rely 
on a measured, or assumed, distribution of pressure in 
the region surrounding the jet at exit. (The presence of 
a post-exit thrust implies that the internal flow from the 
exit affects the pressure distribution on the external 
surfaces.) 

Some consideration must also be given to the con- 
ditions of the issuing jet at exit. When measurements 
are made at this position in flight or on a test-bed, the 


*The term “intrinsic thrust” is preferred to “internal thrust ” 

because throughout this Report “internal” means “ relating 
to the fluid that flows through the body,” which embraces the 
fluid ahead of and behind the body as well as that actually 
within the body. The word “intrinsic” is thought to be 
etymologically suitable in this context. 


+The term “ram drag” is due to Jakobssen“™) and is to be 
preferred to the term “intake drag,” which is also used to 
denote the aerodynamic drag of the actual intake. Ram 
drag is the same as the “sink drag” which occurs with 
suction aerofoils. 


magnitude and direction of the resultant velocity must 
be found, at all points of the cross section, in order to 
obtain the velocity component V.cos¢%,. The engine 
manufacturer, in preparing brochures from which 
estimates can subsequently be made of the engine thrust 
in given conditions, envisages a somewhat idealised 
picture of the jet flow at the exit and assumes that the 
gas velocity is entirely in the direction of the undisturbed 
stream. An “effective” exit area A; (at some station 
“f*) is then quoted by the manufacturer, which differs 
from the real exit area A., to compensate for this simpli- 


fication. These two viewpoints will be borne in mind ‘ 


when deriving expressions for thrust to be used in 
practice. A further discussion of the relation between 


“ 


quantities at stations “e” and “f” is given in 


Appendix III. 


4. Practical Definitions of Thrust* 
4.1. THE STANDARD THRUST 

If the post-exit thrust can be neglected, the thrust is 
given by the sum of the pre-entry and intrinsic thrusts, 
This is called the net standard thrust and is the thrust 
normally specified in engine brochures. In terms of 
the conditions at the duct exit station the net standard 
thrust is defined by 


(pV? COS? 6) 


When an effective area is given (at station “f”), the net 
standard thrust is 


| (pV? +P;—P..)dA,— mV... . 


In equation (7), cos’ ¢ does not appear because at station 
“f” the flow is assumed to be in the direction of the 
undisturbed stream. The omission of the second term 
(the ram drag) in either of the above equations gives the 
gross standard thrust (Ox). 


4.2. THE JONES THRUST 

If the jet is subsonic at exit, and exhausts into a 
region where the pressure differs from P., the simplest 
assumption that can be made for calculating the post- 
exit thrust is that there is no transfer of momentum or 
energy. during the subsequent flow to infinity where the 
pressure is P,. It is recommended that the thrust 
defined by this model of the flow, in which the fluid 
downstream of the exit is assumed to be expanded or 
compressed adiabatically and isentropically to the 
undisturbed pressure P,., should be called the net Jones 
thrust, since it represents an extension of B. M. Jones’ 
original method”) of obtaining the drag of a body in 
incompressible flow, to the problem of determining 
thrust. 

Using the condition of continuity of mass-flow in the 
wake and the energy equation for the compressible flow 
of a perfect gas, it may be shown (Appendix II) that the 


*See Appendix II for a more detatled development of the 
equations given in this Section, 
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net Jones thrust is given, in terms of conditions at the 
exit plane, by 


2 


x {1- cos dA, — mV... 


The ratio of the specific heats of the fluid at exit (y.) is 
assumed to remain constant in the subsequent flow to 
infinity, but may differ from the value of y appropriate 
to conditions ahead of the body. 

Substituting suffix “f” for suffix “e” and putting 
cos@ equal to unity in equation (8) would give the net 
Jones thrust in terms of conditions at the idealised flow 
plane “f,” where the flow is assumed to be everywhere 
parallel to the undisturbed stream. 

As before, the omission of the term mV. leads to 
the gross Jones thrust (Q5,). 


“ 


43. THE PEARSON THRUST 

When there is a difference of pressure across the 
solid boundary between the internal and external flows 
at an infinitesimal distance upstream of the exit station, 
as there may be when the jet is supersonic, a more 
realistic assumption is that the issuing fluid either 
expands or is compressed until it attains the local static 
pressure of the surrounding fluid P,, (assumed uniform), 
by an adiabatic process in which the quantity 

(pV? cos*¢+P—P,)dA. ‘ (9) 
remains constant along each streamtube. It is assumed 
that when the pressures in the internal and external 
flows have become equal, the flow is parallel to the 
direction of the undisturbed stream (station “g”) and 
there is then no transfer of momentum or energy during 
the subsequent .flow to infinity where the pressure is 

Applying the energy equation to both parts of the 
flow (between the jet exit and station “g,” and between 
station “g” and infinity downstream) and using the 
condition that the quantity given by (9) remains constant 
between the exit and station “g,.” the thrust obtained 
is (see Appendix IT) 


+ 


cos? | dA.—mV,.. 


(10a) 


Alternatively, in terms of conditions at the idealised 
flow station “f,” this can be written 


b 


It is recommended that the thrust 9), be called the 
net Pearson thrust, after the member of the Panel who 
derived the equations; the gross Pearson thrust Op, 
results from the exclusion of the ram-drag component 
mV... 


5. Comparison of the Three Practical 
Definitions of Thrust 


If the jet discharges into a region where the static 
pressure P,, is equal to the undisturbed pressure P,,, the 
post-exit thrust is zero and the standard thrust is 
appropriate for all exit velocities and pressures; the 
Pearson and standard thrusts are then identical. The 
former may therefore be regarded as a generalisation 
of the standard thrust, to be used when the jet exhausts 
into a region whose local static pressure differs from 
that of the free stream. 

Where P,, is not equal to P,,, there may be a post-exit 
thrust. For a subsonic jet, there can be no difference 
of pressure across the boundary between the internal 
and external flows at the exit, and if the flow inclination 
at the exit is small, so that cos@. ~ 1, the Pearson and 
Jones thrusts are then equal. Hence the Pearson 
thrust can be regarded as a generalisation of the Jones 
thrust to the case where a difference of static pressure 
exists between the internal and external flow at exit, as 
in a choked or supersonic jet. 

Finally, a comparison may be made of the Jones and 
standard thrusts. For this purpose the flow inclination 
at exit ¢ will again be assumed to be small, so that 
cos¢.=1. Expanding the radical in equation (8) in 
powers of (P.—P..)/p-V.” then gives* 


= [pve [1+ 


] 
dA. 
whereas, with the same assumption that cos ¢.=1, 


Thus (if ¢ is small) the standard thrust is equal to 
the first two terms only of the expanded formula for 
the Jones thrust. The difference between the two thrusts 
is of course equal to the post-exit thrust, and is shown in 
Fig. 2. 

The condition of zero momentum or energy transfer 
in the post-exit flow, assumed in deriving the Jones 
thrust, will not be obtained in practice, although with a 
subsonic jet the errors involved are unlikely to be large. 
In certain cases, however, the use of the Jones thrust is 
clearly wrong. For example, when the total pressure of 
the jet at exit is less than the static pressure of the 
undisturbed stream, the fluid flowing through the body 
cannot attain the pressure of the undisturbed stream 
again, unless there is some transfer of momentum or 
energyt Alternatively, if the velocity at exit is sonic 


*The subsequent equations could equally well be written in 
terms of conditions at station “ f.” 


+This type of flow may occur for example in the flow through 
a piston-engine cooling duct as the exit flap is opened (see 
also Appendix V), 
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Figure 2. Curves of the difference between the Jones thrust 
and the standard thrust. 


or supersonic, there can be a difference of static pressure 
at exit between the internal and external flows. With an 
ideal fluid and no mixing, the wake could exhibit a 
periodic flow pattern extending to infinity and there 
would be no station at which the internal static pressure 
was at once uniform and equal to that of the 
undisturbed stream. 

The latter case, which is the more important in 
thrust problems, can be dealt with by using the Pearson 
thrust. The assumption made in deriving equations 
10 (a) and (b) probably represent the best simple 
model of the flow that it is possible to obtain, though 
this model may not always correspond with reality. If, 
after the initial non-isentropic process, there is sub- 
stantial mixing of the internal and external flows before 
the jet has reached the pressure P,, the actual thrust 
obtained will be affected and the assumption made by 
Pearson violated. In addition the external flow itself 
may not be isentropic; for example shock waves may be 
present. Moreover the Pearson thrust depends on the 
value of P,; in practice this is unlikely to be uniform 
and some mean pressure must be used in the calcula- 
tions. It is felt, however, that points such as these do 
not seriously affect the usefulness of the flow model and 
resulting thrust given by Pearson. 


6. Internal Drag and Propulsive Thrust 

In problems of internal drag, similar concepts to the 
three types of thrust given in Section 4 are required. It 
is therefore recommended that the quantities defined by 
the negative of the three types of net thrust should be 
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called the standard internal drag, the Jones internai drag 
and the Pearson internal drag, if such quantities are 
required. The adjective net is unnecessary because the 
concept of gross internal drag has no significance. Some 
formulae for internal drag are given in Appendix V. 
The overall longitudinal force on a nacelle housing 
an engine, treated as a whole, is felt to be worthy of a 
name and it is suggested that this be called the propul- 
sive thrust, as it is the thrust available for propelling the 


rest of the aircraft. 


7. Concluding Remarks and 
Recommendations 

The foregoing investigation has shown that the thrust 
can be divided into pre-entry thrust, intrinsic thrust and 
post-exit thrust. In the absence of the last component, 
the standard thrust is obtained. Two idealised models 
of the flow downstream from the exit were put forward 
and enabled an estimate to be made of the post-exit 
thrust in certain cases. This led to the concepts of Jones 
thrust and Pearson thrust, which can be used to assess 
the amount of overall longitudinal force that it is fair to 
attribute to the internal flow. The flow models are not 
claimed to be exact representations of the actual flow 
but it is felt that they are sufficiently realistic to be worth ! 
serious consideration. At this stage it is not possible to _ 
estimate with any accuracy their possible usefulness. 
There are related internal drag forces corresponding to | 
each of the three types of thrust. 

It is suggested that, of the thrusts considered 
here, the Pearson thrust represents the best measure of 
the thrust that it is fair to ascribe to a jet engine in the 
general case. When the internal static pressure at the exit 
station is equal to the local static pressure surrounding 
the jet exit (subsonic, or unchoked, jet), the Pearson 
thrust reduces to the Jones thrust. When the local <* suc 
pressure surrounding the jet exit is equal to the 
undisturbed static pressure, the Pearson thrust reduces p 
to the standard thrust; in this case the post-exit thrust 
is zero and the standard thrust is the correct measure ; 


of the thrust. 
For ease of reference the various thrusts are collated 


and briefly defined in Appendix I. 


APPENDIX I 


DEFINITIONS OF SOME QUANTITIES USED IN THIS REPORT | 
1. Internal Flow.—The flow of fluid which passes 


through the duct. 
2. External Flow.—The flow of fluid which passes — 


around the outside of the duct. 
3. Pre-entry Stream Tube.—The stream tube, extend- 


ing from infinity upstream of the body to the duct entry, | 

whose boundaries separate the internal and external flows. | 
4. Equivalent Post-exit Stream Tube.—The stream 

tube, extending from the duct exit to infinity downstream, — 


whose boundaries are assumed to separate the internal and _ 


external flows in the absence of mixing. 

5. Ram Drag.—The force, acting in the downstream 
direction which would arise if the momentum of the 
internal flow at infinity ahead of the body were destroyed. 
(It is equal to the product of the rate of mass flow at the | 
duct entry and the velocity of the undisturbed relative | 


airstream, and is the effective sink drag of the air intake.) | 
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JET-ENGINE THRUST AND DRAG DEFINITIONS 


6. Thrust (associated with the internal flow).—The 
force corresponding to the rate of increase of momentum, 
in a direction parallel to the undisturbed relative airstream, 
of the internal flow through the duct, this increase being 
calculated between stations at infinite distances upstream 
and downstream of the body. The thrust can be divided 
conveniently into three parts, the pre-entry thrust, the 
intrinsic thrust, and the post-exit thrust. 

7. Pre-entry Thrust.—The force, in a direction parallel 
to the undisturbed relative airstream, arising from the 
pressure forces acting on the internal surface of the pre- 
entry stream tube. 

8. Intrinsic Thrust—The force, in a direction parallel 
to the undisturbed relative airstream, arising from the 
pressure and friction forces acting on the internal surface 
of the duct. 

9. Post-exit Thrust——The force, in a direction parallel 
to the undisturbed relative airstream, arising from the 
pressure forces acting on the internal surface of the equi- 
valent post-exit stream tube. 

10. Standard Thrust.—The net standard thrust is the 
sum of the pre-entry thrust and the intrinsic thrust. The 
gross standard thrust is the arithmetic sum of the net 
standard thrust and the ram drag. 

11. Jones Thrust——The net Jones thrust is the sum 
of the pre-entry thrust, the intrinsic thrust and the post- 
exit thrust, the last being calculated by assuming that the 
fluid that has passed through the duct is compressed or 
expanded adiabatically and isentropically downstream of 
the duct exit until the undisturbed static pressure of the 
relative airstream is attained. The gross Jones thrust is the 
arithmetic sum of the net Jones thrust and the ram drag. 

12. Pearson Thrust.—The net Pearson thrust is the 
sum of the pre-entry thrust, the intrinsic thrust and the 
post-exit thrust, the last being calculated by assuming that 
the flow downstream of the duct exit has two stages. In 
the first, the fluid is compressed or expanded adiabatically 
and irreversibly until the local static pressure of the sur- 
rounding fluid (assumed uniform) is attained, the quantity 
(pV? cos?¢+P-—P,)dA remaining constant along each 
streamline; the flow is then assumed to be parallel to the 
direction of the undisturbed relative airstream. In the 
second stage the fluid is compressed or expanded adiabati- 
cally and isentropically until the undisturbed static pressure 
of the relative airstream is attained. The gross Pearson 
thrust is the arithmetic sum of the net Pearson thrust and 
the ram drag. 

13. Propulsive Thrust—The overall force exerted on 
a nacelle housing an engine, in a direction parallel to the 
undisturbed relative airstream; it is considered positive 
when acting in the upstream direction. 

14. Standard Internal Drag.—The negative of the net 
standard thrust. 

15. Jones Internal Drag.—The negative of the net 
Jones thrust. 

16. Pearson Internal Drag.—The negative of the net 
Pearson thrust. 


APPENDIX IL 

DERIVATION OF THE THRUST AND RELATED QUANTITIES 

Figure 1 represents the “internal” flow ahead of, 
through and to the rear of a stationary ducted body. 
Consider the control volume bounded by the intake and 
exit stations “i” and “e” (taken for convenience as 
normal to the free-stream direction) and the internal sur- 
faces of the body. Momentum considerations show that 
the sum of the pressure and friction forces acting on the 


fluid at the boundaries is equal to the rate of change of 
momentum of the fluid flowing through the control volume. 
Resolving in the direction of the free stream gives 


| (P F eos 9) ds 


i e int 


= | cos” Pe dA, pV? cos” pi dA,, 


e 1 


where F is the local friction force per unit area, dS is an 
element of area of the internal surface, and @ is considered 
to be positive when the duct is diverging. (The third 
integral on the left-hand side represents the force in the 
free-stream direction exerted on the fluid by the internal 
surfaces of the duct.) 


Hence | (P sin — F cos @) dS = | (pV? cos? @, + P.) dA, 


int e 


five cos” dA,. 


The integration of a constant pressure P,, (the undisturbed 
Static pressure) round the closed boundary must give zero 
force, 


ie. | P dA, - | P.dA,+ Py sin dS =0. 
1 e int 


Hence by subtraction 


| | (P — P.) sin F cos dS = 


Int 


| (p.Vez C08? +P, —P..) dA, - 


e 


- cos* +P,;- dA,. 


The forces that are of significance are those associated 
with pressures measured relative to the undisturbed static 
pressure. Thus the left-hand side of the above equation 
gives the force in the free-stream direction exerted on the 
fluid by the internal surfaces of the duct, measured 
positively in the downstream direction. This is clearly 
equal to the effective force in the direction of the free- 
stream exerted by the fluid on the internal surfaces, 
measured positively in the upstream (thrust) direction, and 
this quantity will be defined as the “ Intrinsic Thrust,” 90 


Thus 
= | cos? +P. - P.) dA, 


e 


int’ 


| (pV? cos? +P,- P.) dA,. (11) 
i 
Now consider the flow between stations at infinity 


upstream (co) and the exit (e); in similar manner a 
“ thrust ” can be derived which is given by 


| (p.V.2 cos? +P, Py) dA,— (pV Py) dAQ, 
e ~ 


(12) 


and is equal to the sum of the longitudinal forces on the 
boundaries of the pre-entry stream tube and the internal 
surfaces. 
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By subtraction, then : — 
(* Thrust ” given by eqn. (12)) -©j,,—axial force on pre- 
entry stream tube 


| (p,V.2 cos? 9, +P, - P..) dA,— 


i 
This will be defined as the “ Pre-entry Thrust,” Oye. 
Thus 
cos” +P, dA,— mV.» (13) 
i 
since V. is constant over the whole area A, and 
M,= py VA 


a" 

Similarly, if the flow is considered between stations at 
the exit and at infinity downstream (w), a thrust is intro- 
duced which is defined as the “ Post-exit Thrust’ and 
given by 


w 


pre — 


since P,. =P, and @,,=0, although V,, is not uniform. 
Since both the pre-entry and post-exit thrusts are com- 
municated to the body, the total thrust is given by 


5 | dAy MV 


The net standard thrust is defined as the sum of the 
pre-entry and intrinsic thrusts and is given by 


or a similar equation in terms of station * f.” 

In the Jones model of the flow downstream from station 
“e” or “f,” the expansion or compression from the local 
pressure P,, to the pressure of the undisturbed stream P, 
at station “w” is adiabatic and isentropic, no mixing 
occurring between the external and internal flows. 

The energy equation for the compressible flow of a 
perfect gas can be written 


P. 
Ye- lp, Ye- 1 Pw 


Using the condition of isentropic flow, this equation 
becomes 


Substituting this value of V,. into equation (15) and using 
the condition of continuity, 


PwV =pV, cos 


gives the net Jones thrust 


1 +- G.-)M2 1 - cos o.dA, 


-mV,. (18a) 
If the station “f” is used, @, is zero and 


2 
V 


(185) 


The Pearson model of the flow has for its initial stage 
a process in which the fluid flowing internally expands or 
contracts adiabatically and irreversibly from the internal 
static pressure at the exit to the local external pressure P,, 
the process being complete at station “ g” where the flow 
is assumed to be parallel to the undisturbed stream. During 
this process, the pressure at the dividing surface between 
the internal and external flows is assumed to have the 
uniform value P,, and for each stream tube the quantity 
(pV? cos? ¢+P-—P,)dA is assumed to remain constant 
between the exit and station “ g.” 

The gross Pearson thrust ©), is given by the general 
equation 


@= pyV y2dA wy 


with conditions at station ‘* w”’ related to those at the exit 
by means of the flow model suggested by Pearson. By 


continuity, 
PwV pd Ay = pV, Cos 


Thus | PV Vy COS GdA,. . (19) 


For the isentropic and adiabatic flow between stations * g” 
and “ w,” the energy equation for a compressible fluid can 
be used, as in considering the Jones thrust, to give 


(It is assumed that y does not change downstream of the 
exit.) Between stations ““e” and “ g,” 

pg "aA, cos? P,— Py) dA,=9. 
But by continuity 
pV dA,=pV, cos A,, 
pV? cos* P, Py 


pV. COS 


and hence Vs 


It remains now to eliminate M, from equation (20). 
This can be done by applying the energy equation to the 
adiabatic flow between “e” and “ g,” giving the following 
equation 


Substituting (21) and (22) in (20) and simplifying the result, 
leads to the following expression for the net Pearson thrust 
when the ram drag term is included 


eV. cos? P,—P, |? 


Ye 

(234) 


In terms of conditions at station “f,” where the flow 
is assumed to be parallel to the direction of the free} 
stream (i.e. ¢;=0) the net Pearson thrust is 


VP +P, — YE 
| pV? P, 


P 2 $ 


(236), 
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APPENDIX III 


SOME REMARKS ON REAL JET ENGINES AND THE USE OF 
STATIONS “e” AND “f” 


In general, the discharge coefficients of engine exhaust 
nozzles are less than unity. The engine manufacturers 
assume that the engine behaves as if the flow leaving the 
final nozzle converges and passes through some area 
(known as the effective area and denoted by A; in the 
present Report) in which the streamlines are all parallel 
to the undisturbed stream. This area is calculated from 
4 knowledge of the thrust, mass-flow and exhaust temper- 
ature, measured on a test-bed, or better still in an altitude 
test-cell. For estimating flight performance, either the 
effective area is assumed to be independent of altitude, or 
an assumed variation is used, or the results obtained from 
an altitude test-cell are employed. In the present context 
this variation and the method of obtaining it are 
immaterial; it is only necessary to allow for the existence 
of an effective area when constructing the expressions for 
the various thrusts, whether this area remains constant or 
not. 


In the case of a real nozzle from which the flow 
contracts and passes through an “effective area,” the 
contracting jet (which is assumed not to mix with the 
surrounding flow) can support a pressure from the external 
flow and thus contribute to the total thrust. The problem 
is unaltered if the boundary of the jet is replaced by a solid 
surface, so that for the purpose of defining the thrust an 
idealised picture of the jet may be used, in which the new 
exit area, at which the flow is parallel to the free-stream, 
is the effective exit area A,;. As was said earlier, the exit 
area given by the engine maker is this effective area A;; 
expressions for the various thrusts in terms of conditions 
at station “ f ” are therefore required. 


The important case of flight measurements of thrust 
must now be considered. These measurements are made 
in the jet, near the engine outlet, and a method is required 
for calculating the thrust from them. Now it may be 
readily shown that it is impossible to separate the internal 
thrust from the external drag unless the wake surveys are 
carried out in the plane of the final nozzle (station “e”’), 


and a complete picture obtained of the flow at this position. 


It is therefore of great importance to have alternative 
expressions for the thrust in terms of quantities which can 
be determined in the actual exit plane. 


One method that can be adopted is to use the results 
of the measurements at station “e” to calculate the 
conditions at station “f,” and this method will be used 
where a knowledge of the effective area is required. But 
in the case that is probably of the greatest importance, a 
choked jet exhausting into a region whose pressure is equal 
to the undisturbed pressure, a simpler approach can be 
used.* In the absence of any mixing with the external 
flow the pressure and velocity must be constant along the 
bounding stream tube of the jet, the pressure being equal 


*The external pressure surrounding a jet exit under static 
conditions is not quite equal to the pressure of the undisturbed 
Stream. This is because the jet induces an external flow around 
the nacelle which changes the external local pressure, parti- 
cularly at the jet exit. The thrust obtained from engine 
brochures is ultimately based on test-bed measurements 


_ obtained under such conditions, and thus includes any post- 


exit thrust due to P, not being equal to P.,. The importance 
of this effect is not known. 


to P,.. Between stations “e” and “f,” the following 
relation holds 


(p.V.2 cos? P,.) dA, - (pV dA,—0. 


The difference between the two integrals given above is 
equal to the resolved component of the force due to the 
pressure difference (P — P,.) acting on the boundary of the 
jet between the two stations. This difference must be 
zero, since on the boundary stream tube P is equal to P,.. 
Thus the thrust, given in this case by the standard thrust, 
is obtainable from the first integral in the equation entirely 
in terms of quantities in the exit plane, where measurements 
of static and total pressure, temperature and angle of flow 
are required. 

If the external pressure is not equal to the free-stream 
value, the relation given above must be modified. Assuming 
that the external pressure has the constant value P,, between 
stations “e” and “ f,” 


| 00s" 6.4 P,—P,) dA. = | Pi) dy 


so that the velocity V, in the intermediate plane “ g” can 
be related to conditions at station “ e.” 

Thus, for the purposes of the definition and calculation 
of thrust, when conditions at station “e” are known, it 
is often not necessary to calculate the conditions at station 
“f." The choice between stations e” and “ f” is mainly 
one of convenience. In general, in the absence of measure- 
ments, it is almost impossible to estimate conditions at the 
exit plane, but it is relatively easy for the engine manu- 
facturer to estimate the conditions at station “ f.” 

In many cases, too, the velocity and pressure can be 
assumed to be uniform across station * f ” when estimating 
thrust from data supplied by the manufacturer, and there 
is then no need to integrate over the jet area. To define 
the thrusts in terms of uniform flow at station “ f” would 
be unnecessarily restrictive, however, since the exact 
correspondence between points in stations “e” and “f” 
connected by a common stream tube would be lost. 


APPENDIX IV 


SOME CAUTIONARY REMARKS REGARDING CENTRAL BODIES AND 
SIDE INTAKES 


All the thrusts (or internal drags) defined in Section 4 
are based on the rate of change of momentum of the 
fluid that flows through the body between a station at 
infinity upstream and some downstream station. If, there- 
fore, the ducted body has a central body protruding ahead 
of or behind the main body (see Fig. 3), the drag of the 
portion of the central body protruding ahead of the entry 
is automatically included in the internal thrust or drag, 
just as are portions entirely within the main body. The 
drag of any portion of the central body protruding behind 
the exit is not necessarily included, however; the method 
to be used in allowing for this drag depends on the 
definition of thrust that is adopted. 

When the central body protrudes a long way ahead of 
the entry, it is often more convenient to treat the central 
body as the main body, and to treat the ducted part as a 
“side intake.” This is particularly necessary when the 
intake does not entirely surround the central body. In 
such cases the internal drag of the duct, under both 
definitions, will include the drag of the portions of the 
main body ahead of the entry that are wetted by the fluid 
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Surfaces of the body which contribute to 
the internal drag are shown by heavy lines. 


tit 
Diagrammatic representation of a ducted body 
having a central body within it. 


FiGure 3. 


that flows through the duct. However, it is usual to include 
the drag of the portions of the main body wetted by the 
internal flow in the external drag, and to correct for this 
by reducing the pre-entry drag of the intake by the same 
amount, to yield, say, a “ modified pre-entry drag.” That 
such a procedure is rational is apparent when the ducted 
parts of the body become so detached from the main body 
that they can be treated as separate nacelles. 


APPENDIX V 
FORMULAE FOR INTERNAL DRAG 

One of the earliest attempts to define the internal drag 
of a ducted body was made by Meredith,'*’ but he con- 
sidered only the case when the exit static pressure was equal 
to the undisturbed pressure. In this case all the definitions 
compared in Section 5 are equivalent to one another and 
to Meredith’s definition. 

Other definitions (unpublished, 1936) covered the case 
in which the exit static pressure was different from the 
undisturbed pressure and were equivalent to equation (1) 
of this report, but in order to apply the definition the 
velocity V,, was related to the exit velocity V, and the exit 
static pressure P, by the assumption of no transfer of 
momentum or energy in the wake. This method is now 
commonly called Jones’s method, because Jones made the 
same assumptions in obtaining the profile drag of an 
unducted body from measurements in the wake.) 
Taylor?) has pointed out, however, that the effect of 
turbulent mixing is not always negligible, particularly when 
the exit static pressure is appreciably different from the 
undisturbed pressure. Jones’s formula, for constant mass 
flow through the duct, neglecting is 


After it had been observed that the definition of internal 
drag given by equation (1) could be expressed in terms of 
measurable quantities by the use of what was, in effect, 
Jones’s formula, it was natural that Betz’s formula‘*) should 
also have been considered, for it had been shown that the 
two formulae were equivalent to the first order when used 
to estimate the profile drag of unducted bodies.’ Betz’s 
formula, which applies to incompressible flow, is 


and on expanding the radical in powers of (P, — P,.)/4 pV’ 
this becomes 


ov.[v.- {1 + = 
1 Pe)? 


By expanding the radical in equation (24) in powers of 
(P, — P,.)/+ pV.”, the equivalent form of Jones’s formula 
is found to be 


D,= | pV. V, {14+ 


Thus y eR ws and (27) differ only in the terms of 
and higher. 


order 


Betz’s formula has the apparent advantage over Jones’s } 


that it gives a real answer when the total pressure at exit 
is less than the undisturbed static pressure, for which case 
Jones’s formula gives a complex answer. Betz’s formula 
has, in fact, been introduced into internal drag problems 
(e.g. in Ref. 10) in order to cope with this “ weakness ” of 
Jones’s formula, although Squire in 1936 pointed out that 
Betz’s formula should not be used when Jones’s would give 
a complex answer. Betz’s formula was in fact developed 
for determining the profile drag of closed (unducted) bodies, 
and has never been extended to give a formula for the 
internal drag of a ducted body. 


It is recommended therefore that Jones’s formula should 
be used if an assessment of the internal drag of a ducted 
body is required, and that Betz’s formula should not be 
used for this purpose. 
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~«| The Stress Analysis of the Circular Conical 
mula 
Fuselage with Flexible Frames 
by 
(27) 
W. J GOODEY, M.A., A.F.R.Ae.S. 
“< Summary: This paper is concerned with the problem of the skin-and-stringer-covered conical 
fuselage of circular cross section, having flexible frames and being subjected to an arbitrary 
distribution of loads applied in the planes of these frames or on the end cross sections. The 
; skin thickness and stringer cross section may vary from bay to bay but are assumed to be 
ness constant in a bay, i.e. between any two consecutive frames. The stringers are assumed to be 
t exit uniformly spaced round the circumference. The frames are assumed to be of uniform cross 
| Case section (circumferentially) but their stiffnesses are arbitrary, as is also their spacing along the 
‘mula fuselage. 
blems This paper is a sequel to one previously published in the JOURNAL), but the method of 
5” of analysis is now by minimum strain energy instead of by deflections. Owing to the generality 
t that of the investigation it is impossible to present the solution—in terms of the known loads and 
i oi stiffnesses—in any explicit form, and attention has therefore been concentrated on the technique 
give of numerical computation, A scheme of tabulation has been evolved, and its use is illustrated 
— by a fully-worked numerical example of a moderately complex type. 
odies, 
sie. 1. Introduction (iii) A part of the fuselage consisting of a frame and 
the structure between this frame and the next 
ome (vay from the ofthe con) wil be 
ot be sectional dimensions small ared referred “Dey.” covering: 
th th length, is liable to be very different that 
P cumference. Thus, the stringers lie along 
of Jet title generating lines of the cone. The skin thickness 
erence Ne ee and stringer cross section may vary from bay 
Memo. } 4 Circular Fuselage ”*. to bay 
{2 In spite of its title it must now be admitted that this (iv) The stringers will be assumed to be sufficiently 
paper did little more than demonstrate the existence of closely spaced to be regarded as a continuously 
oe the si — it may have served a useful distributed reinforcing of uniform thickness. 
, ie purpose if it did only this. Bending and torsional stiffness of the individual 
149. The present investigation is of a much more general stringers (in association with a piece of skin) 
character, and is concerned with the problem of the skin- will be neglected. 
oni and-stringer-covered conical fuselage of circular cross (v) As already mentioned, the applied loads are 
elias section, having flexible frames and being subjected to an either in the planes of the frames, or on the end 
arbitrary distribution of loads applied in the planes of cross sections, but are otherwise arbitrary. The 
BS these frames or on the end cross sections. A typical analysis could in fact be made to cover prac- 
problem of this type is illustrated diagrammatically in tically any type of loading, but it is convenient 
anil Fig. 1. } ; ; to impose this restriction, which involves only 
eg The basic assumptions governing the analysis are a slight loss in generality. 
as follows. Any further assumptions of a secondary nature 
‘ircraft (i) The surface of the fuselage forms part of a right which it may be desirable to make will be brought in as 
Glycol circular cone, the frames being in planes the occasion arises. 
anuary normal to the axis of the cone. 
(ii) The cross section of a frame is uniform, but 
by the the cross-sectional dimensions, and therefore 
roan the stiffness, may vary from frame to frame. 
' The stiffness which is of primary concern is the 
sect resistance of a frame to distortion in its own 
i plane; any resistance to bending or twisting out 
of the plane will be neglected. The spacing of 
Fluid : : 
i? the frames along the axis of the fuselage is 
arbitrary. 
nique | Originally received November 1953. 
ests. in *See also Ref. 2 which includes an interesting review of current 
1947. literature on the subject. Ficure 1. 
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total cross section of end-load carrying 
material 

in Section 10, symmetrical matrix of struc- 
tural coefficients 

modulus of rigidity 

depth between the centroids of the two 
flanges of frame j 

effective depth of web of frame j 
Young’s modulus 

defined in equation (21) 

defined in equation (21) 


offset of the neutral axis of the cross 
section of frame j from the skin line 


moment of inertia of cross section of 
frame j, for bending in its own plane 
defined in equation (20) 

bay j consisting of frame j and skin and 
stringers between frames j and j+1 
defined in equation (20) 

length of bay j=2;,, Z; 

in Section 10, column of loading terms 
term in strain energy equation nj due to 
the applied loading 

bending moment in a frame, at the cross 
section defined by angle 4 (see Fig. 3) 
value of M when applied loads are 
balanced by skin loading of elementary 
bending and torsion theory 


externally applied moment loading, posi- 
tive anti-clockwise looking along Oz away 
from O, on frame j (see Fig. 3) 

arbitrary constants defined in equation (29) 
direct stress loading (i.e. tension per unit 
width) along the generating lines, taken by 
the skin and stringers 


externally applied radial loading, positive 
outwards, on frame j 

defined in equation (18) 

defined in equation (19) 

defined in equation (18) 

defined in equation (19) 


shear loading (shear per unit length) in the 
skin at section z, parallel and normal to 
the generating lines, positive as in Fig. 2 


shear loading applied to frame j by the 
skin, positive when acting in the direction 
of increasing 6, i.e. anti-clockwise looking 
along z-axis away from origin (see Fig. 3) 
externally applied tangential loading (at 
skin radius Rj) on frame j 

radius to the skin surface at section z 
value of R at frame j 

shear in a frame at the cross section 
defined by angle 6 (see Fig. 3) 

value of S corresponding to M 

end load in a frame at the cross section 
defined by angle 6 (see Fig. 3) 

value of T corresponding to M 


in Section 10, lower triangular matrix and 
its transpose 


t; skin thickness in bay j 
ty; effective thickness of end load carrying 
material in bay j 
tw; web thickness of frame j 
U strain energy of complete structure 
U; strain energy of bay j 
U,; total strain energy of frame j 
Uy,r; strain energy of bending of frame j 
Uy; strain energy due to loading N in bay j 
Usp; shear strain energy in frame j 
Us; shear strain energy in the skin of bay j 
U;,p; strain energy of end load in frame j 

v tangential deflection of the point on the 
neutral axis of the frame defined by angle 
6, taken positive in direction of increas- 
ing 0 

W load applied to a frame, in Section (8) 

w radial deflection of the point on the 
neutral axis of the frame defined by angle 
#, taken positive outwards 

x distance of section z from the frame just 
behind it (i.e. looking towards O). For 
bay j, 

Y column of subsidiary variables 

z distance of any cross section normal to the 
axis of the cone from its apex O 

a semi-angle of the cone=tan~' R/z 

y; Yatio of total flange area of frame j to 
effective web area=/)dy;/ 

6 angle between plane zOA and a fixed 
plane through zO (positive anti-clockwise, 
measuring from the “ downward vertical,” 
and looking in the direction of z 
increasing) (see Fig. 3) 

Aj =(Rj.. Ri) /(Ri.. +R) 

p; Yadius of gyration of the cross section of 
frame j 

7; torque about the z-axis, in bay j 


defined in equation (7) 


VYni 


2. The Equations of Equilibrium 


The position of any point A on the skin surface may 
be specified by the co-ordinate z and the angle 6 between 
the plane zOA and a fixed plane through zO. This fixed 
plane can be conveniently regarded as the “ vertical” 
plane. 

Referring to Fig. 2, consider now the equilibrium of 
the element ABCD of the shell bounded by sections z 
and z+6z and the generating lines defined by the angles 
6 and 6+66. 

Resolving along the normal to the element there is 
no circumferential stress in the skin, since g and WN are 
both in the tangent plane. Bending stiffness of the 
stringers (+skin), which would make possible the exis- 
tence of shear forces normal to ABCD, is neglected, 
while normal pressure on the skin is assumed to be zero. 

Resolving along the generating line AB, 


a 
az NR)+ sp 


and taking moments about O, 
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(The same equation comes from resolving along AD, 
but the derivation is less elegant.) 
Thus, qgRz is independent of z or, since z/R is 
constant, 
is independent of z_. (3) 


This only applies between one frame and the next, 
since there may be tangential loads applied to the skin 
at a frame, so equation (2) is not valid if the element 
ABCD includes a portion of a frame. 

Hence qgR? is a function of the angle @ and of the 
bay number j. 

Equation (1) may now be written in the form 


ny 1 oO 
(NR cos 2)= — R (GR), 


which can then be integrated with respect to z. How- 
ever, to make the formulae applicable to the limiting 
case of a circular cylinder (when the point O moves 
away to infinity and R becomes constant) it is convenient 
to replace oz in this last equation by ox. It is then in- 
tegrated with respect to x, adjusting the “constant” of 
integration (a function of 6) to make N=N, at frame j 
(x=0). The result is that 


ta 
NR cos z=N;,R;cos (4) 
0 
z R 
1 dR\R, R RR;’ 
j 
since 
R-R; 
Therefore 


x 
NR cos 2=N;Rj; cos 2 - RR, (aR ) (5) 


\ >+—-—! 
B 
\ 


(“verticat”) 


(Note. When z2=0, R=constant, so that N then varies 
linearly with x between frames.) 
When x=L;, N=N;,,, and so from (5) 


= (W,R,c08 2 — C08 2) 
(6) 
Given N; and N;,, this equation enables gR? (for bay 
j) to be calculated, by integration with respect to 6, the 
constant of integration being adjusted to give the correct 
torque about Oz. 


_ 3. Types of Stress Distribution 
Let N,R;cos z= Yv,;cosn6+ Xo,;sinnd . (7) 
n=0 n=1 


6 is taken to lie in the range - = to +7. The upper 
limit of 7 is omitted from the summation signs because, 
although in theory it will generally be infinite, in practice 
only enough terms are retained to give the required 
degree of accuracy. The bending moment about the 
axis XX (see Fig. 2), at frame j, is equal to 

+z 

| (N;R; cos z) Rj sin@d6==Rjo, ;, 
positive moment giving tension on the side 6=0 to 7. 
Similarly, the bending moment about YY 


= W,R,cos 2) R, cos 


positive moment giving tension in the range 
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and the end load (in the z-direction) 
= Jur, cos 2)d#=2zv, ; (tension positive). 


The cosine terms (including y, ;) in equation (7) re- 
present a loading (N;) which is symmetrical with respect 
to the plane XX, the terms with odd multiples of 4 being 
anti-symmetrical with respect to the plane YY (since 
sinn(@-=/2)=+cosné when n is odd). 

The sine terms give a loading N; which is anti- 
symmetrical with respect to the plane XX, the terms 
with odd multiples of 6 being symmetrical with respect 
to the plane YY (since cosn(6- =/2)= +sinn@ when n 
is odd). The terms (v, ;+v,,;cos4+¢,,;sin 4) in equa- 
tion (7) clearly represent stress distributions of the 
“elementary bending theory” type, the first being a 
uniformly distributed end load, the second a bending 
moment about YY, and the third a bending moment 
about XX. All the other terms are self-equilibrating 
distributions, and can be regarded as corrections to the 
elementary theory. These self-equilibrating distribu- 
tions are essentially of the nature of redundancies, and 
the coefficients ;...Unj...3 Which deter- 
mine the amount of each distribution present, may be 
determined by the principle of minimum strain energy. 
In practice, calculations will be limited to finding only 
the first few terms in the series, which is usually fairly 
rapidly convergent. 

But first a formula must be derived for the skin 
shear loading in terms of the v’s and ¢’s, using equations 
(6) and (7). From these, 


xg (R*)= Yn, 1-1) COS + 
+ On, 4-1) Sin 
whence, on integration with respect to 4, 
4 
n=1 
cosn4 T 


where 7; is independent of 6, being equal to { qR*dé, and 
is therefore the torque about zO. 

Note that the coefficient (v,, ;- v,, ;.,) is omitted from 
equation (8). This coefficient would give a term 


Yo, 141) 9, 

which would make q discontinuous at 6= +7. Now by 
implication the end load at frame j is continuous (in the 
z-direction), since otherwise the quantity v, ; would be 
two-valued, ¥, ;+, say, on one side of the frame and 
v,,;- on the other side. This possibility is excluded by 
assumption (v) of Section 1. Thus, either 
=0, or there must have been a change of end load 
between frames j and j+1. This change could only be 


brought about by the application of load along a finite 
number of generating lines of the surface, since the 
equations of equilibrium imply that there is no distri. | 
buted loading of the skin-friction type. But loading ) 
along generating lines has also been excluded by 
assumption (v) of Section 1, and hence ;=¥,, 

Equations (7) and (8) give the most general form of 
stress distribution in the skin and stringers on the 
assumptions made in Section 1, and the strain energy of 
the various parts of the structure may now be calcu. | 
lated. For this purpose, and for the calculation of the ? 
strain energy derivatives with respect to the unknowns 
Vn) and >,; (n=2, 3, etc.), it has been found very con- 
venient to use matrix notation. Only the most elemen- 
tary parts of the theory of matrices are required, up to 
the differentiation of quadratic forms. 


4. The Strain Energy of Shear in the Skin 

(It is assumed that the skin, stringers and frames are 
made of the same material throughout. If not, the 
appropriate modifications to the formulae must be made, 
though probably the simplest way of dealing with this is 
to use equivalent thicknesses and cross-sectional areas 
all in the same material.) 


Then 6EU 


1,R sec zd0dx 


_ eve | (qR*)? dé dx 
~ 


0 -z 


The factor 6E is aed in for convenience in numerical 


2R?? 


L; 
1 
=cot z jar _ _ 
0 R, 
Therefore 


3EL;(R)+R;.,) sec 2)2 


6EU s;= 


which using (8) 


3E=(R;+R;,,) sec z 


2C1,L; [32 n - Vn, i4 iv 


+ | +a multiple of 


Cr = 85 mee 


The n=1 terms and 7; are known from the applied 
loads, and so they will disappear on differentiation 
with respect to any of the unknowns. Thus, only the un- 
known, or variable, part of 6EUs; need be considered. 
Moreover, expressions need only be obtained for the 
partial differential coefficients of 6EU,; with respect to | | 
the v’s, since those with respect to the ’s will be obtain- 
able therefrom, simply by replacing each v by the corres- 
ponding 9. 
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L 
Now 
ty; R sec dx 
n,i+1 
L; 
3Ex(R)+R;,,) sec z 
Jove cos 2)? dé dx. 
0 


+a similar expression with ’s instead of ¥’s 
+a multiple of 


Differentiating with respect to y,,; and v,,;., in turn, 
and denoting the pair of differential coefficients—for a 
particular value of n—by 6E0Us;/0¥,, 


Ws; 3Ex(R)+Rj,,) sec z [ ] 
n,j+1 


1 


1 1 
n n Vn, 


(n=2, 3, and so on) 


These differential coefficients could have been written 
down easily enough without using matrix notation, 
since the expression for the strain energy in the skin is 
a particularly simple one. The advantage of using this 
notation will, however, become more evident in the more 
complicated expressions for the strain energy of the 
stringers and the frames. Note that the top row of (9) 
is the differential coefficient with respect to v,,;, while the 
bottom row is the differential coefficient with respect to 
Vvija. Note also that the form of (9) is independent of 
the applied loading, and for this reason the strain energy 
terms—i.e. the coefficients of the wW’s—will be called 
“structural terms.” Terms whose form depends on the 
applied loading will be called “loading terms.” 


5. The Strain Energy of the Stringers 


Let Uy; be the strain energy due to the loading N, 
in bay j. This loading is along the generating lines of 
the cylindrical shell, and is carried as end load in the 
stringers. Each stringer may usually be assumed to have 
a piece of the skin, to which it is attached, acting with it 
as an end-load carrying member, thus giving an effective 
stringer area somewhat greater than the area of the 
stringer alone. Assuming these effective stringer areas 
to be spread uniformly round the circumference, let ty; 
be the effective thickness of end-load carrying material 
so distributed. 

Since, in practice, the number of stringers and the 
effective cross-sectional area of each stringer will be 


_ fixed quantities in a bay, the value of ty; will depend on 


R. Thus, if A; denote the total cross section of end- 
load carrying material in bay j, then 


A; 


tyj= (10) 


Now from equations (5), (7) and (8) 


NR cos (1 - (vn; cos nO + Sin nO) + 
j 


Hence 6EUy;= 
(12) 
where 


L 


i 
x\? R? X)RRin 
a= — ax, Fac dx, 


The values of these integrals may be found by ex- 
pressing the integrands in terms of R, and integrating 
between the limits R; and R;,,. Omitting the inter- 
mediate steps and merely noting that 


dx = L;dR/(Rj,, Ri) 


and x/L;=(R—R))/(Ri.. — Ri), 

the are results are obtained: 
a; = +R; lo Ris, 

and 
1 R; 2R; lo R;. 
L; ~ (Ri - RY Ri. — R; R, 


These expressions are not very convenient for prac- 
tical use, particularly if (R;.,- R;) is small compared 
with Rj, since they take the indeterminate form 0/0 
when R;,,=R;. In order to avoid this difficulty, and to 
obtain approximate values for a;, b; and c; which may be 
used when—as will usually be the case—R;., differs only 
slightly from R;, we write 


Ri, RB; 
(13) 
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The logarithmic term may then be expanded as a 
power series in A;, and it is found that 


1 


In practice the square and higher powers of A; may 
be ignored, and it will be sufficient to take 


L; 3 3(R;+ R;.,) 
2 
Equation (12) then takes the form 
6EUy;= 
sec® 
= = Xj) + Vain, +(1 +Aj) 


A; n=1 
+a similar expression with ¢’s instead of v’s 


The variable part of 6EUx; may be expressed in 
matrix notation as 


+a similar expression with @’s instead of w’s, 


and hence 

6E 1 || | 
1 2(1 J LYn, 


6. Bending Moment, Shear, and End Load in 
Frame j 


It is usual in investigations of this nature to take 
into account only the strain energy of bending of the 
frames in their own planes, since it is argued that. unless 
the depth of the frame cross section is appreciable com- 
pared with the radius (R,), the strain energy of end-load 
and shear is likely to be small compared with the bend- 
ing strain energy. 

It was originally decided to make this assumption 
with the refinement that allowance would be made for 
any offset of the neutral axis of the frame cross section 
from the skin line, since this may have an appreciable 
effect on the frame bending moment. Further investi- 
gation, however, showed that if the offset of the neutral 
axis was allowed for, the strain energy of end-load in 
the frame should be taken into account at the same time, 
for otherwise the computed stress distribution might be 
seriously in error. Finally, it was decided for the sake 


of completeness to include the shear strain energy of the 
frame. 

It is worth noting at this stage that, although frame 
bending moments will generally be less than would be 
obtained by the application of the “elementary” 
theory (i.e. assuming a sinusoidal distribution of skin 
shear), the skin shear—and therefore the end load and 
shear in a frame—will probably be considerably higher 
than the “elementary” theory would predict in the 
neighbourhood of a concentrated applied load. Thus, 
the bending strain energy will tend to be a lower pro- 
portion, and the end load and shear strain energy a 
higher proportion, of the total strain energy of the frame. 


6.1. THE LOADING APPLIED TO THE FRAME 
The loading (q,;) applied to frame j by the skin is 
given by 


(bay j) gR* (bay j 1) 
Substituting for the gR*’s from equation (8), and 
dividing throughout by R;, 
_R; (Rx 


sin né 
Vn, i+ 


Rj,, cos né 
On, i+1 


L; n 
=(say) sin n#+Q,,; cos né) . (8) 


In addition to the skin loading there may be external 
loads applied to the frame. These may be resolved into 
radial and tangential components, and moments in the 
plane of the frame, and it will be assumed that each of 
these types of loading may be expressed in the form of a 
trigonometrical series of the type in equation (18). As 
will be seen later, even point loads or moments may be 
dealt with in this manner. The term “loading” means 
the amount of load or moment applied per unit length 
of circumference measured at the skin line. 


Let S (Py sin nd + cos n#) (19) 
n=0 

n=0 

and = (21) 


The total tangential loading applied to the frame is 
(qr) + 


6.2. THE EQUATIONS OF EQUILIBRIUM OF THE FRAME 
The frame (Fig. 3) is assumed to have a finite depth 
of cross section which is fairly small compared 
with R). 
Consider an element of the frame between the radii © 
6 and 64-66. 
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TR 


VIEW LOOKING ALONG Oz (AWAY FROM 0) 


FiGureE 3. 


Resolving tangentially, and proceeding to the limit 
gives 


d 
Resolving radially gives 
dS 
-T=0 (23) 
and taking moments about C gives 
dT dM 
My, Rj + + Gri) Ri? + (Rj - hj) (24) 
From (22) and (24), eliminating dT / dé, 
dM — 
(R\-h)S= + (qr) + Gri) Rihj (25) 
and then, from (23) and (25) 
aM d 
(R,;-h) T= We + do Gri) +R; 
+ Rj(R;- Nr; (26) 


These last two equations enable § and T to be found 
when M is known. 
Further, from (24) and (26), eliminating 7, 


dM dM ad = 
at + Gri) Ri + Rjh; de (dei + Gri) + 
dNr; d?Mp; 
+ R;(R;-h,) + (M+ R;=0 (27) 


6.3. THE SOLUTION OF EQUATION (27). 
BENDING MOMENT 
Substituting from equations (18), (19), (20) and (21), 
the result may be put in the form 


+2 + Rj - n(Rj- h,) K,; 


1 
— (P,jh; + + sin 
+E (Qui + Qui + Hn) Ry +n (Rj hy) 


(Q,.;h; H,,,R))) cos né = 0. 
(28) 


THE FRAME 


By considering the equilibrium of the complete 
frame it can be shown that the n=0 and n=1 terms in 
this last equation are missing. 

Thus, resolving the load on the frame in the direc- 
tion XX (see Fig. 2), 


{(4r; +4ri) sin 6 Nr; cos 6} R,;d6=0, 


whence ;+P,,;- K,,;=0, while the coefficient of sin 
in equation (28) is (P,,;+P,,;- K,,;)(R;- h)). 
Similarly, resolving in the direction YY, 


® 


{(qri + Gri) cos 4 + Nr, sin 6} R,d6=0, 


whence 0, :+0,,;+,,;=0, 


while the coefficient of cos @ in equation (28) is 
(R; hi). 
Finally, taking moments about C (see Fig. 3) gives 


| + Gri) Rj +M;;} R,d6=0, 
+0,; +H, ;=0, 
while the constant term in equation (28) is 
The general solution of (28) is therefore 


M=-% {(Poj + + Ri - n(R;- h) 


n=e 


whence 


{Qui + Qaj + Hj) Ry+n(Ry+ hy) In — 
= sin né 
-n (Q,,;h; +Q,,;h; H,,;R;)} n(n? —1) 
+(m, +m, sin cos 4) . (29) 


The values of m,, m, and m, may be found by con- 
sidering the deflection of the frame. The equation for 
the radial deflection (w) is 


EI (ws 


j dw 
(R; hy 


and hence from equation (29), 


d’w (R; - 
—— +w= — 6 6 
+w EL, (m, +m, sin 6+ m, cos 6+ 


+ sinn#, cos né terms), 
which has the general solution 


(R; ( 1 
El, my COS 6 + 


+ m,6 sin 6+ sin cos né terms ) 


+ a complementary function of the form 


(w, sin cos 4) 
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which merely represents a “rigid-body ” displacement 
of the frame. 

Now w must be a periodic function of 6, of period 
2, which leads to the conclusion that both m, and m, 
must be zero. It is assumed here that the frame is 
initially unstrained, for otherwise a cut made at any 
cross section would give a discontinuity in w at that 
section, from the unstrained shape. Alternatively, w is 
defined as the radial deflection from the unloaded state, 
rather than from the unstrained state. The quantity M 
is then the bending moment due to the applied loads, 
and does not include any initial bending moment which 
may be “built in” to the frame during manufacture. 

Finally, the value of m, may be found by considera- 
tion of the circumferential strain in the frame, along its 
neutral axis. 

This is given by (w+dv/d6)/(R;—h;), while it is 
also equal to T/E divided by the cross-sectional area 
of the frame. The cross-sectional area is most con- 
veniently expressed here as the moment of inertia (/;) 
divided by the square of the radius of gyration (;) of the 
section, so that 

dv _(R;-h)) 
wt EL, 
Now « must be a periodic function of 6, of period 2= 
(assuming v to be the tangential deflection from the 
unloaded state) so dv /d6 will have no constant term. 

Thus, the constant term in w must be equal to the 

constant term in 
(Rj - p?T 
EI; 
whence 
= pT 
n,=constant term in (R;-h) 
and from equation (26) this 
= pF RiNr; 
=constant term in (R,-h)’ 
Hence, from equation (20), 


m,= - (30) 

The frame bending moment is then given by equation 
(29), with this value of m,, and with m,=m,=0. In 
passing, it may be noted that point loads applied to a 
frame will give rise to a bending moment of the form 
(m, sin 6+ m, cos 9), and the reason why these terms are 
absent is that the applied loads have been assumed to 
be continuously distributed round the frame. Point 
loads may be dealt with by the present method, but 
only as a limiting case of a distributed loading. 

A useful alternative expression for the frame bend- 
ing moment, presenting a much less cumbersome appear- 
ance than equation (29), may be obtained as follows. If 
P,;=Q,;=0 for n=2, 3, etc., the bending moment in the 
frame will evidently be that due to the given applied 
loads, balanced by skin loading of the form 


;+Q,,, cos 9), 


i.e. by the loading of elementary bending and torsion 


theory. Denoting this bending moment by M, equation 
(29) then takes the form 


— (R; n*hj)(P,,; cos Q,,; sin 
M=M n(n? - 1) 


Equation (31) is often in practice more suitable for 
calculating M than equation (29), since the series for 


(31) 


= 


(M - M) is usually very rapidly convergent, much more § 


so than the series of (29), while M is given by standard 
formulae for the usual types of loading, or may be 


derived by a simple strain energy calculation on a | 


frame. 

Equation (29), however, is more convenient for the 
calculation of the strain energy terms which are required 
in the present investigations. 


6.4. THE FRAME END LOAD 

The tension in the frame is given by equation (26). 
Substituting for M from (29), and for the other quantities 
from equations (18), (19), (20) and (21) the following 
equation is obtained after some simplification : 


R; 


Kuit h, 


((G,,;+K,, Cos 6 J,,;) sin 6} + 


6 


sin né 


The n=1 term is put in this form by using the 
equilibrium conditions 
j= 


The alternative form for the frame end load, cor- 
responding to equation (31), is immediately derivable 
from the latter with the aid of (26) and (18), and is 


T=T+ (Pui cos n9 Qn, sinné) (33) 
where T is the tension in the frame as given by the 
“elementary ” theory. 

This also follows from (32), since all the “ barred” 


terms make up T. 


6.5. THE SHEAR IN THE FRAME 

The shear in the frame may be similarly obtained, 
using either of equations (22) and (25). The n=1 term 
is again simplified by bringing in the conditions of 
overall equilibrium. 

Leaving out intermediate steps, the required result is 


| 


sin cos 


+ (Pay + Proj NK 

The “barred” terms make up the shear (S) which would 
be obtained from the “elementary” theory, whence is 


derived the alternative form (corresponding to equations 
(31) and (33) ). 
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is: n° = 1)? J; mel hy) (Q,; O,,;)+ 


The quantities P,,; and Q,; are defined by (18), and 
are linear functions of Un,j., and Of on 
nis On,j.15 Tespectively. In what follows only the ¥-terms 
(and their associated loading terms) will be considered 
in detail, since the ¢-terms may be written down imme- 
diately from these by comparison of equations (36) 
and (37). 

Using matrix notation 


R; R 
Ynint (38) 
(n=2, 3, etc.) 
The braces { } are used here to signify a column. 


War } 

L; 


ttt the a side of this equation by 
UW yp;/dv,, and using (36) and (38), 

7; 


L,/°’ IJ; 


ick 


n* (n? - 1)? 1; 
(Rj 1°h;) Pn; n(R; - Ky; (n? — 


Ri, 


6E- 


1) R,G,,} 


=, (#2, Ra} 
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7. The Strain Energy of Bending, End Load, _ 6% (Ri 
ie. 6E = 
and Shear of Frame j 1) 1; 
7,1. THE STRAIN ENERGY OF BENDING R...R, 
Let Uy»; be the strain energy of bending of frame j. on 
3 (R;- 3(R;- hj) | 2 Ri, Ris 1 Ries 
Then 6EU M* do, , ( Ry - ~ 
which using (29) * 
n,j+1 
1 (R; (Structural Terms) 
+R {(R; n?hj) (Proj + Prj) - n° - 1, ~ Pas 
hj) Ky; (n? - 1) + 
=P (n? - 1) {( n \(Q +Q )+ (Ris 
Therefore Rin 
n? (n? 1) I, {(R; n hj) +P,;) (39) 
K,,; -(n? - 1) R,G,,)} . (36) The asterisks in the square matrix indicate that the 
and matrix is symmetrical about its leading diagonal. 


The complete set of partial differential coefficients 
for frame j is given by equation (39), taking n suc- 
cessively equal to 2, 3,.... ete. 


7.2. THE STRAIN ENERGY OF END LOAD 
Let U;,; denote the strain energy of end load in 
frame j. 


2 
Then 6EU = | T*d6, 


which, using (32) 
p? + MP Ky)? + 2Qnj + Ini)? 


I; (n? = 
+n=0, | terms. 
Therefore 
6E 67 (R; hj) pi (n*Pp; + = nK,,;) 
: (40) 
and 
7; 
‘ (41) 


(for n> 1) 
As before, attention is confined to the v-terms, i.e. 
to equation (40), whence with the aid of (38) this result 
is derived : 


Terms) 
\. 
(42) 


(Loading Terms) 
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This equation may be written in the same way as 
equation (39). Clearly the only difference is in the values 
of the coefficients in front of the matrices, and to obtain 
the strain energy of bending and end load the corres- 
ponding coefficients of (39) and (42) are simply added 
together. 

In passing it should be noted that the “ simple beam 
theory ” expressions were used for the values of U,; and 
U;y;. These are good enough if the radius (Rj) is large 
compared with the depth of the cross section, and are 
probably acceptable for most of the frames likely to be 
encountered in practice. 


7.3. THE SHEAR STRAIN ENERGY 

Let Us,; denote the shear strain energy in frame j. 
This is the strain energy in the web of the frame. It is 
assumed here that the cross section of the frame consists 
essentially of two flanges which carry the end load and 
bending moment, connected by a web for taking the 
shear. The mean shear loading in the web will be equal 
to S, divided by the depth between the centroids of the 
two flanges. 

Then the element of strain energy in the web is 
approximately 


taking (R,;—h;) to be approximately the mean radius, 


(Rj hj) dé 
2Cty id? 


The quantity d)*ty;/d,; has the dimensions of an 
area, and may be regarded as the effective web cross- 
sectional area for the purpose of strain energy calcula- 
tion. For the idealised cross section consisting of two 
point flanges joined by a web, it is equal to the actual 
web cross section, since then d,,;=d). 

It is convenient to denote the ratio of the total flange 
area (=/;/p,") to the effective web area by the symbol »,. 


tw 


i.e. aU 


Therefore 


3Ep;?n, (R; - hj) 


6EU CI 


S°dé, 


which from (34) 


“Nj (R; - (P,,; a nK + On 
CI; ant - 1) 


+n=0, | terms. 


Therefore 
Wer; OT (Rj hy) (Pry + Pas — 
CI, (n? - 1 
and 
(R; (Q,,; + Qnj+nJnj) (44) 


CI, (n? 


—= 


Again confining attention to the v-terms, the follow. 
ing equation is obtained from equations (43) and (38) 


OU sr; 67 (Ri— hj) etc } 
CI, n? (n? - 19 
etc. | Vnijs Vn, ia} 
(Structural Terms) 


6E 


Cl n(n? - 
(Loading Terms) 


nK,,) Ri-, } 
L, 


A similar expression, for 6E s,;/99,, derives from 


(44). As with equation (42), equation (45), when written 


out more extensively, will take the same form as (39) | 


but with different coefficients in front of the matrices. 


7.4. THE TOTAL STRAIN ENERGY OF THE FRAME 


Let U,; denote the total strain energy of frame j, so | 


that 
Ur; = Up; Urr; 


Then the expression for 6E dU;;/0¥,, is obtained by 5 
the addition of equations (39), (42) and (45). writing the j 


last two of these in the same form as (39). It is con- 
venient for purposes of numerical computation to put 
the coefficients in front of the matrices in the form 
6=(R;-h))/[(n? 1)? multiplied by a polynomial in 
1/n, so that finally 


6E —— = 
dv, 

67 (Ri-hj) { _ 2 
(n? ne 2Rjh; C +h, + pj 
1 Rin Rj-, Ris q 
hak, 

any, -(F= 1 Rj, 
(Structural Terms) 
+ 

(n? nP,,; (or, i TP 

R; 2 ] Ri, 

_ 
L; | 
(Loading Terms) «ap 


3, ec.) 


(45) 
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Equation (46) is necessarily somewhat cumbersome 
in appearance, but it will be seen later that the com- 
putation of the various coefficients can be largely re- 
duced to a matter of routine by suitable tabulation. 

It will often happen in practice that the externally 
applied tangential loads are put into the frame along its 
neutral axis. When this is the case the coefficient for 
the loading terms simplifies to some extent, as follows. 

Referring to Fig. 3 and considering the piece of 
frame between the radii 6 and 6 + 66, 


h; + =0, 
the moment of the applied loads about the neutral axis 


being zero. 
Thus, + My; =9, 
or in the present case 

R; G,,; 0. 


Eliminating G,,, the loading coefficient then reduces to 


To conclude this section, the expression for 
6E0U may be immediately deduced from equation 
(46) by replacing the v’s in the latter by the correspond- 
ing and by replacing P,, by by and 
G,, by H,,,.. The alterations in the loading coefficient are 
obvious by inspection of equations (36), (37), (40), (41), 
(43) and (44). 


8. Loading Terms for Various Types of 
Applied Loading 


8.1. RADIAL LOAD + W aT 6=0 (Case A) 

A radial load at 6=0 may be regarded as the limit- 
ing case of a radial loading distributed over an arc, 
centred on 6=0, when the length of the arc tends to zero, 
keeping the total load finite. For this purpose, consider 
the “triangular” loading defined by the equation 


between =0 and (radians), being zero elsewhere. 


The total radial load is equal to W, while 


the requirement that dN;;/dé must be everywhere finite 
(since it enters into equation (27)) is satisfied except in 
the limit 8—>0. As a matter of interest: this require- 
ment appears to be an unnecessary restriction, since 


as far as Ny; is concerned equation (27) may be immedi- 
_ ately integrated with respect to 6. In fact, a concen- 


_ trated radial load could be treated as the limit of a load 


uniformly distributed over the are — 2 <@<8, as was 
done in Ref. 1. The final result must, of course, be the 
same. 


Returning to the present proposal, N;;R; now has to 


be expressed as a series of the form SK,,;cosn (see 
equation (20)). adi 


Thus 


0 


and K,j= 2 | N,;R; cos nd dé, (n>0) 
0 


B 
= B 


With a fixed value of 8, however small, n8—> 00 as 
n—> ©, and then K,,;—>0. In any case, however, the 
loading terms in equation (39) due to K,,; tend to zero as 
n—> oo, due to the factor (n? — 1)? in the denominator. 
Moreover, the structural terms for the frames (equation 
(39)) and the skin (equation (9) ) tend to zero as n—>©0, 
but, since the structural terms for the stringers (equation 
(16)) are independent of n, the coefficients of the v,,,’s 


remain finite. Hence v,;—>0 as n—> ©, and so, for a 
radial load W at 6=0, 


W W 
RADIAL LOADS +4W AT 6= +7 
This may be deduced from case A, since 


(48) 
8.2. 


scosno- yh. due 6=y, + 


} 
W 


W 
= — + — XScosnycos ng. 


Therefore 


= W W cos ny 


K,(n>0)= (49) 


- 


From (48) and (49) we immediately deduce case C. 


8.3. RADIAL LOADS +W aT 6=0 AND —W aT 6=7 
(Case C) 
This is case A -— case B with y=7. 
Hence 
. 
W 2W — 
K,,,;=0, K,,(n>0)=— (1 cos nz)= (50) 


Thus K,, ;=0 if nis even, and =2W/= if nis odd. Note 
that for this case both loads would be “downwards ” 
in Fig. 2. 


8.4. RADIAL LOADS +W aT 6=y AND - W aT 
(Case D) 
This loading may be expressed by a series of the type 
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Jui sin né (equation (20)), from which 


(51) 


as 2W sin ny 


sin n6dé = 


0 


This result may be obtained by a limiting process 
similar to that employed for case A, but is very simply, 
if somewhat uncritically, deduced from the consideration 
that, since all the load is concentrated at 6=y, the value 
of the integral is 


sin ny sin ny. 


0 


8.5. TANGENTIAL LOAD +W aT 6=0 (Case E) 
For this case take 
geiR; = cos no, 
n=0 


and by comparison with case A it is seen that 


(52) 


8.6. TANGENTIAL LOADS +4W aT 6=+y (Case F) 
By comparison with case B 


(53) 


WwW 
x’ 


With y=7/2 the loads are equivalent to a torque of 
amount WR, and 0,,; (n>0)=W cos (nz/2)/=, which is 
zero when n is odd and equal to (- 1)" * W/= when n 
is even. 


8.7. TANGENTIAL LOADS + W aT 6=y and —W aT 
6=—y (Case G) 
For this case take 
P.. nj sin no, 
n=1 
then from case D, 


(54) 


2Wsinny 
For y==/2 this becomes P,,,=2W sin (n=/2)/=, which is 
zero when n is even and equal to ( - 1)!" 2W/= when 
n is odd. 

Cases A to G cover the principal types of applied 
loading that are likely to be met with in practice. Cases 
of concentrated applied moments may be treated simi- 
larly, and it is not necessary to go into detail here. 


9. The Strain Energy Equations 


Expressions have now been obtained for the deriva- 
tives of the strain energy of each part of the structure in 
terms of the various unknowns, so that the formation 
of the strain energy equations can now be discussed. 
Let U; denote the strain energy of bay j (including 
frame j), and U the strain energy of the complete 
structure. 

Thus, 
U;=Us;+ Uy, +U and Uj. 


Assume that j takes the values 0, 1, 2,....J, J +1, 
the end frames being numbered 0 and J+1. The un. 
knowns are then the y,,;’s for j=1, 2,.. . . J; n=2, 3, 


etc. 

The remaining y’s (v,,;, ¥,,; for j=1 to J and y,,,, 
Un.s+, for n=0, 1, 2, 3, etc.) are known quantities, since 
the total end load and bending moment at every cross 
section, and the distribution of load on the end cross 
sections, are known. The case of a “rigidly built-in” 


end cross section will not be considered here since it is ' 


hardly likely to arise in practice. 

Similar remarks apply to the ¢-system of variables, 
and since in fact the whole of the following discussion 
will apply equally to the ¢- or Y-systems only the latter 
will be discussed as in the previous sections. The 
modifications to the loading terms in passing from one 
system of variables to the other are trivial, and are 
referred to at the end of Section 7. 

To continue 


_ Us 
6E W, (eqn. 9) + 
+6E 16)+6E (eqn. 46) . (55) 


Referring to the equations mentioned here, the most 
important point to note is that the strain energy deriva- 
tive with respect to any particular ¥ contains only v’s 
having the same value of n. This is a consequence of 
the representation of the stress distributions in the form 
of trigonometrical series, and of the uniformity of the 
structure in the circumferential direction. For example, 
the derivative 6E 9U;/0v,,;, which is obtained by taking 
the top row of (9), the top row of (16) and the middle 
row of (46), contains ¥,, j-,, Ynj, Un,j-, and terms depend- 
ing on the applied loads. 

The typical strain energy equation is 


oU 


6E =0, . . (56) 
and, in view of what has just been said, it is seen that 
6E 0U _ 6E0U;_, ,.6E0U; 6E9U;., 

Wri 


(in general, but some of these items will be missing at 
each end of the structure.) The expression 6E 0U;_, / 0¥,; 
is derived from equations (9), (16) and (46) by replacing 
j by j— 1, and taking the sum (bottom row of (9) + bot- 
tom row of (16) + bottom row of (46)), while the 
expression 6E0U;_, / 0€,; derives from equation (46) only, 
being given by the top row of the latter after replacing 
j by j+1. 

When all the various items are added together a re- 
lation is obtained between five consecutive v’s (i.e. five 
consecutive values of j) which may be conveniently re- 
presented in the form 


CE —— = Cn, Yn, + 
+ + +Cnj Vn, (57) 


The first five terms in this equation are the structural — 
terms, and the relationships between the coefficients of 
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and Yn, j+1, aNd Of Yq, and follow from the Similarly, taking j=2, equation (57) may be put 


symmetry of the matrix of the structural coefficients of 
the complete set of equations about its leading diagonal. 
This, in turn, is due to the circumstance that the strain 
energy (U) is a quadratic function of the variables. 


NOTE ON THE END CONDITIONS 
Equation (57) requires further consideration when 


j=l, 2, N-1 or N, since in these cases the loading 
applied to the end cross sections (j=0 and N +1) must 
be taken into account. The loading coefficient /,; con- 
tains only rib loading terms, being derived entirely from 
equation (46), applied to the three consecutive ribs j- 1, 
jand j+1. 

It is clearly sufficient to consider only the end j=0, 
and putting j=1 in equation (57) the two left-hand terms 
become Yn,-1 + Bn, Un, While for j=2 the left-hand 
term of (57) is Cn, 0 

The value of y,,, is known, but it is not clear what 
significance should be attached to y,,_,, since there is no 
“bay — 1.” wn»,-, May be regarded as a coefficient of the 
direct stress loading (Nj, ;-_,) at the remote end of a 
fictitious ““bay - 1” which is so long that it can have 
no constraining effect on the cross section at frame 0, 
ie. its length L_,= 90. Reference to equation (46) applied 
to rib 0, then shows that c,,_, =0, c,,; being the coefficient 
of vy,j.. In the top row of (46) (applied to rib j +1). 
Thus, for j=1 the two left-hand terms of equation (57) 
reduce to b, Un, 9 (note that b,,, will be affected by the 
assumption that L.,=°00), and the equation then 
becomes 


+ 1 + (ln, 1 + Dn, o Ye, 


the item b,,, Y,,. in effect constituting an addition to the 
loading term. 


in the form 


Dn, 1 Vn, Dn,» Yn,s +Cn, 4 +(ln,2 +Cn, 0 Vn, 


A particularly simple case arises when y,,,=0 
(n> 1), that is, when there is no direct stress loading on 
the end cross section, or when such loading is of the 
“simple beam-theory ” type. 


10. The Computation of the Strain Energy 
Coefficients a,; 5,; c,; and 
10.1. THE STRUCTURAL COEFFICIENTS @j, b,j AND Cn; 

The quantity a,,; is the coefficient of v,,;, in the equa- 
tion 6E dU /ov,,;=0, which will be called “ strain energy 
equation nj.” 

The left-hand side of this equation is made up 
of items from bays j- 1, j and j +1, each bay contribut- 
ing One or more terms which may be obtained from 
equations (9), (16) and (46). Thus, it is found that a,, 
contains skin and stringer terms from bays j- 1 and j, 
and terms from frames j—1, j and j+1. 

Similarly, 5,; contains skin and stringer terms from 
bay j, and terms from frames j and j+1, while c,,; is a 
single term from frame j+ 1. 

The algebraic expressions for a, and b,; are evidently 
very complicated, particularly when the strain energy 
of end load and shear in the frames is taken into account. 
It can be imagined that the numerical computation of 
these coefficients directly from such formulae could be 
a very troublesome business unless tackled systema- 
tically, bearing in mind that values are required over the 
range j=1 to J, and for n=2, 3, etc. (the highest value 
of n will depend on the type of loading and the degree of 
accuracy required). To this end, a method of tabulation 
has been devised which reduces the work largely to a 
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FIGURE 4. 


BASIC DATA: Skin thickness =0-048 in. throughout, E/C taken=2°5. No. of stringers =32; 
effective stringer cross section (including 1:5 in. of skin)=0-144 sq. in., giving A;=4-608 in.’. 
Tangential loads (at frames 2 and 8) applied on neutral axes of frames (see expression 47) 
¥,,,=B.M. about axis YY+2R,;=(, —9,402, —17,462, —16,977, — 12,732, —8,988, —5,659, 


—2,681, 0) for j=0 to 8. (% ;=%,;=9). 


For other basic data, see Table I. 
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matter of routine. The scheme is set out in Table I, the 
numerical quantities tabulated being for the example 
illustrated in Fig. 4, and it deals with the general case 
where the offset of frame neutral axis (A;), strain energy 
of frame end load, and shear strain energy of the frames, 
are all taken into account—in addition, of course, to the 
bending strain energy of the frames, the shear strain 
energy of the skin and the strain energy of the stringers. 
The scheme is (it is hoped) more or less fool-proof, and 
only a few words of explanation appear to be necessary. 
The chief point to note is that each coefficient (@,;, b,j, 
Cnj) is the sum of as many as five products, each of 
which may be resolved into two factors, one factor being 
dependent only on the structural data and the other 
being a function of n; one of these functions is, inci- 
dentally, a constant, which is given the value unity. 

The upper part of Table I is concerned with the 
evaluation of the structural factors, which are then tabu- 
lated in the lower part of the Table, to the left of the 
block of “multipliers.” The latter are the functions of 
n referred to above, their numerical values being given 
for n=3, 5, 7, 9 and 11. (The nature of the applied 
loading is such that ¥,;=0 for even values of n, but 
when the applied loading is unsymmetrically disposed 
with respect to the plane YOz the tabulation must in- 
clude both odd and even values of n.) 

The lower part of Table I is essentially a convenient 
device for the “mass-production” of the required co- 
efficients. The computation of the individual coefficients 
will be clear from the example 
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10° sin 

10°n sin 
( 

10° sin 


15°572 


5880:2 


a3) (n=3,j=1) 
(2:1043) (192-9) — (0:2919) (1736-1) + 
+ (0-0278) (15625) + (2673) (0-11111)+ 
+ (418-7) (1) = 1049-2. 


The use of a calculating machine with automatic 
multiplication and accumulation of binary products may 
be regarded as essential in all this work. With its aid, 
problems which would otherwise probably never be 
attempted may be solved in quite a reasonable time, 
particularly after practice has given the user confidence 
in the accuracy of his work. The calculation of the co- 
efficient a,, given above will be recognised by the reader 
familiar with matrix or vector theory as the formation 
of a scalar product. This is just the sort of work which 
the calculating machine revels in, and it is not being 
used to the best advantage if it is not given plenty of 
this to do. 


TABLE IL (Continued) 


70-131 | 


3 | 
— 2479-0 
109:64 
13-791 
3365-9 
— 131°493 


1,3 = (— 21486) (— 578:70) + (95-02) (— 5208-3) + (— 11-952) (— 46,875) + (760-7) (1736: 1) +(— 18-994) (15,625) = 2,333 x 10% 


76 
25-743 
61°362 


10.2. THE LOADING COEFFICIENTS 1,,; 


The quantity /,; is the term in strain energy equation 
nj due to the applied loading, and it derives from frames 
j-1, j and j+1, in the general case when all three 
frames are loaded. Thus, referring to equation (46), 
I,,=bottom loading term for rib j- 1+middle loading 

term for rib j + top loading term for rib j+ 1. 


In the present case, since the given tangential loads 
are applied on the neutral axes of the frames, the factor 
(47) is used in place of that given in (46). 

The calculation of the coefficients /,,; may, as in the 
case of the structural coefficients a,,;, b,; and c,;, be very 
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much simplified—one might almost say made possible 
—by appropriate tabulation of the various quantities 
involved, and the scheme shown in Table II has been 
designed with this end in view. The tabulation is based 
on the consideration that /,,; may be expressed as a scalar 
product of a set of quantities depending on the applied 
loads, with a set of “ multipliers,” the latter being func- 
tions of n depending partly on the type of frame loading. 
as evidenced by the presence of the factors sin n=/2 and 
sin? nz /2. 

The quantities depending on the applied loads are 
computed in the upper part of Table II. These are then 
transferred to the left-hand side of the lower part of the 
Table, and, with the block of “multipliers” on the 
right-hand side, the “ mass-production” of the loading 
coefficients may now be undertaken. 

It should be noted that the particular arrangement of 
Table II has been devised for the two types of applied 
loading (case G with y=7/2, and case C) to which 
the structure illustrated in Fig. 4 is subjected, but a 
similar scheme appropriate to any other type of loading 
should present no particular difficulty. An important 
feature of this table, which the given system of applied 
loads does not fully emphasise, is that it is particularly 
suitable when there are a larger number of loaded 
frames, as for example when three or more consecutive 
frames are subjected to tangential point-loads. The 
device of “diagonal summation” then comes into action, 
as shown by the arrows in rows 106, 107 and 108, the 
“diagonal sum” for columns j-1, j and j+1 being 
entered in column j. 


10.3. THE SOLUTION OF THE STRAIN ENERGY EQUATIONS 


As is shown in Section 9, a set of equations (j=1 to 
J) is obtained for each value of n, of the form given in 
(57). This makes it possible to solve for the unknowns 
v,») (n> 1) for each value of n separately, and is a very 
important feature of the method, since each set of 
equations is of a simple character readily admitting of 
speedy solution. 

The method adopted for the solution of the 
equations is due to Choleski, and the layout is 
exhibited in the left-hand part of Table III. An 
account of Choleski’s method is presented in Ref. 3, 
but for the sake of completeness a brief description of 
the method will be given here. The proof is most 
economically set out in matrix notation. 


Let the system of equations be written in the form 
Av+L=0, 


where A is the symmetrical matrix of structural co- 
efficients, and L the column of loading terms. 

Put A=TT’, where T is a lower triangular matrix, 
and T” is its transpose. As is seen below, the resolution 
of A into a product of this form can in fact be made, 
and the details of the process make it evident that the 
result is unique. 


Thus, TT’J+L=0, 
and writing (58) 
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Y is a column of subsidiary variables defined by (59), 
the original variables (¥) being derived therefrom by 
the solution of (58). Both (58) and (59) are triangular 
sets of equations which may be solved by direct substitu- 
tion, starting from the “tip” of the triangle in each | 
case, and the matrix of coefficients for one set is simply [ 
the transpose of the matrix for the other set. The | 
solution of the original equations is thus made to depend | 
on the resolution of the matrix A into the product TT’, 
i.e. on the determination of the matrix T. 


The coefficients of T and the terms in Y may be 
found by certain easily remembered rules whose deriva- ¢ 
tion is most simply illustrated by taking a moderately 
simple example. For this purpose, assume that there 
are only four unknowns, for which 


My, 
T= 
y, 
Ui 45 U3 45 u,,| 
Biss ae L, 
Also, let A= |, mai= 
As 
* * * Ay 


The asterisks in A signify a symmetrical matrix, and 
represent terms (the transpose of those above the lead- — 
ing diagonal) which in practice need not actually be 
written down. Then 
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Further, 

+ Ag, =Uy3 Uy 4 + Ugg + Ug, 
and Ly Y, Vo + Ug3 Yas 

so that - Us,” = 


33 
Finally, + +U,,? 
and L,=U,,Y, Yo + V3 + V4, 


A very economical lay-out is obtained if the co- 
efficients of T and Y are written below the leading 
diagonal of A, separating them from those of the latter 
by a thick zig-zag line, the loading column (L) being put 
at the right-hand side of A. This is shown in Fig. 5, 
which also includes a column e and a row s whose 
purpose will shortly be explained, and a column v for 
the final solution. 


From (60) we see that the column below A,, (not in- 
cluding s,) is the first row of the matrix [A, L], divided 
by (A,,)!. 

The column below A,, then follows from equation 
(61), the column below A,, from equation (62), and so 
on. The extension to any number of variables will be 
obvious. 


With the modern electrically-operated desk calcula- 


- tor the computations are, with practice, rapidly done, 
since all the intermediate operations for any particular 


coefficient are carried in the machine, and only the final 


_ answer need be written down. 


(61) 


The column e and the row s are used for checking 
purposes, the terms in s being the sums of the items in 
the corresponding columns (up to the zig-zag line). 


STRESSES IN A CIRCULAR FUSELAGE _ 545 


Thus, U,,S,=U,,(U,,+ ... 
from (60) 
or . (64) 


Similarly, from (60) and (61) 


Sy + Ug2 (Uy, + + +Y,)+ 


+ Ago +Ao,+Ao,+L., 
Or Aj, + Ags t+ Ly (Uy2 + S2)=0 
(65) 


Proceeding in this fashion, from (60), (61) and (62) 


+ + + Ya) +Uss + +¥a) 
= Aj5 + Aos.+ Ass 


or - 
— (U3 8, + Ss 5,)=0, . (66) 
and so on. 


The residual values of the left-hand sides of (64), 
(65) and (66), etc. are entered in column e as @,, @,, @;, 
etc., as each successive check is carried out. These 
entries are really an expression of satisfaction on 
the part of the computor at the completion of his labours 
on each column, since the check will not succeed if there 
are any errors in the column, and until such errors are 
eliminated the entry in column e cannot be made. 


Having completed the matrix T and the row Y’, the 
original variables v,, v., ¥, and ¥, now follow from 
equation (58). Thus, referring again to Fig. 5, 


Uy, Vv, =0, whence Was 


giving v,, and so on. 


An overall check on the ¥’s may be obtained by 
summing each row of 7 and forming the scalar product 
of these sums with the column of ’s. The result should 
be equal and opposite to X y, since 


Uy, + (Uy, + Uae) + (Uys + Ugg + + 
+ + + + Ug) t+ Yet Ys +y,)=0. 


Alternatively, a similar check may be made using the 
coefficients of the original equations: in this case it is 
verified that 


+(A,, + Age + Aas + 
+ + Ans + Ags + Asa) Ang t+ Aa at 
+(L,+L,+L,+L,)=0. 
The lay-out—corresponding to Fig. 5—for a typical 
set of strain energy equations (57) is shown in Fig. 6. 
This should be compared with the left-hand part of 


Table III, where five sets of equations have been solved 
according to this plan. 
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2 3 4 5 6 7 |Lx 10-4 10-9! nj 
~ 282-0 48-9 2333-2 | 
~282-7 62:0 4587-9 —12487 31 
1203-3. 493-6 151-4 5300-5] 29958 32 
~784:1 198-6 ~2021-5| —3°6457| 33 
2005°8 —1006-0 257-1 0:2352 34 
40-308 2385°8 —1281-9 0:7288 | 35 
5-239 22-489 43-042 2743-8 |-2980-4| 11144) 36 
6378  —26°450 44-761 1:5386 37 
72-032 —125-931 122472 -O-911 —14:129 —7:271 —68-868 | 
97-227 —101-078 146363 23245. 9-321  —24-107 | 
559°4 26:1 46 —80°7 | 
564-7 24:8 5-4 214-2} 51 
579-5 7-0 116 110-6) —0°3959 52 
~143 13:7 —2-9| 02064 53 
24-805 640°8 —27-0 15-9 00062 54 
0-282 -0°582 25-303 662°7 —40°5 0:0010| 55 
0-552. —1-054 25-715 680:°7| 140-0; —9°0127 56 
—1°549 26037 —0:2064 57 
—3-412 9-182 0-145 0-091 0-007 5-375 
21:538 34-183 19846 24630 24-968 24-173 31-412 
69:4 13 316 
68:2 16 ~58-7| —0-:0852 71 
61:8 35 11-0) 01382 72 
505°3 54:4 4-2 14-5} —0-0381) 73 
22-303 514-9 49-6 4-9 —0-0247 74 
0-159 2419 22-562 44°5 0-0028 75 
0-188 2-178 22-773 530°5| —41-0) —9:0067 76 
0-217 1-933 22-950 0:0778 77 
1-438 —2-897 0-907  —0-065 0-002 1-786 
26'634 22-119 25-843. 24892 24-708 21164 
84-8 0:5 
06 40-2) 0-0591 91 
80:8 14 -0-1165 
77-4 1:7 7:5, 0:0966 93 
21-267 471-5 75-0 2-0 —0:0334 94 
0-066 3-628 21-409 475°8 72:5 0-045 95 
0-080 3490-21532 479-6 | 19:3, 00057 
0-093 3-352 21-642 —0-0411 97 
0-801 2-064 —1910 0693 —O112 0-016 0-890 | 
24°537 27095 23-077 24880 24-900 22:532 
918 0-2 78 
21-036 443-9 91-1 0-3 —0-:0265 111 
4-364 20°612 89-5 06 2-8] 112 
0-010 4-418 449-2 87-7 0:8 3-8) —0-0138 | 113 
0-015 20:747 451°8 86:3 0-9 | —0-0062 114 
0-029 4-221 20:832 454-4 84-9 0-0021 | 115 
0-039 4-135 20-912 4566 | —10:8) —09:0050 116 
0-043 4-051 20981 0:0246 117 
0-371 —0°821 0-311 0-119 —0-025 0-005 —0°516 
25-781 24-224 25-322, 25-126 24-985 24-968 20:465 
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TABLE Ill (Continued) 


0 0 en) 0 0 0 
1 —9402|—1249] 161} —85| 59] —27 
(Ib.) 2 |—17462| 2996}—396| 138/—117| 43 
(See 3 |—16977\-3646] 206; —38| 97| —14 
Fig.4 | 4 |—12732| 235 6| -25| 
for 5 | —g9gs} 729) 1 3 5 2 
6 | —5659} 11141 —13} —7 
7 | —2681| 1539}-206| —41| 25 
8 0 o| =O 0 0 0 
0 9402} —32} —7 2 
1 8060|—1415} —32| 20] 
—(Ynj- | 2 —485| 2214/-120| 25| —24 5 
_y | 4245-1294) 40} -2) 14) -1 
: 5 | —3329] —128 3 1 0 1 
(Ib) | 6 | —2078) -142} 30) —3 
7 | —2681| 513) —41 2 
nP,; | 1 —969) 42]; —-9 4; 
| 2 2016] 12) +7 2 
(Ib.) | 3 | —2019} 50} 
Se | 4 | 649} —13} 0 
equa- | 5 | 13 1 1 1 0 
(38) 7 | 500) —33 -2 1 
8 ~3 1 0 
(q R*);L; 
| N, cos (Ib.) (Ib.) 
e\G| 2 3 7 i 2 3 
0 | —14798| —20372\— 1286 0 0 0 
5 | —14805| —20218|—1322| 374 482| —680 
10 | — 14824] — 19752-1422} 761 935) — 1342 
15 —14846| — 18969|—1570| 1174) 1333) — 1967 
20 | —14864| — 17869] -—1752| 1624| 1646! —2534 
25 | ~14874| —16480|—1960/ 2120 | 1849, —3020 
30 | — 14862) — 14860/—2189| 2670} 1921) —3409 
35 | — 14787) — 13096] —2429| 3278 | 1846| —3690 
40 | —14588 3943| 1627) —3862 
45 | —14202) —9521/—2830| 4653 | 1276| —3934 
50 —13590) —7849|—2933) 5389 —3922 
55 | —12755| —6295|—2953| 6132 284] —3842 
60 | —11718| —4863|—2890| 6862 | —299| —3712 
65 | —10485| —3564/—2741| 7562 | -—902| —3547 
70 —9018) —2431]/—2487) 8213 | — 1495] —3365 
715 | ~7248 —1508]~2092/ 8786 | —2037/ —3186 
80 | —5121) —821|/—1531| 9243 | —2479| —3034 
85 | —2662| -—347| —813| 9540| —2771| —2931 
90 0 0 0| 9644 | —2873| —2894 


11. Numerical Results for the Example of 


Figure 4 


The derivation of the unknowns (v,,;) in the strain 
energy equations is an essential step towards the 
complete solution of the problem, but the shear and end 
load distributions in the skin and stringers, and the 
frame bending moments, end loads and shears have still 
to be calculated. 

In view of the large amount of information which 
can now be obtained, a systematic tabulation of the 
required quantities is almost a necessity on the grounds 
of speed and accuracy. A few notes on this subject, 
illustrated by the example of Fig. 4, will therefore not 
be out of place here. 

First note that, in Table III, the loading column L 
has been divided by 1,000. This is done to reduce the 
loading terms to the same order of magnitude as the 
structural coefficients, the purpose being to keep the 
number of significant figures needed for the checking 
process down to a minimum. The column on the right 
of L x then gives v,; x 107°. 

The full values of v,; (rounded off to the nearest 
whole number) are then tabulated in the upper block on 
the right-hand side of Table III, and below these are 
entered the quantities (v,;—,,;.,)/n. These latter are 
the coefficients in the sine series for gR.L;/(R; Rj~,) (see 
equation (8)), while in the present example ¢,;=7;=0. 

The third block on the right-hand side of Table III 
gives nP,,;/(n* — 1), P,,; being derived from equation (38), 
and is required for the calculation of frame bending 
moments (31), end loads (33) and shears (35). 

The shear and end load distributions in the skin and 
stringers may now be computed from equations (7) and 
(8) respectively. For this purpose tables are needed of 
sines and cosines of multiples of 6 (6, 26, etc., up to 
116 or 126, say), over the range 6=0 to 180°, an interval 
of 5° or even 10°, being quite good enough for prac- 
tical use. Such tables do not seem to be readily avail- 
able in print, although they must surely be used 
extensively in harmonic analysis, but they can very 
quickly be made up if, as suggested, a fairly large 
interval of tabulation is acceptable. 

The calculation of any particular function, for a 
range of values of 4, is facilitated by writing the co- 
efficients of the series on a strip of paper so that, when 
placed upon these home-made trigonometrical tables, 
the proper factor will appear against each cos n# or 
sin n@, as the case may be. The sum of products is then 
obtained on the calculating machine and the answer 
tabulated or plotted as desired. 

In the present example the most interesting distri- 
butions are N;R; cos z at frames 2, 3 and 7, and 
(qR*); L;/(R; Rj+,) in bays 1, 2 and 3, and these are 
therefore put in Table III, below nP,;/(n?—1). The 
complete set of distributions was afterwards computed, 
and plotted as Figs. 7 and 8. The first of these diagrams 
shows N;R; cos z divided by v,,;. while the second dia- 
gram shows (qR°); L;/(R; R;.,) divided by j+,). 
The “elementary theory” values of these functions are 
respectively cos@ and sin#, which are accordingly in- 
cluded for comparison. 
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FRAME BENDING MOMENTS, END LOADS, AND SHEARS 

The most convenient formulae for numerical work 
are (31), (33) and (35), provided simple expressions can 
be obtained for M, T and S. For radial or tangential 
point loads these are more or less standard results, but 
it is interesting to derive them from equations (29), (32) 
and (34), by summation of the trigonometrical series. 
Only the cases required for the analysis of the frames 
in the numerical example will be considered. 


11.1. RADIAL LOADS SYMMETRICAL ABOUT THE PLANE 
XOz 


The general expressions for M, T and S for loading 
symmetrical about the plane XOz (the “ vertical” plane 
in Fig. 2) are obtained from (29), (32) and (34), respec- 


tively, by retaining only the coefficients K,; in these 
equations. Thus 


= 

a cos 4 _y K,,, cos né 67 
P=Ko,i+ (Rj-h,) azz (n?-1) 


The cases of interest here are A and C; for case 
A, K, ;=W/2= and K,,(n>0)=W/>, and hence the 
sum & {cos n#/(n*—1)} is required. Differentiation with 


n=2 


respect to 4 will then give the sum ~ { —nsinn@/(n*— 1)}. 


(n’—1). 
Now X {cos né/(n°?—1)} may be shown to be the 


n=2 


Fourier series for 
{4+14 cos @-4(=- 4)sin 6} 
in the range 0627, 


whence 
n=2 n- 4 2 

The term —p;°K, ;/(R;—h,) will be neglected, as 
being insignificant in practice. 

It is then found that, for case A, 


(68) 


_ taking 4 in the range from 0 to 2. 


Note the discontinuity in § at the two ends of the 
range, due to the applied load. 

For case C, K, ;=0 and K,;—2W/= (n odd). Thus 
the sum {cos né/(n°—1)} is required (odd terms 


only). 


Now this latter series may be expressed in the form 


ly 
n*- 


and it is hence the Fourier series for the function 


4{4+3 cos 6-4(x-4)sin6} — 
-4{4+4 cos (x - 4) 46 sin (= 6)}, 


1 si } 


the range now being O<6<7. 
Differentiation with respect to @ then gives 


i.e. for 


0<6<r. 
Thus for case C, 


Ww ( 2h; 
and R,- 


=| 


cos 6+ (6- sino} 
)cos #- (6- =) sine } 
+ 


taking 6 in the range from 0 ‘io t 


(69) 


11.2. TANGENTIAL LOADS SYMMETRICAL ABOUT THE 
PLANE XOz 
We are here interested in case G, with y==/2, 
giving P,;=2W sin(nz/2)/z. The applied loads are on 
the neutral axis of the frame, whence G,;= - P,,;h,/R; 
(see derivation of (47)). 
Equations (29), (32) and (34) then reduce to 


sin — cos né 
_ 2W(R;-h) 


= as 
nsin cos né 
Ly hy cos 4 2 
R;- hj n=3 2 
-1 
sin sin no 
and §= 2” 2W j - 
R;- h, n=3 


n*-1 


Now > (odd terms only) = 6+ 
+ 5) siné, 


Integrate both sides of this equation with respect to 
4. between the limits 0 and 4, and then replace 6 by 
=/2—6. This leads to the result 
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Differentiation with respect to 6 now gives 


A further differentiation will give the series required 
for T; this is confirmed by equation (23), which is now 
T =dS/d6. 


_ WAR; 


-h) (= ) 
Hence M (5 cose 


T= { + a} 


and S= , 


(70) 
in the range —-7/2<6<7/2. At 6=+72/2, T is dis- 
continuous, being equal to +W/2 and -W/2 on the 
upper and lower cross sections respectively. 

In the range =/2<(6<37/2, M, T and S are obtain- 
able from equations (70) by reflection about the 
horizontal diameter (6==/2), with a reversal of sign in 
the case of the first two of these quantities. 

Returning to the numerical results for the example 
of Fig. 4, formule can now be obtained for the 
bending moment, end loads, and shear in any of the 
frames. The unloaded frames—i.e. those not subjected 
to any externally applied loads—are only lightly stressed, 
as can be seen by reference to the values of nP,,;/(n? — 1) 
in Table III. The simplest criterion to take here is the 
sum of these coefficients, which gives the maximum 
tension in the frame (at 46=0). 

Of the loaded frames, the ones which show the 
greatest departure from elementary theory are frames 2 
and 3, and this section is concluded with a brief dis- 
cussion of the results for frame 3. 


For j=3, 


R,=37°5 in., 


kj=2°9 in., = 3183 lb. 
Equations (69), (31), (33) and (35) then give 


M=( 3183) (34-6) { + (o- >) sino} 


— (1:267) ( 2019) cos 36 + (1-4) (50) cos 


T =(- 3183) { 1-6676 cos (o- >) sin } 


— 2019 cos 36 + 50 cos 54 
and 


S=(- 3183) { 0:6676 mee 5) cos 6 } 


~ (2019) cos 36 + (50) cos 


Terms beyond n=5 have been neglected. 
The maximum (numerical) value of M occurs at 


#=0, and is - 52,440 Ib. in. This is 95 per cent. of the} 


corresponding value of M (- 55,070 Ib. in.) 
The value of T at 6=0 is 


may be compared with the maximum (numerical) value 


of T, which occurs at about 6=24° and is approximately | 


— 6,340 Ib. 
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Some Problems Associated with 


Stress Concentration 


H. L. COX, 


M.A., F.R.Ae.S. 


(Senior Principal Scientific Officer, National Physical Laboratory) 


HE theory of stress concentration is such a 

fascinating study that the whole of this lecture 
could be most happily devoted to discussion of the stress 
distribution round all sorts of special boundaries. At 
the same time the most detailed knowledge of elastic 
stress distribution is of little value unless we can be sure 
how the stresses computed will be reflected in the actual 
strength of the part. For that reason this paper is 
divided into three sections: the first deals with certain 
moderately simple cases of stress and load concentration 
and shows how the stress analysis is confirmed by static 
and fatigue tests: the second describes some investi- 
gations of special two-dimensional boundaries and 
indicates some fairly general theoretical conclusions 
which can be drawn: in the third I shall try to review 
the possible reasons why, in practice, the best found 
conclusions are not always borne out. 


Stress Concentrations 


Figure | shows a straightforward stress concen- 
tration, a hole in a strip. Photoelasticity shows where 
the maximum stress is developed and its magnitude: 
moreover from the pattern of fringes the stress distribu- 
tion can be plotted over the entire field. The concen- 
tration of stress arises simply because the stress lines 
have to flow round the obstacle, the hole; the next case 
of a rivet transmitting a load to a plate (Fig. 2) is 
different because the stress lines have in effect to flow 
into the hole. I usually try to distinguish the two cases, 
because concentration of load is really a separable 
effect. For instance, if we cut away a corner of the 
plate, leaving the rivet supported on one side only, the 
least possible multiplication of stress at the side of the 
hole is by about 10, when the depth of section is about 
twice the rivet diameter; but the stress in a plain beam 
of the same section under the same loading would be 
5:5, so that the real stress concentration due to curva- 
ture of the beam is only 10/5-5 or about 1:8. Fig. 3 
shows both stress and load concentration together: the 
point again is, of course, that a good part of the high 
concentration near the point of application of the load 
is inevitable, whereas the concentration in the fillet is 
evitable or nearly so. 

Whether the concentration is avoidable or not, how 
does it affect the strength? Fig. 4 shows a standard 
B.S.I. sling hook, with the maximum stress a little 
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above the horizontal section; Fig. 5 shows this and 
other hooks after static impact and fatigue tests. Under 
static and impact tests the hooks just opened out and 
the stress concentration is evidenced only by the 
marked extension in these regions; under repeated load- 
ing, however, cracks developed, and in the sling hook 
the crack lay precisely at the point of highest stress. 
(I may add in parenthesis that we did not expect this 
result when we put the test in hand; we designed the 
hook with great care all the way from the bed to the 
horizontal section, but we designed the shank, between 
the horizontal section and the eye, by eye. As you see 


From Frocht's ‘ Photoelasticity '’ published by John Wiley 


Stress pattern of a bar with a central circular hole 
subjected to an axial tension. 


FiGure 1. 
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From Frocht's Photoelasticity,’’ published by John Wiley 
FiGuRE 2. Photoelastic stress pattern of the head of a bar 
subjected to a tension by means of a single central rivet with 
a snug fit. 
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FIGURE 3 
(above). 
Stress 
distribution 
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piece end 
fitting. 
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FIGURE 4 
(left). 
Stress 
distribution 
in B.S.I. 
sling hook. 
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Courtesy of the Institution of Mechanical Engineer; 
Figure 5. Hooks after static, impact and fatigue tests. 
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our eye was not entirely “in” and the stress at the 
point of failure was in fact about 5 per cent. greater 
than the stress at the horizontal section. In the Liver- 
pool hook the failure was in the right place and in both 
hooks it occurred after endurances of the order to be 
expected in relation to the computed stress ranges.) 


Another case from the same field is shown in Fig. 6 
which shows the stresses developed in eyebolts under 
direct lift and oblique loading. Fig. 7 shows one eye- 
bolt broken under repeated oblique loading and another 
broken under repeated direct lift. The computed 
stresses at the crown at the extrados and at the sides at 
the intrados in the latter case were about equal and, as 
can be seen, the eyebolt cracked in all three places. 
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FIGURE 6. Stresses in eyebolts. 
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Courtesy of the Institution of Mechanical Engineers 
FiGure 7. Fatigue cracks in eyebolts. 


Moreover the endurance was again just about right for 
the computed stress range. 

Figure 8 shows the stress distribution in ordinary 
and studded chain links and the discontinuity in 
stress at the end of the parallel portion again illus- 
trates stress concentration superposed on load con- 
centration; the latter is of course due to the bending 
moment associated with the offset load, while the 
former—the sudden jump in stress factor—arises from 
the curvature of the end of the link. We have done 


\o-E > 
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T if | H 
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Figure 8. Variation of stress round open and studded chain 
links of same dimensions. 
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Distribution of shear stress round a rectangular 
section under torsion. 


(a) Straight bar. 
(b) Bar bent into helical spring. 


FIGURE 9, 


quite a number of fatigue tests on chains, and failure 
always occurs at about the expected endurance for the 
computed stress range, either at the crown at the extra- 
dos or at the intrados at the end of the straight side. 
It is to be noticed that we never get fatigue cracks 
at the intrados at the crown, although the stress range 
there is far higher than in the other two places; pre- 
sumably the compressive mean stress is the reason. 

Another example of concentration of stress due to 
curvature is shown in Fig. 9, which represents the dis- 
tribution of shear stress over the section of a rectangular 
section spring. Fig. 10 shows the cracks developed 
under repeated loading; again the cracks occur in all 
three regions of highest stress, but it is to be noticed 
that these cracks propagate on planes of maximum 
shear stress (whereas those in the lifting gear propagate 
on planes of maximum tensile stress) and also that 
they advance only slowly into regions originally 
subjected to lower stresses. 

At this point it may be remarked that much evidence 
has recently been compiled that cracks often start in 
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FicurE 10, Positions of cracks in rectangular section springs. 
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Courtesy of the Institution of Mechanical Engineers 
Ficure 11. Fatigue crack at fillet in a wide thin plate. 


Courtesy of the Institution of Mechanical Engineers 
FiGcure 13. Fatigue crack starting at fillet instead of 
central hole. 
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Fictre 15. Concentration of direct stress due to 
polygonal holes. 


FIGURE 12. 


Ficure 14. Fatigue cracks spreading from small 
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Detail of a fatigue crack near the fillet in a wide 
thin plate. 
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FiGureE 16. Concentration of shear stress due to 


polygonal holes. 
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FiGuRE 17. Forms of deep, narrow groove. 


regions of high stress after endurances which are 
_ reasonable in relation to the computed range of stress; 
, but that the rate of propagation of such cracks depends 
very much upon the stress gradient, that is, the rate at 
which the stress falls off from the point of maximum 
value; and most remarkably that the cracks often seem 
blissfully unaware that they themselves ought to be 
having a tremendous influence on the stress distribution 
in their neighbourhood. 

The unerring way in which cracks find regions of 
high stress is illustrated by Fig. 11 which shows a 
fatigue crack at a fairly generous fillet in a wide thin 
plate. Fig. 12 shows the detail of a similar crack: the 
stress concentration was only about 1:20 but the 
gradient was not very steep, so that this particular crack 
formed away from the edge. Fig. 13 illustrates how this 
gradient effect can favour cracking at a fillet with a 
stress factor of 1:2, instead of at a hole where the stress 
factor was 3 but the gradient was very much steeper. 
The stress gradient effect is illustrated again in Fig. 14 
which shows cracks which have spread only rather 
slowly from small diametral holes. 

From this brief survey of the effect of stress con- 
centration in causing fatigue fracture it will appear that 
we may need to know not only the maximum stress at 
the discontinuity of section, but also something about 
the rate of change of the stress in this neighbourhood. 
That question is bound up with the general problem 
of stress analysis for discontinuities more complex than 
the simple holes and fillets which I have dealt with so 
far. 


Special Two-Dimensional Boundaries 


_ General analysis in two dimensions is not particu- 
larly difficult and some results are shown of an 
_ investigation relating to holes of quite a variety of 
forms in an infinite plate. For holes in plates of finite 
width the photoelastic method would have to be used, 
but the general theory suffices to indicate the effect of 
form of hole. Fig. 15, for instance, shows the concen- 
tration factors under uniaxial tension of a wide range 
of polygonal holes, of which one is illustrated. On this 
figure the line 1+2/(a/p) is marked; this represents 
the stress factor for an elliptical hole, and it will be 


seen that it affords a surprisingly good approximation 
to the true factor, even for holes of quite different form. 
Fig. 16 shows the factors for the same holes under shear 
parallel to the axis of the hole; here it will be seen that 
the formula »,=1+ /(a/p) for an elliptical hole is far 
wider of the mark. There is one special class of poly- 
gonal holes which have p=4a alternating with p=00, 
such as a square with flat sides but rounded corners. If 
the number of sides of such a hole is increased without 
limit, the hole becomes indistinguishable from a true 
circle, yet the formula 1+2,/(a/p) would indicate 
a stress factor 1+2+/2, whereas the formula on which 
the figure is based indicates the obviously correct 
value 3. The important conclusion here is that abrupt 
changes of curvature in themselves do not affect stress 
concentration. 

Figure 17 shows a range of rather deep grooves 
having various curvatures at their ends, and Fig. 18 
shows that again 1+2¥/(a/p) affords a good estimate 
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Figure 18. Stress concentrations due to deep, narrow grooves. 
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Figure 19, Forms of 
square holes. 
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FIGURE 20. Stress concentrations due to square holes and 


90° Vee grooves. 
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of the stress concentration under direct stress, whereas 
1+ V(a/p) for the shear cases is far less accurate. The 
same point is further illustrated by Fig. 19 which shows 
a range of square holes and Fig. 20 which shows the 
corresponding stress factors. The curve of «, for 90° 
Vee grooves shown on this figure was computed by a 
different method; that this curve should lie appreciably 
below that for square holes is surprising, because the 
form of the groove differs very little from that of the 
square hole. This difference is difficult to understand, 
but the V-groove curve agrees very well with values 
obtained by Griffith using the soap film method 
(Fig. 21). 

The formal solution for these special boundaries is 
easy enough but the detailed working is tricky. On the 
other hand the solution of the torsion problem for any 
boundary is entirely straightforward. For instance, 
Fig. 22 shows the stress distribution in a hollow lobed 
shaft, and Fig. 23 the stress in a splined shaft. Such 
examples could be multiplied without limit, but Fig. 24 
shows an application which is at once far less elegant 
and far more useful. This was an early check on the 
general method by application to a square shaft. Using 
a rather simple means of representation I finished up 
with a shaft with wobbly sides and well rounded 
corners. However, I was a bit disappointed when I 
found the stress factor varying from about 1:1 to 1-6 
along the middle of the side, instead of keeping close 
to the value 1-35 which Timoshenko assured me was 
the right answer. Then I had the happy idea to 
“correct” the computed values at each trough by 
dividing by 1+ /(h/r), where h is the depth of the 
local trough and r is its radius of curvature, and the still 
happier idea to divide the value at each crest by 
1- V(h/r), where h is the height of the crest and r its 
radius of (convex) curvature. By this means the factor 


Figure 23. Stress distribution 
round a sharp corner in a six- 
splined shaft. 
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FiGurE 22. Distribution of stress in a hollow lobed shaft. 


was brought into line with Timoshenko. Fig. 25 shows 
the same procedure applied to a round shaft with a 
single flat. In this case there is no exact solution for 
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Ficure 24. Distribution of 
shear stress along the side of a 

-6h square shaft under torsion. 


x 


Distance From centre of shaft 


06 os 10 


Crown Copyright 


55 
Uncorrected stress factors 


Stress factors corrected for local 
Stress concentration 


1824 


Figure 25.  Distribu- 160° 
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a round bar with a flat. 
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in tension. (Figures against points are values of 2a/d.) Figure 30. Fatigue failure of crankshaft starting at oil hole. 
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FiGure 31. Crack propagating from a small central hole in a 
wide, thin plate. 


comparison but I think that it is clear that the corrected 
values must be pretty near the mark. 


Practical Aspects 


One of the main reasons for this general investiga- 
tion of form of hole was to determine the stress con- 
centration due to a small crack or flaw at the root 
of a notch. If the two separate factors for the flaw and 
for the notch were always to be multiplied together the 
presence of flaws would not affect notch sensitivity; but 
if the two were not always multiplicative, the effect 
would be to reduce notch sensitivity because notch and 
flaw together might cause only slightly higher concen- 
tration than the flaw alone. This application of the 
general analysis proved rather difficult but the results 
shown on Fig. 26 for the factor in shear are probably 
not far from the truth. I have plotted the apparent 
stress factor, that due to notch and flaw together divided 
by that due to flaw alone. To explain limitation of 
notch sensitivity in this way it is necessary to postulate 
depths of flaw d appreciable in relation to notch 
depth a. On the other hand, the effect is to render the 
stress factor very insensitive to notch root radius, which 
was the experimental finding I was trying to explain. 
The main objection is of course that such long cracks 
ought to be discernible by some means or other, but 
there are other serious objections to the flaw hypothesis 
which it will not be easy to surmount. However it 
seems quite possible, in relation to recent observations 
on cracks which do not propagate, that the old Griffith 
idea may shortly gain a new lease of life in a modified 
form. 

The theoretical results referred to all relate to holes 
in infinite plates or to cases of shear or torsion of long 
bars of uniform section. For other cases photoelasticity 
has to be used and for truly three-dimensional problems 
the frozen stress technique is needed. My colleague 
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Ficure 32. Fatigue crack changing from transverse plane to 
two inclined planes. 
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FiGure 33. Fatigue crack changing from two inclined planes 
to plane of maximum shear. 


Brown has recently applied this method to two cases 
of round bars with circumferential V notches under 
tension across the notch. The results are shown in 
Fig. 27. It will be seen that the sharp, shaliow notch 
causes the higher stress concentration but that the stress 
falls off far more rapidly, so that the stress distribution 
in the core of the bar under the notch is not very dif- 
ferent in the two cases. The figure shows also the 
radial and hoop stresses under the notch, and the shear 
stress in the core is represented by the difference 
between the p and r curves. It will be seen that the 
ratio of p to the difference (p-—r) is greater for the 
deeper notch. This ratio—if we can judge its average 
value—is important because if p reaches the limit- 
ing value for cleavage before (p—r) reaches the limiting 
value for plastic deformation, brittle fracture will result. 

Since photoelastic investigation by the frozen stress 
method is a long laborious task, Brown took a short cut 
to the desired answer by testing a number of notched 
pieces of mild steel. In Fig. 28 he has defined the 
triaxiality ratio average p/average (p—r) as the ratio of 
the one per cent. proof stress for the notched bar to 
that for the plain bar, which is as good a definition of 
triaxiality ratio as any other. The effect in these tests 
of the root radius of the notch was quite unimportant; 
it was the depth of the notch that mattered. It will be 
seen that a notch leaving a core only 20 per cent. of 
the diameter of the bar gave a triaxiality ratio of 2°5. 
This final composite diagram (Fig. 29) shows the con- 
clusions drawn both from the tests on mild steel and 
from the photoelastic analyses with a_ satisfactory 
measure of agreement between the two approaches. 
Neuber’s results appear rather different, but comparison 
between the infinitely deep notch and notches of finite 
depth is very difficult. 

In conclusion I want to return once more to the 
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FicurRE 34. Similar change of mode of crack propagation at 
opposite side of test piece (see Fig. 31). 
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most interesting problem of crack propagation. Fig. 30 
shows a fatigue crack developed from an oil hole in a 
crankshaft. The crack propagated along the plane of 
maximum tension as a result of alternating torsional 
stresses arising probably from torsional vibration. This 
might be described as a classic example, in the sense 
that everyone who feels that his upbringing ought to 
have been classical, will pretend that he completely 
understands it. In contrast, Fig. 31 is a crack develop- 
ing from a small hole in a wide thin plate of aluminium 
alloy under repeated tension. The cracks propagated 
almost equally on the two sides—a most encouraging 
advertisement of the testing technique—suddenly both 
advancing cracks changed their mode of propagation, 
speeded up their rate of progress and cut the picce in 
two. This final stage could be described as failure under 
the static mean load, but still I believe that even in this 
stage the range of load was playing a part—a compro- 
mise description might be “static failure under the 
maximum tensile load in the cycle.” 

However, the change of mode of propagation 
apparent in this figure is not the first. There was 
another change earlier which appears in Fig. 32. The 
crack from the side of the hole starts to propagate on 
the plane of maximum tension, procceds only a few 
millimetres in that manner and then transfers rather 
abruptly to a mode, which pays some, but not entire, 
regard to shear. The cattle shed roof is quite character- 
istic and the angle of pitch is about +30° from the 
horizontal or transverse plane; precise measurement of 


the angles is difficult but they are certainly far removed 
from 45°. The second change of mode of propagation 
apparent on the previous figure is shown in Fig. 33 and 
again in Fig. 34. The change here is very abrupt indeed 
from cattle shed roof to a single plane of shear at 
ebout 45°. 

There is certainly no need at the present time to 
urge that more attention be paid to crack propagation— 
we are all watching cracks extend or not extend for all 
we are worth; but if I may draw my own moral from 
my own talk I should urge that more attention should 
be paid to stress distribution round the advancing crack. 
At present the common practice is just to remark that 
the stress at the tip is very high or infinite and to leave 
it at that: we should do better. 


ACKNOWLEDGMENTS 


Acknowledgment is made by the author to the 
Controller, Her Majesty’s Stationery Office for per- 
mission to reproduce Figs. 15-26, which are Crown 
Copyright and are reproduced from R. & M. 2704; and 
to the Institution of Mechanical Engineers for Figs. 5-8, 
11-14; the Institute of Physics for Figs. 27-29, which are 
from a paper by A. F. C. Brown published in the 
“British Journal of Applied Physics” 1954; and the 
Institution of Mechanical Engineers for Figs. 31-34. 
Figs. 1 and 2 are reproduced from “ Photoelasticity ” by 
Max Frocht, published by John Wiley & Sons Inc., to 
whom acknowledgment is also made. 


ers 
ines 4 
ses 
der 

in 
tch 
ion 
dif- 
the 
ear 
nce 
the | 

the 
age 
nit- 
Ling 
ult. 
ress 

cut 
hed 

the 
of 
r to 
1 of 
ests 
ant; } 
| be 
. of 
3 
and 
tory 
hes. 
ison 
nite 

the 


562 VOL.. “59 


JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY 


AUGUST 1955 


TECHNICAL 


NOTES 


Contributions to this Section of the JouRNAL will be eligible for Journal Premium Awards 
and will normally be published within two months of being received. 


The Effect of Curvature on the Centre of Shear 


W. 


JOHNSON 


(Department of Mechanical Engineering, University of Sheffield) 


CONSIDERABLE AMOUNT of attention has 

been devoted to the study of the position of the 
centre of shear for straight beams but there appears to 
have been little or no consideration of this position for 
beams of which the centre line is either slightly or 
heavily curved. The reasons are probably twofold. 
Firstly, no demand for such an investigation has been 
made by designers, and secondly, and probably not 
least, any very general inquiry into this problem will 
lead to complex expressions which would not only be 
diffiult and tedious to evaluate but will also be of 
uncertain accuracy. 


This problem of taking beam curvature into account 
starts from the point where it normally ends for most 
mechanical engineers, who usually require to estimate 
the bending stresses in members of large curvature such 
as hooks and chain-links. Now the classical theory of 
elasticity applied to these latter members has yet to be 
seen to conform with experimental results to such a high 
degree as makes the added labour attending such 
methods worth while. It has been found that “ engineers’ 
theory” applied to highly curved bars, i.e. such as 
Winkler’s Theory, when applied to these particular 
cases, provides results that are quite satisfactory, and it 
is on the knowledge of this that the approach in 
Section 2 is justified. No claim is made for an exact 
prediction of the results of possible experimental 
investigation, but it is believed that both the trend of 
the theeoretical results and their numerical values will 
provide a useful indication of what may be expected. 


Two simple thin-walled sections are considered, the 
channel and the semi-circle. The thickness of the 
material is assumed to be negligible in comparison with 
other dimensions and the influence of possible altera- 
tions in the shape of the sections under the action of the 
applied loads is also neglected, as is the influence of 
lateral pressure between fibres in the case of the semi- 
circular section. From these and other less important 
assumptions, expressions are derived for the position of 
the centre of shear when a moment acts in a plane 
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parallel to the centre line of the member and when 
normal and shear forces are also acting across the 
section. The particular cases of a plane vertical 
quadrant encastré at one end and carrying a simple 
vertical concentrated load at the other is then investi- 
gated for each type of section and a comparison of the 
results is made. 


NOTATION 
M_ moment on section about the point where the 
centre-line cuts it 
H_ normal force on section 
V shear force across section 
Pp normal stress on section due to bending 
tensile or compressive bending stress 
q_ shear stress 
t thickness of material 
cross-sectional area 
Ak? “modified area” of section 
e shear centre distance 


1. CHANNEL SECTION 
1.1. General Case 


Consider a channel cross section having its flanges 
unequal in width and thickness, as shown in Fig. 1, and 
a centre-line which is the arc of a circle of radius R, see 
Fig. 2. Across a section such as AB the moment and 
forces are as indicated in Fig. 3. 

For force equilibrium normal to the section, neglect- 
ing the bending stresses in the web 


p, at,-p.bt,=H . (1) 
and for moment equilibrium 
h(p, at, +p.bt.)=M ; . (2) 
Substituting from (1) in (2) gives 
M 
2p.0t,= -H 
and thus 
2hbt hH 3) 
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| p+2P 36 
CENTRE P 00 p 
LINE 
h 
Ficure 4. R, 
Ficure 1. Channel section. 
Consider next an element from the flange of width x, H 
see Fig. 1, whose length is R50, compare Fig. 2. For 
force equilibrium, see Fig. 4, a shear force on the back inna M 
face of the block is required, given by ; i 
op 
36 Axt = qt06R,, B 
; The moments of the complementary shear stresses about 
or q= x Op P, which are equidistant from the flanges, give 
R,, 06 
a 
In particular M.=k, = ht, and 
x Op. 
: . (4a) 
and The transverse shear force V is required to act at a 
x Op, distance e from P such that 


The complementary shear stresses to these quantities 
give rise to the ‘shear flow’ round the section. 


Utilising (3), 
x 
FR “hot, — hH)= xk, say (Sa) 
and 
= = (M+hH)= xk, say (Sb) 


FIGURE 2. 


Ve=M,+M.= = +k.b*t,}. 


> 


On substituting and rearranging, 


e= ql Ge *x 


1.2. A Particular Case 


In particular the case shown in Fig. 5 is investigated 
when a=b, and for the loading shown. Here, 


b 6) 


CENTRE 
LINE 


FiGureE 5. 


f 
4 
4 
/a\ 
‘\\s0 4/ 
Ww 
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2-0 | 
| | 
| | | 
| 
| | 
° 4 
% | 
0 os o6 0-4 0-3 
| 
eM | \ 
and 36 =WR sin 6 
H=Wcosé and Wsiné 
V =W 
Substituting in (6) we eventually obtain on reduction, 
that orl o2 03 0-4 0-5 O6 
: 
Ficure 7. 
where Thus 
When x=1,a straight beam is being considered and then rt dé AR 06" A 06 * Ak* 06 ' 
e=a/2=e,, say. Fig. 6 shows the variation of e/e, ) ROM = 
with x. Alternatively in terms of R and h, if h/R=v, 
(ME tan 9/2) (9) 
(7) ~ AK? a’) 
Now, also by considerations similar to those applied 
The variation of e/e, with v is shown in Fig. 7. in connection with Fig. 4, and except that forces and | 
not stresses have to be considered immediately, the shear 
2. SEMI-CIRCULAR SECTION stress across section OQ is given by i 
2.1. Figure 8 shows the section of radius r. Again dF =qt(R +) d6 
the curvature of the member is assumed to be constant 
along its length. Now it can be shown from ordinary or i< (R+y) (10) 7 
engineering theory’ that rt dé r 
y H t 
The total force on that portion of the section covered é 
by angle ¢ is F and 
F= ( pt+do 
0 Ficure 8. I 
1(M Semi-circular section. 
R+rcos¢ R | 
(= M _H _MR 2 
Ak? ! 


| 
Fy 
3 
ae 
= 
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Equating (9) and (10) 


= A(R +rcosq) 


(R?—P) 


The centre of shear is now given by 


ap 


aM 
| 


where 


e= do, . (12) 


Ns 


where e is measured from the centre of the semi-circle. 
Substituting from (11) in (12) and simplifying, 


___2R_ 0M 
h? /(R?—r?) 
where - | r ” 
Ros? 
and 1,= | r r 2 
R av(i-(4) ) 


Writing r°t.k(v) for Ak*®, where v=r/R, substituting 
for J, and simplifying further, 


lv og 
(13) 
It may be shown that, since 
2_ RY 
Ak? = dA 
1] = 
then - =k(v) 


and, if v is small, 


If r is finite and R — © then v —> 0 and Ak? —>=/2r'°t, 
ie. Ak® becomes the second moment of area of the 
section about the centre-line. 


2.2. Applying the result of (13) to the particular 
Situation treated above for the channel section, shown 
in Fig. 5, on reduction, 


0-2 


0-2 0:3 0-4 0-5 


FIGURE 9. 


» is shown in Fig. 9. It 


may be noted that when v — 0, Say. 


The variation of e/r with wv 


3. CONCLUSION 
In Fig. 9 a comparison is also made of the variation 


of e/r with ~ for both sections, when / and a for the 
channel sections are each equal to r. Two points 
emerge: 


(i) that in both cases the tendency is for the centre 
of shear to move from the outside to the inside 
of the section as the radius of curvature 
decreases: 


and (ii) the rate of this movement in the case of the 
semi-circular section is much larger than in the 
case of the channel section. 
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Buckling of Rectangular Plates under Bi-Axial Compression 


J. S. PRZEMIENIECKI, B.Sc., Grad.R.Ae.S. 
(Structures Department, The Bristol Aeroplane Co. Ltd.) 


HE problem of the buckling under bi-axial com- 

pression is considered for flat rectangular isotropic 
plates with simple edge conditions and no lateral 
restraint. The buckling stress coefficient is plotted 
against the side ratio for various conditions of edge 
restraint and a known compression or tension on the 
sides of the plate. It is found that there is some reduc- 
tion in the buckling stress if the sides of the plate are 
subjected to compressive stresses and, conversely, there 
is an increase for tensile stresses. Furthermore, for 
plates with large side ratios, there is a rapid decrease of 
the spanwise buckling stress as the chordwise com- 
pressive stress approaches its appropriate Euler 
instability value*. 


INTRODUCTION 

The analysis of rectangular plates under bi-axial 
compression presented by Wittrick''’ shows that for 
plates with sides simply-supported it is still possible to 
use the standard curves for uni-axial compression *’, 
provided the general formula for the buckling stress is 
modified to include the compressive or tensile stress on 
the sides of the plate. For sides built-in, however, the 
use of uni-axial buckling stress data is not possible. The 
analysis presented here is basically similar, except that 
the buckling stress coefficient is given in a form directly 
applicable for design purposes. In addition, the results 
for built-in edges on the sides of the plate are also 
derived. 

If a rectangular plate is subjected to a uni-axial 
compressive stress f,, it will buckle when f, reaches a 
critical value given by Refs. 2 and 4 


f,=KE(t/bP . (1) 


where K, the buckling stress coefficient, is a function of 
the side ratio a/b and the conditions of edge restraint. 
A similar expression can be found for bi-axial com- 
pressive stresses f, and f,, except that it is more con- 
venient to introduce the buckling stress as 


where f,=47°E’(t/b)? and k is the buckling stress 


coefficient which in this case is a function of the edge 
restraint, the side ratio a/b and the stress ratio f./f,. 


NOTATION 


x, y Cartesian co-ordinates, spanwise and 
chordwise directions respectively 


*Since this paper was written the author's attention has been 
drawn to the fact that results similar to those given here have 
also been published in an advanced restricted N.A.C.A. report 
(see Ref. 5). 
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f,. f. direct compressive stresses in the x. 
and y-directions respectively 
E Young’s modulus 


E’=E/(1-v’) 


Poisson’s ratio 


K, k_ buckling stress coefficients, equations 


(1) and (2) 


a_ plate length 

b plate width 

t plate thickness 

w lateral deflection 

Y function of y only, equation (4) 

m_ integer 

A, B, C, D_ constants of integration 
2,, %, parameters defined by equations (7) 

and (8) 


DERIVATION OF BASIC EQUATIONS 
Consider a flat rectangular isotropic plate subjected 
to compressive stresses f, and f, as shown in Fig. 1, 
where x and y are directions parallel to the sides and 
ends of the plate respectively. The ends of the plate 
are simply-supported while the sides may have various 
edge conditions. The compressive stresses are taken 
here as positive and the stress f, is assumed to be 
known; it may be either positive or negative (tension). 
The lateral deflection w of a rectangular plate under 
bi-axial compressive stresses f, and f, must satisfy the 
differential equation”:* 


+f, By? =0. 
(3) 


Assuming now that under the action of compressive 
stresses the plate buckles into m sinusoidal half-waves, 
the lateral deflection w may be taken as 


w=Ysinmzx/a. (4) 


where Y is a function of y only. Substituting equation 
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Figure 2. Buckling stress coefficient for flat plates under bi-axial 
compression: ends and sides simply-supported. 


(4) into (3) the function Y must satisfy the following 
differential equation 
(f2/f.)} Y2+ 
+(mz/ay (f,/f)} Y=O (5) 
where f,=47°E’ (t/b)* and the suffixes denote ordinary 
derivatives with respect to y. The general solution of 
this equation is 
Y=Acosh z,7y/b+Bcosh z,7y/b+ 
+Csinhz,7y/b+Dsinhz,7y/b . (6) 
where 
a,°=(mb/ay’ —2 (f2/fo)+2 {(mb/a)? + 
a,*=(mb/a)’ —2 (f2/fo)— 2 {(mb/a)? /fo—f2/ fo) + 
+(f2/f.)7}* (8) 


3-0 | 
f, 
20 
[ SIDES FIXED 
a 
J SIDES SIMPLY SUPPORTED 
TENSION [COMPRESSION] 


° O-4 0-6 o8 
e 


Ficure 4. Variation of buckling stress coefficient with S39, 
for 
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Ficure 3. Buckling stress coefficient for flat plates under 
bi-axial compression: ends simply-supported, sides fixed. 


and A, B, C, D are the arbitrary constants of integra- 
tion. The first two terms in equation (6) represent 
symmetrical deformation about the x-axis and the last 
two terms represent antisymmetrical deformation. 


ENDS AND SIDES SIMPLY-SUPPORTED 
Consider first symmetrical deformation about the 
x-axis so that the deflection w is given by 


w=(A cosh z,7y/b+Bcosh z,zy/b) sin mzx/a. 
(9) 
If the sides are simply-supported, then the correspond- 
ing boundary conditions at the edges become 
y=+b/2 
y=+b/2 


w= 0?w/dy?=0 at 
or at 
which gives 
Acosh2z,7/2+ Bcoshz,7/2 =0 
cosh 2,7/2+ 2z,*Bcosh z,z/2=0. 
The condition for non-zero values for A and B repre- 
sents the condition for buckling. This requires that the 


determinant formed by the coefficients with A and B is 
identically equal to zero. Hence 


cosh cosh z, 7/2 


2,7 cosh 2, 7/2 2,7 cosh 2, =/2 


or (z,”—2,°)cosh z, 7/2 cosh z,7/2=0 which, in virtue 
that z,?+2,?, becomes or z,7=-—1. Either 
condition reduces finally to the same result giving 


(f,/f.)+(a/mb) (f./f,)— 4 (a/mb + mb/ay’ =0 


which determines the buckling stress coefficient k =f, /f, 
in terms of f./f,. The function (10) is represented in 
Fig. 2 as a system of curves relating f,/f, to a/b for a 
series of values of f./f,. 

In can be shown that the anti-symmetrical deforma- 
tions produce higher buckling stresses and therefore can 
be neglected in the present analysis. 
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ENDS SIMPLY-SUPPORTED, SIDES FIXED 
The appropriate boundary conditions in this case are 


w=0w/dy=0 at 
or Y=Y,=0 at 


Considering again only symmetrical modes of buckling, 
the condition for non-zero values for A and B becomes 


cosha,z/2 | 


cosh 2, 7/2 | =o. 


| a, sinh 2,7/2 a, sinh a, 7/2 
Hence the equation for the buckling stress coefficient 
f,/f, becomes simply 

(11) 


The function (11) is represented in Fig. 3 as a system of 
curves relating f,/f, to a/b for a series of values of f../f,. 


a, tanh 2, 7/2 — 2, tanh z,7/2=0. 


CONCLUSIONS 

The curves relating f,/f, to a/b for various values of 
the chordwise stress parameter f,/f, show that there is 
some reduction in the buckling stress coefficient for 
positive (compressive) values of f,/f, and that there is 
an increase in the buckling coefficient for negative 
(tensile) values of f./f,. 

For convenience the buckling stress coefficient f,/f, 
is also plotted against the stress parameter f,/f, for 
(a/b)=00 (see Fig. 4). This shows that there is a rapid 
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decrease in the buckling coefficient as f, approaches jts 
appropriate Euler instability stress 


f.=0-25f, for sides simply-supported 
and f.=1-0 f, for sides fixed. 


For a given value of a/b and f./f, the number of 
half waves in which a plate buckles along its length js 
given by a number of troughs to the left of the appro. 
priate side ratio a/b; part of a trough being counted as} ‘0 
a complete one. 
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A Chart for Steady Circling Flight 


by 
Dr-Ing E. SPONDER 
(Madras Institute of Technology) 


T OR horizontal, or nearly horizontal, steady straight 

flight there exist some charts from which can be 
determined the speed at any altitude, the inclination of 
the flight path, and other quantities. This can be 
readily done, for instance, by a graphical comparison of 
the forces of drag and thrust, or by a corresponding 
representation of the powers, required and available. 
These charts, however, always involve the polar of the 
particular aeroplane, and are thus valid only for that 
aeroplane; they must therefore be replotted for each 
example. To avoid this and obtain a generally valid 
chart, the effect of the aerodynamic characteristics of the 
particular aeroplane must be eliminated; it is thus 
possible to prepare a chart even for the steady circling 
flight at any inclination of the flight path. 

The three equations for the description of the flight 
path may be derived on the simple assumption that the 
thrust acts almost along the direction of flight, i.e., that 
the angles of attack and of sideslip are negligibly small. 


NOTATION 
W weight of the aeroplane (Ib.) 
g specific gravitational force (= 32:2 Ib./slug at 
sea level, corresponding to 32:2 ft./sec.*) 
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S wing area (or other reference area, ft.’) 
T total thrust (Ib.) 
T, see Fig. 1 
V speed of flight (ft. /sec.) 
p density of air (slug/ft.*) 
p, density of air at sea level under Standard con- 


ditions ( = 70 slug / ft. ) 
h_ altitude of flight (thousands of ft.) 


o=, relative density of air 


0 


66—h 
+. valid up to 36,000 ft] 


q= a dynamic pressure of flight (Ib./ft.*) 
y angle of inclination of the flight path, positive 


in climb (°) ral 

® angle of bank to right or to left (°) aly 
T: 


+This approximate formula for « is derived from an essay of 
Knoller in Flugsport, 1924, p. 309, and reproduces the 
Standard Atmosphere up to the altitude mentioned with less 


— , valid down te 


than 1 per cent. error; its inverse is h~66 


o—0°3. 
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+b / 2 
y=x 
+b 
y=x 
| 
| 
ain “4 
4 
Wee, 
} 
; 
| 
° 


TECHNICAL NOTES—E. SPONDER 


= 
—— 


oplane 
paper, ¢ 
1dvice, 


Some 
Plates 
> Aero- 


2s Data 


Stress 
Royal ¢ 
134-638, 


s. First 
940, 

nations 
yr Flat 


1 con- 


0 ft) 


sitive 


say of 
>s the 
th less 


ywn te 


T 
oT 
Ficure 1. Substitution of T by T, + 4 
R_ radius of curvature of the projection of the 
flight path on the horizontal plane (ft.) 
angular velocity of the radius vector R 
(rad. /sec.) 
t= = . time required for one complete turn (sec.) 
C= lift coefficient of the aeroplane 
C= ert . drag coefficient of the aeroplane 


Sq 
oT 

C= a / S , slope coefficient of the curve of T 

against q 


C 
r= over-all lift/drag ratio 


/Q- rs) 


n= cua . load factor 


The three reference axes fixed in space, for which is 


set up the equilibrium of forces in steady circling flight, 
are: 


the direction of flight itself; 
a horizontal direction to the right perpendicular 
to the direction of flight; 
the upward direction perpendicular to the two 
others; the lift makes with it the angle of 
bank 9. 
Corresponding to these three directions, the follow- 
ing three equations result 


C)Sq= T —W siny, (1) 
C,Sq sin ¢= = QV cos + 
8g Q) 
or = cos y)* 
gR 
C,Sqcos¢=Weosy. . (3) 


In (1), T is strongly dependent on q; but for the 
range of speeds of flight in question one can nearly 
always apply with sufficient accuracy the approximation 
T=T,+(0T/¢q)q (Fig. 1), in which T, as well as 
eT /@q are constants. 

In addition, introducing C;=(@T /0q)/S as the slope 
coefficient of the curve of T against g, (1) assumes the 
more convenient form 


(Cy —Cy) Sq=T,,— W sin y. (4) 


It follows from this that one may replace the variable 
T with a constant T,, provided that one diminishes the 
Cy-values by Cy, therefore, according to the positive or 
negative sign of Cy, by simply displacing the polar to 
the left or to the right. An over-all lift/drag ratio 
r,—C,/Cy is then, for the displaced polar, modified to 


= ae =r, [(1- 


It is assumed for the following that the lift coefficient 
C,, and the corresponding over-all lift/drag ratio r, are 
known, regardless of the conditions which govern their 
choice, be it the question of the fastest and the sharpest 
turns* or any other one. It is then the problem to bring 
the relations (1), (2) and (3) into such a form that they 
can be applied in a generally valid manner—i.e. without 
being bound to any particular polar—to make a quick 
analysis of the steady circling flight at that C,- and 
r,-values. For that purpose, dividing (1) by (3): 


(2 ~siny) | cosy. 


The essential variables here are r,, 9, T/W and y, 
the relationship between which is represented in Fig. 2. 

Proceeding from the y-scale at the left, one finds at 
the inter-section of the T/W-curve in question as a sub- 
sidiary result the value of r,cos@ on the scale at the 
bottom; more important, by moving along the ordinate 
through this point to the given value of r, on the right- 
hand scale, one can interpolate the associated value of 
» from the family of ¢-curves. These curves are all 
parallel, which facilitates the interpolation; besides, 
they are so arranged that the load factor n can be found 
with their help, as one has from (3): 


C.Sq 
Ww 


(5) 


(6) 


This simple relation permits a convenient graphical 
determination: one need only follow the same ¢-curve 
to the intersection with y and read off then directly under 
it the value of n, as indicated in Fig. 2 for the 
subsequent example. For the inverse problem change 
suitably the marked directions of the arrows. 

The only aerodynamic characteristic of the aero- 
plane involved in this chart is the over-all lift/drag ratio 
r,. Which of course must be known but may have any 


*To determine this, it is necessary to know also the nature of the 
polar between C,, for ry aNd Cy Which region is not 
covered by the usual parabolic formulae. It is found in many 
cases that the entire polar of the aeroplane can be represented 
by a quarter of an ellipse, which can be conveniently trans- 
formed into a quarter of a circle by suitably proportioning 
the scales for C, and C,. Given only C,, (at C,=0) and 
Ta max: and using a judicious estimate of C,,,,,- it is possible 
by choosing the ratio of the units for the C,- and C,-scales 


amax 
“circular polar,” which forms a simple, sufficiently accurate, 
and reasonable basis for numerical computations (see, for 
instance, Fig. 3). 


, to determine completely this 
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value. As special cases, both steady straight flight QV 
(»=0) and gliding flight (7 /W =0) are included; its chief © otang 
purpose is, however, the general treatment of steady r . (7) 
circling flight, especially of the fastest and sharpest or = ee. hence R= ey 


turns. In this connection it is important to note that if 
a polar corrected for Mach number is to be used, the 
corrections should be based on the actual speed of flight 
round the curve. 

As to the scales of the chart, it may be remarked that 
one must have r, cos y=n(r, cos ¢)=constant; here, the 
constant is 5, but any suitable value can be chosen. 

Now, the speed of flight V=/(2q/c) is deter- 
mined by equation (6) and may serve for a verification 
and, if necessary, for a correction of the total thrust T 
introduced in Fig. 2. Finally, by (1) or (4) too, the result 
can be checked. The speed of climb during the circling 
flight is then V sin y. 

The questions of the time required for one complete 
turn ¢, and the radius of curvature R of the projection 


The best results with regard to the quickness of 
manoeuvre in circling flight will generally be obtained 
at over-all lift/drag ratios which are only a little smaller | 
than Famax- Thus, with the knowledge of Fraimax. one! 
would have already a very useful indication of the} 
values of r, that must be introduced in the chart. It j 
possible to give average values of amax for the various | 
types of aeroplanes, which may be of good assistance 
in the absence of better data: Table I* gives th? 
necessary information. 

As an example, the performances of a jet fighter 
which executes its fastest and sharpest turns in climbing 


*Chiefly after Jane's All The World’s Aircraft 1950-1951 and 


of the flight path are treated by dividing (2) by (3): 1952-1953, 
= 30 
20 
410 
\ \\ 49 
\ 
\ Ficure 2. Chart for steady 
= | 
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/\ / / / \ 
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TECHNICAL NOTES—E. SPONDER 
flight with y= 10° at 13,200 ft. (e=2/3) may now be 4°5 
considered. The constructional data available are: oe 
W =11,300 Ib., Lmax 
the power plant is a Rolls-Royce Nene II, 1947. C 
TABLE I 
Over-all lift/ drag ratio 40 
7? Ta max (average value) Rmin (sharpest turn) 

Single-seat fighter 13 
Four-engined transport aeroplane | 15 (fa atest turn) 
Sporting and touring aeroplane a 

landing gear retractable 13 

landing gear not retractable | 11 
Sailplane up to 32+ 05 

As no data about the horizontal-flight performance 

are accessible, the polar can be only guessed at. The 
Cp,-value (at C;,=0) may be about 0-021; a large error 
here would of course lead to very doubtful results for 

e maximum speed, but would be of much less 0 = 
0 Cp, 0-4 0-2 C, 03 


importance in the calculations for circling flight. With 
remax 13 (for single-seat fighters) and a likewise esti- 
mated Cy, max= 1:4, one has all the necessary data for 


TABLE II 
Flight Fastest "Sharpest computed 
characteristics turn turn from 
| 
C, 0:77 0-90 any method to 
10-3 9-67 find these 
optimum values 
41° 36° Fig. 2 
n 1:30 P22 [or equation (5)] 
q 70:5 Ib. /ft.? 56:6 Ib. /ft.* 
V 298 ft./sec. 267 ft./sec. 
=203 m.p.h. = 182 m.p.h. equation (6) 
V sin y 51°7 ft./sec. 46-4 ft./sec. 
= 3,100 ft./min. =2,780 ft. /min. 
R 3,150 ft. 2.980 ft. 
(minimum) equations (7) 
67°3 sec. sec. 


(minimum) 


+According to Aero-Revue S7. 


FicureE 3. The probable polar of the jet fighter. 


plotting a circular polar as shown in Fig. 3. 

From the characteristic curves of the jet engine, the 
relation T=3,390—2:7q may be taken which covers 
fairly well the range of the expected speeds of flight. 


Hence 7, =3,390 Ib., 
aT 
and S=Cy=-0-:010; =0-30. 


The chart for steady circling flight can be used, as 
indicated by the arrows marked for the sharpest turn 
in Fig. 2; the results are given in Table II. 

It can be seen that the results for the fastest and 
sharpest turns differ only a little from each other. 
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Effect of Variations in weight or Normal Acceleration on the Ceiling 
of Jet-Propelled Aircraft — Correspondence 


R. J. WOODHAMS* 

HE TECHNICAL NOTE by Mr. P. L. Sutcliffe in 

the June JoURNAL (p. 435), apart from repeating 
well-known proofs that T/o is constant above the 
tropopause and for loss of height due to g in turns, seeks 
to show that rocket boost is necessary for high altitude 
fighters. 

The Note is based on the false assumption that the 
minimum drag speeds, or speed at 1g ceiling, for 


' bomber and fighter are identical, whereas in fact such 


*Received 27th June 1955. 


cannot be the case due to the much higher wing area of 
the bomber. Further (and from the Note this holds at 
all altitudes) the fighter must have a high rate of climb 
at the bomber |g altitude in order to get there in a 
reasonable time. Normally this well covers the power 
required for turns. 


By and large, one can say that, above the tropopause 
as the bomber 1g ceiling height increases, the ratio of 
excess thrust required for climb or turns to sea level 
rated thrust for the fighter decreases. Given the 
necessary sea level thrust rating to ensure rapid climb to 
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bomber |g ceiling, the fighter design problem would which follows from this calculation is of course 
become easier with increasing bomber altitude. The 14,500 ft. below the | g absolute ceiling). 
fighter thrust/weight ratio at sea level will of course Not only are the maximum altitude reductions 
increase with increasing bomber altitude. This during interception manoeuvres likely to be less than 4 
increased thrust/weight ratio of fighter relative to those implied by (7a), but also these reductions are f 
bomber is required at all altitudes in order to produce achieved theoretically only after an infinite time jp 
the necessary minimum time of climb to bomber ceiling. manoeuvre. The rate of descent is a function of the 
thrust-drag unbalance which tends to zero at the new 
“ steady-state” altitude. Practical problems concerning } 
_ The derivation of the expressions in the Note is interception manoeuvres do not show height losses as |_| 
independent of the actual minimum drag speed of the large as those obtained from (7a). I 
aircraft under consideration. The significant relation- Figure 2 refers to the effect of weight change on the { : 
ship is that given in equation (5), and it is worth repeat- absolute ceiling and it will be obvious from the fore. |, 
ing that this is quite independent of the shape of the going that it should not be interpreted as applying also 4 ¢ 
aircraft, the thrust of its engines, its maximum altitude to cruising ceilings. The optimum design condition on 
at a given value of n and W, or its flight Mach number. a long-range aircraft such as a bomber, can be atan L/D | | 
With these facts in mind it is quite possible to com- away from maximum, even when speed-limited. Weight | | 
pare, for instance, the supersonic fighter and a subsonic changes from such a design condition have an effec |! 
bomber on exactly the same basis, the relationship much smaller than that indicated in Fig. 2. Q 
between the various ceilings being taken directly from t 
the graphs given in the Note. P. L. SUTCLIFFE a 
is he aaa It should be emphasised from the outset that the 
iShuet-Brechers anil Maslead Led) Note only considers flight in longitudinal equilibrium; | 
Attention should be drawn to the fact that the Tech- climbing and diving flight or flight with longitudinal}, 
nical Note by Mr. P. L. Sutcliffe applies to a particular acceleration are not covered by the equations. oe 
set of conditions. The purpose of this letter is to _ The maximum altitudes, whether 1g or otherwise, } 
emphasise this point because the findings are not gen- discussed in the Note are related to a given Mach num. 4 f 
eral, and therefore are liable to misinterpretation. ber and are not necessarily the absolute ceilings defined {| * 
The simplifying assumptions in question are that by the maximum lift/drag ratio and, as quoted in reply | | 
L/D is constant during a manoeuvre or with change of to Mr. Woodhams, the relationship given in equation (5) : 
weight, and that initial conditions apply to the absolute is quite independent of the actual minimum drag speed. : 
lg ceiling. These assumptions restrict the aircraft to The validity of the assumption of constant lift/drag } 
operating throughout at the lift coefficient for maximum ratio when flying at a given Mach number and engine > 
L/D. setting (T/«r), irrespective of altitude, can be shown as}, 
In the general case, (4) should read follows:— 
For the equilibrium condition of thrust equalling drag 
so that (5) becomes: T/o=Cy}p 
te, (5a) where a=speed of sound in the stratosphere and the} 
nW, Po, other symbols are defined in the Note. P 
Now if (T/«) and M are fixed it is clear that, for any 
so that (7) may be written particular aircraft, C, must be completely specified. | , 
n,W, This fixes the value of the lift coefficient and conse-) 
h, —h, = 20,784 [ tos, ~ log. P ] (7a) quently the value of the lift/drag ratio, irrespective off « 
the altitude. e 
Using (7a) we see that the correction to (7) is This being so, Mr. Brailsford’s equation reduces to} t 
Ah= 20,784 log. P that of the Note, i. 
= 4,630 ft. for P=1-25 ( nW = e 
constant 
= 8,450 ft. for P=1-5 
= 14,390 ft. for P=2:-0. ‘ (8) and in his terminology P=1-0, Log.P=0, Ah=0. It d 
The effect of (8) can be shown readily for an inter- should be noted that values of P other than unity} » 
ception at constant Mach No., involving a 2g man- indicate violation of the conditions of M and (T/«) con- tt 
oeuvre, where the terminal condition is at the C;, for stant and steady level flight. ; 
maximum L/D. In this example, P= 1-25, so that the Furthermore, in the case of the long range aircraft, : 
maximum reduction in altitude is 14,500—4,630— the optimum flight plan for either range or altitude) 
9.870 ft. and not 14,500 ft. as might be inferred from requires both the engine setting and the Mach number} 
Mr. Sutcliffe’s Note. (The “steady-state” altitude (and inherently the lift/drag ratio) to be fixed. Thus} 
= eee the effect of weight changes during a ceiling cruise are}; 
*Received 7th July 1955. correctly given by the equations of the Note. vi 
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Reviews 


GAS TURBINES FOR AIRCRAFT. /van H. Driggs and Otis 
E. Lancaster, Ronald Press Co.: New York. 349 pp. Illus- 
trated. $10.00. 

This book, by the Chief Scientist, U.S. Naval Develop- 
ment Centre, and an Assistant Director, Bureau of Aero- 
nautics, raises expectations which are, on balance, well 
realised. The inclusive title is perhaps rather misleading, 
for until the reader realises that the main purpose of the 
early chapters is to provide background and introduction to 
Chapters 8, 9 and 10, he may feel some parts of them to be 
sketchy and inadequate. The main justification for the 
book lies in fact in these later chapters, which deal with 
Gas Turbine Performance, Control Problems, and Aircraft 
Performance Analysis. In these chapters one feels that the 
authors are on their own ground, and it is on them that the 
book should be judged. They might well have been 
amplified, even at the expense of more cursory treatment 
earlier. They would justify, in the reviewer’s opinion, the 
more suitable title of ‘“‘ Aircraft Gas Turbine Performance.” 

The earlier chapters deal with Perspective (a useful 
historical survey), Fundamentals, Cycle Analysis, Com- 
pressible Flow, Compressors, Turbines, and Combustion, 
and are on the whole well done. The authors’ treatment is 
fundamental, and directed to the “ why” rather than the 
“how.” The reviewer has no quarrel with this in principle, 
but confesses (is this insular prejudice?) to a preference for 
evaluating cycle efficiencies and output on the “traditional” 
basis of maximum temperature, component efficiencies, and 
mean specific heat, rather than through changes in enthalpy 
and entropy. Such a treatment may lend itself to the rapid 
evaluation of ideal cycles (using the charts given in the 
book) but it is apt to overlook the essential cycle limitations 
with which the designer is concerned. 

The description of flow phenomenon in diffusers and 
nozzles is lucid and worthwhile. The principles of centri- 
fugal and axial compressors are discussed briefly, although 
in view of the extensive use which is made later of 
“ dimensionless * performance parameters, it is surprising 
that this method of plotting should receive little more than 
passing reference. The principles of turbines are covered 
in ten pages, after which one feels that the authors turn 
with relief to a review of blade and disc materials and 
cooling methods. Chapter 7 on Combustion is rather super- 
ficial; the terms “* combustion chamber,” “ combustor ” and 
“burner” are used variously and interchangeably. How- 
ever, in view of the consistent use of enthalpy rather than 
temperature throughout the book, it is refreshing to find 
here the admission that: ‘“ Unfortunately the strength of 
the metal parts responds to temperatures rather than 
enthalpies 

Chapters 8, 9 and 10 are really worth while and repay 
careful reading. One might only have asked for more 
detailed treatment. (For example, the discussion of what 
we should call equilibrium running conditions for the 
turbo-jet is rather over-simplified by the assumption of an 
infinitely-variable propelling nozzle, and the uniqueness of 
the engine running line (with a choked fixed nozzle) thereby 
missed.) The basic problems of turbo-jet control are well 
illustrated by a number of proposed operational schedules, 
and acceleration limits are clearly if briefly discussed. 
Methods of estimating engine weight are described and 
formulae derived for total engine and fuel weight for 
various operating conditions. In Chapter 10 the authors 


describe a novel and rapid method devised by them (and 
used by the Bureau of Aeronautics) for predicting the 
complete aircraft performance on the basis of generalised 
“thrust required” curves. The advanced student will 
undoubtedly find these chapters stimulating and valuable. 
The book is admirably printed, illustrated and produced.— 
E. GLAISTER. 


FLUID DYNAMICS OF JETS. Shih-I Pai. D. Van Nostrand 
Co. Inc., New York, 1954. 227 pp. Illustrated. 37s. 6d. 

The first impression given by the title of this book is 
that jet flow cannot be satisfactorily treated in isolation, 
as it is only a particular section of general fluid dynamics. 
If this view is correct it follows that many parts of the 
general subject, for example, conformal transformation, 
have to be cursorily described in order to explain the 
application of the general theory to the flow of a jet. The 
author of the present book has failed to resolve this diffi- 
culty, because in fact no solution is possible within the 
chosen framework. Nevertheless Dr. Pai has written a 
valuable book which collects a lot of information on recent 
theoretical investigations of jet flow; experimental results 
are given in some cases where a comparison with theory 
is possible, but the emphasis is on the theories. 

The first three chapters deal with inviscid flow at 
subsonic and supersonic speeds: this ground has been 
adequately covered in general works on fluid mechanics. 
The next four chapters deal with laminar and turbulent 
jet flow, including the case of the mixing of two streams 
of different velocities. Wakes which, as the author points 
out, are similar in many respects to jets, are, however, not 
discussed in these chapters. The temperature distribution 


* in heated jets is not thoroughly discussed and this omission 


is regrettable. An account of the mixing of streams and 
jets of different gases is included and this has not been 
presented previously. 

The last two chapters deal with the stability of jet 
flows of incompressible and of compressible fluid. The 
author has done research on this subject himself, and he 
has given an up-to-date, though unavoidably incomplete, 
account of it. It would have been useful if he had indicated 
the conditions under which laminar jets and stability are 
important: any jet which is not fully turbulent is, in fact, 
a rare occurrence. 

Dr. Pai writes clearly and has given us a useful book, 
but it is a book which will be mainly used as a reference 
work and is not suitable for use as a text book.—H.B:S. 


AIRPLANES OF THE WORLD. Douglas Rolfe and Alexis 
Dawydoff. Simon and Schuster Inc., New York, 1954. 320 
pp. Illustrated. $2.95. 

In five years time, when my small son is old enough to 
view an Avro 504 as one of the great aeroplanes of history, 
instead of as an old-fashioned stick-and-string monstrosity, 
this is the book I shall take from the shelf for an evening’s 
reminiscing. 

It contains drawings by Douglas Rolfe of 843 aircraft 
spanning five centuries, from Leonardo da Vinci’s man- 
powered ornithopter to the delta-wing Avro Vulcan and 
supersonic research aircraft of today. The selection is first- 
class, covering almost all of the aeroplanes one looks for 
in a historical review, and auite a lot of the lesser known, 
but no less interesting, experimental types. 
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Most notable omissions are Ellehammer’s pioneer air- 
craft of the 1904-7 era, Fabre’s seaplane of 1910, the 
Heinkel He 178, and Russian types such as Mozhaisky’s 
1882 monoplane—claimed by Moscow as the first powered 
aeroplane to fly—and the modern MIG-15 and IL-28. 

The style of the drawings varies considerably, giving the 
impression that they were produced over a period of several 
years for different magazines, and then gathered together 
for publication in this book. A few of them will offend the 
enthusiasts who insist on perfect technical accuracy; but, 
particularly in the pre-1914 era, they display the aircrafts’ 
“character ” more clearly than would a faded or heavily- 
retouched photograph. 

Each of the nine sections of the book is introduced by 
a brief review of aviation development in the period 
covered by the drawings that follow. Written by Alexis 
Dawydoff, this narrative is readable and far more compre- 
hensive than one would expect in a total of only 10,000 
words. It contains some rather remarkable Americanisms, 
and a number of mis-spellings and factual errors. For 
example, Sperry invented the first automatic pilot in 1912, 
not 1932; the R.F.C. was established in 1912, not 1911; the 
designer of the Antoinette was Levavasseur, not Levasseur; 
the “ Bloater” was the B.E.8, not B.E.3; the Pfalz D.III 
was more, not less, rugged than the Albatros D.V; and the 
three-seater of 1913 was not the first Sopwith tractor design. 

These are mere pin-pricks in a book that will enable 
many an old-timer to relax in an armchair, turn its pages 
slowly and recall the days when flying really was flying.— 
JOHN W. R. TAYLOR. 


STROMUNG DURCH ROHRE; UMSTROMUNG VON 
KROPEN BEI ZWEIDIMENSIONALER STROMUNG; 
UMSTROMUNG VON KORPEN BEI RAUMLICHER 
STROMUNG. (Flow through tubes; flow past two-dimensional 
bodies; flow past three-dimensional bodies.) Schiller. 
Springer-Verlag, Berlin. 134 pp. Diagrams. In German. 

This is a reprint of a section of the first part of Vol. I 
of the Randolt-Bornstein Numerical Values and Functions 
for Physics, Chemistry, Astronomy, Geophysics and 
Technology. Tables and curves together with numerous 
references are given for the aercdynamic properties of pipes 
and bends, radiators, flat plates, surface excrescences, aero- 
foils and wings, bodies of revolution, undercarriages and 
the numerous bits and pieces associated with aircraft. This 
volume forms a useful companion to the well known work 
of Hoerner (* Aerodynamic Drag”) although it is to be 
regretted that so much data is given with so little discussion. 
—J.A.D. 


JAHRBUCH 1953 DER WISSENSCHAFTLICHEN GESELL- 
SCHAFT FUR LUFTFAHRT. Edited by H. Blenk. Vieweg 
and Sohn, Braunschweig, 1954, 224 pp. Illustrated. DM.28. 
In German. 

Previous volumes of Jahrbiicher in the library cover the 
years 1921-1926 and 1929 and while the style may have 


THE OBSERVER’S BOOK OF WEATHER. Reginald mM. 
Lester. Frederick Warne, 1955. 151 pp. Illustrated. Ss. 
This small book gives an interesting description of the 


main weather phenomena in terms which can be readily ¥ 


understood by anyone who is interested in the subject. It 
is not, neither is it intended to be, a book for the expert, 
although even he may find here a few simple rules to 
enable him to fix more accurately the dates for his picnics, 
The colour plates add to the beauty of the volume but are 
not always particularly clear.—Dp.c.s 


THE CONTINUUM AND OTHER TYPES OF SERIAL 
ORDER. E. V. Huntington, Second edition. 1955. 80 pp, 
1 dollar. 


ALMOST PERIODIC FUNCTIONS. 4. S. Besicovitch. 1954, 
180 pp. 1 dollar 75. 


TREATISE ON THERMODYNAMICS. Max Planck (Trans- 
lated by Alexander Ogg) Third Ed. 1955. 297 pp. 1 dollar 75, 


ADVANCED DYNAMICS OF A SYSTEM OF RIGID 
BODIES. E. J. Routh. Sixth Ed. 1955. 484 pp. 1 dollar 95, 


THE DYNAMICAL THEORY OF GASES. Sir James Jeans, 
Fourth Ed. Dover Publications, New York, 1955. 444 pp. 
2 dollars. 


The publication dates given are of the special reprints 
that Dover Publications have prepared of these old classics. 
The first was originally produced by Harvard in 1917 and 
the second by C.U.P. in 1932, and Dover Publications are 
performing a great service by reprinting these titles which 
students must have spent much time seeking in secondhand 
bookshops. None of the books needs reviewing in the 
usual sense because each is already well known. It is 
hoped that this firm will continue the good work and this 
reviewer would like to suggest Todhunter and Pearson 
History of the Theory of Elasticity as a future title — 
F.H.S. 


MANUFACTURE AND APPLICATION OF LUBRICAT- 
ING GREASES. C. J. Boner, Reinhold Publishing Corpora- 
tion, New York; Chapman & Hall Ltd., London, 1954. 977 pp. 
Illustrated. 148s. 

During the past fifteen years or so, a vast amount of 
effort has been directed to the study of methods of manu- 
facture and application of lubricating greases. This has 
resulted in grease technology, in all its aspects, being 
placed on a more scientific basis and becoming more 
intelligible to the workers in the field. There was, how- 
ever, an urgent need for a reference book on the subject 
and the present volume amply fulfils this requirement. A 
very extensive survey of all available results has been made 
by the author (himself a prominent worker in the field) and 
an exhaustive bibliography is included. 

The book is stated to be for the use of “ manufacturers, 
research chemists, lubricating engineers, raw materials 


changed, the excellence of the papers is still what one has suppliers and all users of grease lubricants.” It will 
learnt to expect from Germany. The present Year Book probably be found that most benefit will be gained by 
covers the annual conference of W.G.L. held at Géttingen workers in the productive side of the industry, since this | 
in 1953 and is a collection of authoritative papers on all aspect is very well covered. The chapter on raw materials . 
aspects of aeronautics, but mostly on aerodynamics; one is very compact, yet provides the essential information. t 
is able to gauge their quality by such well-known names Manufacturing processes are dealt with for all the common § 
as Professor Young and Dr. Lachmann. There is a soap thickeners, including “mixed” and “complex” / 
summary of each paper in French and English. soaps, with comprehensive details of equipment used, and . 
A short while ago a daily paper’s air correspondent the latest types of non-soap thickeners are discussed. : 
expressed his disbelief in the possibility of a renaissance of The book will prove to be valuable also to research ; 


German aeronautics in, I think, 20 years. He should 


glance at this annual.—F.H.s. 


workers, particularly the new entrants into the field, who 
are now provided with an up-to-date picture of most of the 
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significant developments. For lubricating engineers and 
general users of greases, it is felt that the volume will 
prove to be somewhat less informative. Although details 
are provided of most of the routine tests, the significance 
of the latter is not always stressed sufficiently and the 
chapter on application of lubricating greases is more or 
less restricted to methods of dispensing. 

The author recognises that future developments in all 
aspects of the grease industry will depend largely upon 
further advances in our knowledge of grease structure. 
The section of this book dealing with structure and theory, 
while it covers the published work, deals largely with 
actual fibre structure, as determined by X-ray and electron 
microscope studies. The results, however, provide an 
incomplete picture and for this reason the succeeding 
sections of the book are still somewhat empirical. More 
recent research involving studies of the state of aggrega- 
tion of the soap fibres is producing a much clearer picture 
of the effect of processing upon grease structure 
and, equally important, the effect of structure upon 
performance.—D. EVANS. 


LECTURES ON PARTIAL DIFFERENTIAL EQUATIONS. 
I. G. Petrovsky. Translated from the Russian by A. Shenitzer. 
Interscience Publishers, New York, 1954. 245 pp. 41s. 

The book contains supplemented notes of lectures 
given by Professor Petrovsky to students of mathematics 
in the Department of Mathematics and Mechanics of the 
Moscow State University. 

The first thing about the book which impresses itself 
upon the reviewer is the excellent way in which the author 
has included both rigorous proof and physical explana- 
tion. For example, as early as page 2 the author discusses 
in clear physical terms the derivation of the equation of 
heat conduction. He then discusses physically what 
boundary and initial conditions would suffice to determine 
completely the solution to the problem, distinguishing 
clearly between the different types of possible boundary 
and initial conditions. Similarly the author considers 
vibrations of a membrane. This is immediately followed 
by rigorous considerations of the conditions under which 
the Cauchy problem has a unique analytic solution. This 
happy marriage, achieved so early in the book, continues 
to the end, making it both readable and complete. 

The book is divided clearly and conveniently into 
four chapters. Chapter 1 deals with the classification of 
equations. The conditions for the existence of a unique 
solution are discussed rigorously, and it is shown how 
various types of equation may be reduced to canonical 
form. 

Chapter 2 deals with hyperbolic equations, and con- 
tains such widely varied physical applications as (i) 
Kirchhoff’s formula for the solution of the wave equation 
(ii) the Lorentz transformation and special relativity (iii) 
vibrations of strings and membranes. It also includes such 
important mathematical concepts as (a) the “ Fourier” 
method of separation of variables and (b) general pro- 
perties of eigenfunctions and eigenvalues. 

Chapter 3 considers elliptic equations, and begins by 
defining the Dirichlet problem for a circle and solving in 
the form of Poisson’s integral. The author examines with 
great vigour the existence of a solution of the general 
Dirichlet problem, and discusses how it may be solved 
approximately by the method of finite differences. 
Physical applications deal with potential theory. 

Chapter 4 deals briefly with parabolic equations, with 
application to the conduction of heat. 

Although apparently written primarily for graduate 
students it will be valuable also to undergraduate students 


of applied mathematics, provided they are willing to pass 
over a few of the more difficult passages and concentrate 
mainly upon the applications to specific physical problems. 
The reviewer highly recommends the book to all who wish 
to further their understanding of the solution of partial 
differential equations.—N. CURLE. 


ANALYSIS OF DEFORMATION. Vol. II. K. Swainger. 
Chapman and Hall, London, 1954. 365 pp. Diagrams. 70s. 

In Volume I of this work Mr. Swainger formulated a 
linear theory to analyse the finite deformation of all 
amorphous continuous substances. General inferences 
were drawn without solving equations and boundary con- 
ditions for particular cases. Volume II is now concerned 
with the application of the theory and with the examination 
of such experimental evidence as is available for the 
equilibrium deformation of various solids having various 
boundary conditions. 

That there is at present little experimental evidence on 
the behaviour of engineering metals under complex stress 
is a point much laboured in the author’s preface. Some 
of the tests undertaken by the author to add to the experi- 
mental evidence were apparently the cause of some con- 
troversy in the case of values of the plastic transverse 
contraction ratio. Whatever the rights and wrongs of the 
controversy there remains the point that the experimental 
evidence is, as yet, sparse and there is a school of thought 
in engineering which would not venture very far with the 
application of the theory until the experimental corro- 
boration was more complete. 

Beginning with the consideration of one-stress 
definitions and tests the author proceeds through the 
various stress conditions giving theoretical examples and 
experimental results where appropriate. There are special 
chapters on rubber deformation and its theoretical inter- 
pretation, couples applied to massive bodies and the shear- 
ing of thin plane laminae. 

Vector analysis, complex variables, finite differences 
and Fourier series are conveniently covered in Appendices 
for a better understanding of the mathematics involved in 
the author’s treatise.—A.J.B. 


THE FIRST AND THE LAST. Adolf Galland. Methuen 
& Co. Ltd., 1955, 368 pp. Illustrated. 18s. 

Books about the last War, as seen from the German 
side, have been appearing in increasing numbers translated 
into English in the last year or two. Of those dealing with 
the War in the Air, “ The First and the Last ” by Lieutenant 
General Adolf Galland, the outstanding and miuch- 
decorated German fighter pilot who was the Luftwaffe’s 
“General of the Fighters Arm” from 1941 to 1945, is 
undoubtedly outstanding. 

It is of interest not only as a dramatic record of an 
operational pilot’s experiences but also as an exceptionally 
well-informed account of the higher policies of the Luft- 
waffe as they were known to a key senior officer who had 
unusual qualifications to see and interpret the problems and 
decisions of the High Command while retaining his own 
contacts with operational flying. As a result of the 
vicissitudes of high office in Germany under the National 
Socialist régime, Galland—who had started the War as a 
Lieutenant in command of a fighter squadron—finished 
the War again leading a fighter squadron into action, in 
spite of the fact that he was by then a Lieutenant General! 
This first-hand experience of the Luftwaffe’s actual job 
both at the first and the last stages of the War enables him 
to see unusually clearly the operational significance of the 
many striking technical advances in the design and 
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manufacture of aircraft and equipment with which the 
Germans sought to change the tide of war in its later stages. 

This book has claims to historical importance as a 
personal record of the Lutwaffe during the Second World 
War. But, much more than this, it is of value as an 
interpretation of the effect of technical advance on the 
technique of air warfare as it was practised during those 
years by the German Air Force. Not everybody will agree 
with all Galland’s assessments and conclusions but it is 
clear that he states a German point of view which, if it 
does not always agree with the contemporary official 
German viewpoint—indeed, is often strongly critical of it— 
is yet one which is well-informed and which, if it had 
prevailed at the time, might have had far-reaching effects 
on the progress, perhaps even on the outcome, of the War. 

The reader finds, however, that there are not a great 
many facts in this book which interested people do not 
already know. Indeed, in one or two places Galland 
himself appears to have got his facts wrong. Undoubtedly 
of most value is his outline of German aircraft development 
during those fateful years. He relates this to the contem- 
porary political and operational scene and to the practical 
experiences of an outstanding pilot so as to give old facts 
a fresh and fascinating new perspective. The story is told 
straightforwardly, without undue recrimination or excuses, 
in a way which makes it absorbingly interesting.— 
P. W. BROOKS. 


SPONTANEOUS IGNITION OF LIQUID FLUIDS. B. P. 
Mullins. Butterworths Scientific Publications, London, 1955. 
117 pp. Diagrams. 20s. 

This monograph (Agardograph No. 4) was written at 
the request of the Combustion Panel Advisory Group for 
Aeronautical Research and Development within N.A.T.O. 
and its publication has been sponsored by this group. The 
author, who is chairman of the Combustion Panel, is Head 
of the Chemical Physics Department, National Gas Turbine 
Establishment of the Ministry of Supply and has for many 
years been engaged in the study of combustion processes. 
In his preface the author states that the book “ is essentially 
a review and assessment of the state of knowledge in one 
area of the combustion field, but it should serve in addition 
as a source book of the great amount of experimental 
information on spontaneous ignition published during the 
present century.” Unfortunately, perhaps, the author has 
adhered rigidly to the terms of his assignment and presented 
a comprehensive review of the literature on methods of 
determination of spontaneous ignition temperatures with 
only tantalisingly brief reference to the applications of this 
knowledge to such problems as static fire hazards, spark 
ignition and compression-ignition engines. 

The book contains fifteen chapters, the first two deal- 
ing with General Theoretical Considerations, the next 
eight with the determination of spontaneous ignition 
temperature by the Moore, Jentzsch, Furnace, Bomb, 
Adiabatic Compression, Flow, Heated Surface and 
miscellaneous methods. Chapter eleven consists of data on 
spontaneous ignition of 433 substances which might have 
been more conveniently presented in tabular form, but 
which nevertheless represents a most important collection 
of data. The remaining chapters deal respectively with 
the effect of fuel additives, factors affecting ignition delay, 
applications and conclusions. 

The monograph is obviously intended for the specialist, 
for research workers in the fields of combustion, fuel 
properties and engine design, but it will be of considerable 
interest to other workers in related fields. For this reason 
it might have been desirable for the author to have 
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expanded the text a little in certain respects, particularly 
in the chapter on general theoretical considerations and 
in the chapter on applications. This is, however, a minor 
criticism of a most excellent monograph which within its 
one hundred pages contains a critical and at all times 
constructive review of some 1,350 scientific papers, 
Perhaps the most interesting parts of the text are those 
in which the author allows himself to depart a little from 
his position as a reviewer and points out the areas within 
the subject concerning which further knowledge is urgently 
needed. While serving as reference text to all workers in 
the field of combustion, the book will also direct and 
inspire new researches in this field. 

The book is well produced and illustrated with excellent 
line drawings. Despite the volume of work reviewed the 
text is eminently readable and well indexed. The price is 
obviously below cost and the book should find a place in 
all fuel and engine laboratories.—F. MORTON. 


HIGH SPEED AERODYNAMICS AND JET PROPULSION. 
Volume VI. General Theory of High Speed Aerodynamics. 
Edited by W. R. Sears. Oxford University Press, 1955. 758 
pp. Diagrams. 100s, 

This is the first volume to be published of a projected 
twelve volume series, usually referred to as the Princeton 
Volumes and having the general title ““ High Speed Aero- 
dynamics and Jet Propulsion.” From the point of view of 
the theoretical worker, no volume of the twelve could 
have been more welcome as a curtain raiser than this 
volume with its title “General Theory of High Speed 
Aerodynamics.” The editor, Professor W. R. Sears, and 
his very distinguished list of contributors, are to be con- 
gratulated on producing a very stimulating work. 

Books of this nature which are so fascinating to read 
present the reviewer with serious problems. It is difficult 
within the space allotted to a review to give some idea of 
the general scope and content and at the same time to 
avoid sweeping generalisations about the general present- 
ation and outlook which, let it be said at once, are no 
more uniform here than in other comparable composite 
works. 

The volume opens with an introduction by Professor 
von Karman. What more could an editor dream of than 
an introduction by Karman, and here it is fulfilling his 
dreams and mine. This is followed by a section on the 
Mathematical Aspect of Flow Problems of Hyperbole Type 
by K. O. Friedrichs which is short, to the point and very 
well written. The editor himself has written the section 
on “Small Perturbation Theory” and has produced the 
results one would expect of him when batting on his own 
ground. He has collaborated with Y. H. Quo in the section 
on “ Plane Subsonic and Transonic Potential Flows.” The 
longest section, some 218 pages, is devoted to “ Supersonic 
and Transonic Small Perturbation Theory” and has been 
written by M. A. Heaslet and H. Lomax. It is almost a 
book in itself and could, I feel, have been materially 
shortened without serious loss, indeed with some benefit 
to the reader. The section on “ Higher Approximations ” 
has been written by M. J. Lighthill. The editor singles it 
out in his preface and describes it as a piece of work 
unique in aeronautical literature and I would certainly not 
cavil at this description. The last two sections, one on the 
“ Methods of Characteristics” and the other on “ Super- 
sonic Flows and Shock Waves ” have been contributed by 
A. Ferri. I formed the impression that in the former the 
wealth of detail offered, much of which could reasonably 
be expected to be supplied by the reader, made the subject 
as well as the section much more cumbrous than it need be. 
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It will be clear from this summary that the theory 
developed is entirely inviscid. Viscous effects and turbu- 
lence are reserved for Volumes 4 and 5. 

The type of reader for whom the work is written is 
indicated by the remark in the preface that “ The need for 
a treatment of this type has long been felt in research 
institutions, universities, and private industry and _ its 
potential reflected importance in the advanced training of 
nascent aeronautical scientists has also been an important 
motivation in this undertaking.” Certainly it is clear that 
a sound mathematical background is an essential pre- 
requisite. The state of the subject is such that I for one 
would have taken exception to any other basis for a volume 
such as the one under review, though it has not escaped 
criticism on this account elsewhere. 

To one who has had to try to condense the material of 
about five of these volumes into the space occupied by 
one, the space available seems incredibly lavish. Whether 
any two people would agree precisely with the editor’s 
choice of how the space should be filled is exceedingly 
doubtful. It is difficult to avoid a feeling in places that 
the policy has been to insure against omissions by including 
everything possible. For my part I would have welcomed 
a closer sifting of the material with an eye on what is 
likely to stand the test of time, even at the risk of some 
omissions. However, this may be put down as carping 
criticism or even jealousy from one whose problem was 
so different. 

One of the advantages of the spaciousness is that it 
has been possible to say something about problems 
unsolved at the time of writing. For the research worker 
this adds a spice to the pleasures of reading this volume 
and does something to offset the earlier criticisms of 
excessive length. Perhaps I could illustrate by reference 
to the section by Sears and Kuo on “ Plane Subsonic and 
Transonic Potential Flows.” They discuss the breakdown 
of the hodograph method and the possibility of smooth 
(shock-free) transonic flows. Special solutions showing 
such flows have been obtained theoretically, but experi- 
mentally shock waves always seem to occur and Busemann 
has given a possible explanation. Sears and Kuo conclude 
by saying, ‘““The argument presented by Busemann is 
surely plausible from the physical standpoint. Nevertheless 
his conjecture has not been convincingly proved in a 
mathematical way, possibly because the formulation of 
the problem has been inadequate.” What nicer way is 
there of throwing down the gauntlet? 

Other examples may not be so clearly ticketed but they 
are there for all to see none the less. Ferri’s last section 
on “Supersonic Flows with Shock Waves,” for example, 
is a very happy hunting ground. Some discussions such 
as the flow past a yawed cone and in particular the 
associated entropy field are clearly incomplete. Others 
such as the location of the sonic point in transonic flow 
past boundaries with sudden changes of slope raise, at 
least in my opinion, as many problems as the treatment 
given answers. 

The book contains in addition a wealth of references, 
many cited at the appropriate place in the text, others 
additional to the material of the text. They form not the 
least valuable part of the book. 

In short, the book is a welcome addition to the library 
of aeronautical works. It will find, expensive though it 
is, a place on the shelves of those whose spirit is stirred 
either by the present-day achievements of the theory of 
this science or by the promise of things to come.— 
L. HOWARTH. 


FRONTIER TO SPACE. Eric Burgess. Chapman and Hall, 
London, 1955. 174 pp. Illustrated, 21s. 

The recent announcement that Government support is 
to be given for a British programme of high-altitude rocket 
research (at least up to the modest height of 60 miles) 
makes the appearance of this book very topical. Its title 
might suggest that it is yet another of the many current 
books of a popular nature on the possibilities of space 
flight, but actually it is a quite comprehensive and fairly 
technical review of achievements to date in the use of 
rockets for various kinds of physical research in the upper 
atmosphere. As such, of course, its name is none the 
less apt. 

There is one chapter largely devoted to the history and 
particulars of the rocket vehicles themselves (Aerobee, V.2, 
Viking, and so on) and another which discusses at some 
length the future possibilities of more ambitious unmanned 
missiles as artificial satellites, and as expendable “ probes ” 
for one-way interplanetary voyages, for example, to a 
reconnaissance orbit around Mars. In each case, scientific 
data would be telemetered back to the terrestrial base. 

For the rest, however, apart from a considerable number 
of very interesting illustrations, the material of the book 
is likely to be of greater value to physicists concerned 
with its main theme than to rocket engineers, for whom 
it is not principally intended. 

Present knowledge on the structure and phenomena of 
the atmosphere is surveyed and summarised. Cosmic and 
solar radiation, the propagation of radio waves as affected 
by the ionised layers at great heights, meteors, and other 
subjects of a similarly relevant nature are dealt with. So 
also are the questions of the instrumentation and tech- 
niques of telemetry employed in the research rockets which 
have been used to investigate these subjects. 


All these matters have previously been dealt with at 
greater length, by the specialists concerned, in Newell’s 


“ High-Altitude Rocket Research,’ and in “ Rocket 
Research in the Upper Atmosphere *—the proceedings 
of the 1953 Oxford Conference on this subject, 


edited by Prof. Massey. Both books are quoted as 
references in the work now under review, which features 
excellent bibliographies at each chapter ending and for 
the book as a whole. 

Mr. Burgess’s volume might perhaps best be recom- 
mended to those having a general or fringe interest in 
upper atmospheric phenomena, or to those just entering 
this field of study, as an introduction to the authoritative 
sources mentioned above. Viewed in this light, it is felt 
that it can be given a very definite recommendation, as 
the most complete review available at a reasonable price, 
within the knowledge of the present writer.—aA. V. CLEAVER. 


WHO'S WHO IN WORLD AVIATION. Edited by Wayne 
Parrish. American Aviation Publications Inc., 1955, 345 pp. 
$10. 

It is a pity that the price of this book is so high but one 
must bear in mind that it has taken two years to compile. 

Only time will prove the extent of its usefulness and 
find its omissions but 2,400 names must include most of the 
notabilities in aeronautics. One is, perhaps, a little 
surprised to find that the combined Honorary Fellowships 
of the Society and the Institute are not enough to qualify 
the former Secretary of the Society. The reparation of 
such omissions, where they occur, is promised in the next 
edition. 

A straightforward alphabetical arrangement (those 
perennial Mac people coming after Maz for a change) 
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makes the book easy to use and the biographical notes to 
most names are sufficiently detailed. All branches of 
aviation are catered for, government officials, pilot organi- 
sations, air force officers, civil executives and—as the title 
says—the “ coverage” is World-wide.—F.R.s. 


THE AIRPORT VISITOR AND AIR TRAVELLER. Com- 
piled by G. A. Frais. Penman Enterprises Ltd., 1955. 72 pp. 
Illustrated. 2s. 

A new edition of this booklet has just appeared, and 
it has several new features. A “Why and Wherefore ” 
section for air travellers, including sample lists of the 
amount of clothing which can be carried within the 44 lb. 
and 66 Ib. luggage allowance; a short article on inter- 
planetary travel; and an “ Aircraft Check Chart” in the 
form of a family tree, with civil aircraft classified according 
to their landing gear, and the number of their engines and 
fins. The descriptions of aircraft are fuller than before, 
and are classified in the same way, which should make 
identification simpler for the interested airport visitor. 


FLEXURE DEVICES—PIVOTS, MOVEMENTS, SUSPEN- 
SIONS. P. J. Geary. Published by the British Scientific 
Instrument Research Association. 44 pp. Diagrams. 10s. 6d. 

This is the first of a series of bibliographical surveys on 
instruments parts. It comprises a 13 pp. technical 


introduction with 21 diagrams, digests. of 28 selected texts 
of major importance, and annotated lists of 128 other texts 
of lesser or limited importance, together with prefatory 
matter and indexes. The period covered is from 1900 to 
date. The subject matter includes Flexure Pivots (e.g, 
cross-spring pivots); Parallel Movements; Twin Strip 
Mechanisms; Buckling Strip Devices; The Sine Spring. 

The technical introduction makes clear the advantages 
of these devices (e.g. absence of friction, backlash, 
hysteresis, wear and need for lubrication) and indicates 
their limitations. 


CONTINENTAL MILITARY AIRCRAFT. J. W. R. Taylor, 
75 pp. Illustrated. 

BRITISH MILITARY AIRCRAFT. J. W. R. Taylor. 75 pp. 
Illustrated, lan Allan Ltd., London. 2s. 6d. each. 

The Ian Allan abc series really need no review apart 
from an announcement that a new title has appeared. 
Anyone who has browsed at a Smith’s bookstall knows 
their make-up and the scope of each is explained in the 
title. Nevertheless, | could wish for a more satisfactory 
layout. Even with the explanation given in the Intro- 
duction I found it most difficult to follow the alphabetical 
order. And even though it is in service with the Royal 
Navy I must say I resent the presence of a Douglas Sky- 
raider in a book of British Military Aircraft.—F.H.s. 


Additions to the Library 


Biot, M. A. THE DIVERGENCE OF SUPERSONIC WINGS 
INCLUDING CHORDWISE BENDING. (C.A.L. Report 67.) 
Cornell Aeronautical Laboratory. 1954. 

Blough, G. O. MASTERS OF THE AIR. Smithsonian. 1954. 

Brombacher, W. G. and T. W. Lashof. BIBLIOGRAPHY 
AND INDEX ON DyYNAMIC PRESSURE MEASUREMENT. 
U.S.G.P.0. 1955. 

Bureau of Medicine and Surgery. AVIATION MEDICINE 
Practice. U.S.G.P.O. 1955. 

Chapel, C. H. et al. AIRCRAFT POWER PLANTS. McGraw- 
Hill. 1955. 

Goodman, T. R. A THEORY OF SHOCK WAVE BOUNDARY 
LAYER INTERACTION. Cornell Aeronautical Laboratory. 
1953. 

Grover, H. J. et al. THE FATIGUE OF METALS AND 
Structures. NAVAER 00-25-534. U.S.G.P.0O. 1954. 

Jeans, J. H. THe DyNamicat THEORY OF GASES (4th 
edition). Dover 1955. (C.U.P. 1925.) 

Latimer-Needham, C. H. THE FLIGHT REFUELLED 
Bomser. Flight Refuelling. 1955. 

Mahlmann, C. V. and W. M. Murray (Editors). PrRo- 
CEEDINGS: THE SOCIETY FOR EXPERIMENTAL STRESS 
ANALYSIS. Vol. XII, No. 2. S.E.S.A. 1955. 

Marquard, E. VIBRATION DYNAMICS OF THE FAST ROAD- 
VEHICLE (UNDERCARRIAGE APPLICATIONS). Cornell 
Aeronautical Laboratory. 1954. 

N.P.L. Noise MEASUREMENT TECHNIQUES. (NOTES ON 
APPLIED SCIENCE No. 10.) H.M.S.O. 1955. 

Planck, M. TREATISE ON THERMODYNAMICS (3rd edition). 
Dover 1955. (Longmans Green 1922.) 

Routh, E. J. ADVANCED DYNAMICS OF A SYSTEM OF RIGID 
Bopres (6th edition). Dover 1955. (Macmillan 1905.) 

Taylor, J. W. R. A.B.C. oF HELicopters. fan Allan 
Ltd. 1955. 

Templin, R. L. Fatiue oF ALuminiuM. (H. W. Gillett 
Memorial Lecture 1954.) A.S.T.M. 1954. 

U.S. Weather Bureau. FLYING WEATHER FORECASTS: 
How USEFUL ARE THEY? Aviation Series No. 1. 
U.S.G.P.O. 1955. 

U.S. Weather Bureau. ICE ON AIRCRAFT: ITS CAUSES 
AND EFFECTS. Aviation Series No. 2. U.S.G.P.O. 1955. 


U.S. Weather Bureau. THE JET STREAM: A BAND OF 
VERY FAST WINDS FOUND AT HIGH ALTITUDES. Avi- 
ation Series No. 3. U.S.G.P.O. 1955. 

U.S. Weather Bureau. TURBULENCE: ITS CAUSES AND 
Errects. Aviation Series No. 4. U.S.G.P.0O. 1955. 


1.A.S. Preprints 


Papers presented at the Fifth Anglo-American Aeronautical 
Conference, Los Angeles, 20th-24th June 1955. 


550 The Interaction Between Shock 
Waves and Boundary Layers D 

551 Power Control Systems for Aircraft J 

552 Design of Large Helicopters B. Kelley 

553 Power Plants for Supersonic Flight E 

554. Hypersonic Flow L 

555 Some Results of the Princeton 
University Smoke Flow Visuali- 
zation Program D. C. Hazen 

556 Stress Analysis of Multi-Web Boxes W.S. Hemp 

557 On the Behaviour of Boundary 
Layers at Supersonic Speeds 

558 Some Influences of Equipment 
Installations and Systems on Air- 
craft Design Cc 

559 Jet Noise F 

560 Low Consumption Turbine Engines A. A 

561 Combustion for Aircraft Engines W.T. Olson 

P.B 
E 


R. J. Monaghan 


562 Fatigue of Aircraft Pressure Cabins 
563 Design of High Speed Aircraft : 
564 Contributions to Fluid Mechanics I. I. Glass & 

from Shock Tube Research G.N. Patterson 
565 Airline Use of Elementary Statis- 

tical Methods in Aircraft Perform- W.C. Mentzer 

ance Measurement & F.S. Nowlan 
566 Inlet Duct-Engine Flow Compati- 

bility J.S. Alford 
567 Operating Experience with Turbo- 

Prop Aircraft P. G. Masefield 
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Reports 


AERODYNAMICS 
BOUNDARY LAYER 
see INTERNAL FLOw and WINGS AND AEROFOILS 


COMPRESSIBLE FLOW 
See WINGS AND AEROFOILS 


CONTROL SURFACES 
See STABILITY AND CONTROL 


FLuip DYNAMICS 


The fluid flow associated with the impact of liquid drops with 

solid surfaces. P. Savic and G. T. Boult. National Research 

Council of Canada. Report No. MT-26. 1955. 
A theory is presented of the fluid flow associated with the 
impact of liquid drops with a solid surface. The shape of 
the spreading drop is calculated as well as the pressure 
distribution over the impact plane. Experiments carried out 
with high speed spark camera bear out the main features 
predicted by the theory.—(1.4 x 1.9.1). 


Contribution a la théorie cinétique des écoulements unidimen- 
sionnels. J. J. Bernard. O.N.E.R.A, Publication No. 74 (1955). 


(1.4). 


Variation of local liquid-water concentration about an ellipsoid 

of fineness ratio 10 moving in a droplet field. Rinaldo J. Brun 

and Robert G. Dorsch. N.A.C.A. Tech. Note 3410. April 1955. 
Trajectories of water droplets about an ellipsoid of revolu- 
tion (10 per cent. thick) in flight through a droplet field were 
computed with the aid of a differential analyser. Analyses 
indicate that the local concentration of liquid water at 
various points about an ellipsoid varies considerably and 
under some conditions may be several times the free-stream 
concentration. Curves of the local concentration factor as a 
function of spatial position are presented in terms of dimen- 
sionless parameters that describe flight and atmospheric 
conditions.—(1.4). 


INTERNAL FLOW 


Compressor cascade flutter tests. 20° camber blades, medium 

and high stagger cascades. D. A. Kilpatrick and J. Ritchie. 

A.R.C. C.P. No. 187. (December 1953, published 1955). 
The tests show the existence of three main zones of flutter 
at high stress; stalling flutter, shock-stalling flutter and chok- 
ing flutter. These zones are similarly located (with reference 
to the aerodynamic characteristics) for both the medium and 
the high stagger cascades tested, and they extend over a wide 
range of incidence. Good correlation between the zones of 
flutter and the experimentally measured blade force deriva- 
tives, with respect to Mach number and incidence has been 
obtained. More experimental data are, however, required 
before a quantitative analysis of the problem can be 
achieved.—(1.5.4.2). 


Theory of the jet syhon. B. Szczeniowski. N.A.C.A. T.N. 

3385 (May 1955). 
A new approach to the theory of the mixing of two currents 
in an injector is presented dealing with an incompressible 
ideal fluid. The theory shows new potentialities in an 
appropriate shaping of the form of the walls of the mixing 
zone so as to improve the jet-syphon efficiency beyond that 
heretofore predicted theoretically. Some examples of ways 
to improve jet-syphon efficiency are indicated.—(1.5.1). 


An investigation of several N.A.C.A. 1-series nose inlets with 

and without protruding central bodies at high-subsonic Mach 

numbers and at a Mach number of 1:2. R. E. Pendley and 

H. L. Robinson. N.A.C.A., T.N. 3436 (May 1955). 
Measurements of pressure distribution drag, and internal 
flow pressure loss are presented for three nose inlets, two 
of which were fitted with protruded central bodies. Test 
Mach number and inlet-velocity ratio ranged from 0:4 to 1-2 
and from 0 to 1:34 respectively. The nose-inlet pressure 
drag at a Mach number of 1:2 and the central body effects 
on subcritical drag, the supercritical drag rise, and the inlet 
total pressure loss are discussed.—(1.5.1.1). 


Preliminary investigation of a submerged air scoop utilizing 
boundary-layer suction to obtain increased pressure recovery. 
M. R. Nichols and P. K. Pierpont N.A.C.A., T.N. 3437 (April 
1955). 
Results are given of low-speed tests of a submerged inlet 
consisting essentially of a conventional scoop located in a 
dimple in the fuselage surface. Boundary layer control 
systems investigated are shown to provide important 
increases in performance. It appears that the flow instability 
frequently encountered in the case of twin internally coupled 
inlets will be avoided for design high speed inlet-velocity 
ratios as low as 0°5.—(1.5.1 X 1.1.5). 


Loaps 
See STABILITY AND CONTROL and TESTING AND INSTRUMENTS 


STABILITY AND CONTROL 

Low-speed tunnel model tests on tailplane rolling moments in 

sideslip. A. Spence, J. W. Leathers and D. A. Kirby. R. & M. 

2941 (October 1951, published 1955). 
To provide data for stressing purposes, measurements have 
been made of the effect of sideslip on the rolling moment on 
a 41:5° swept-back tailplane mounted at three heights on the 
fin of a model of a single jet aircraft with a 40° swept-back 
wing. Incidence and tailplane setting were varied, and the 
effects of rudder deflection were obtained with the tailplane 
at the top of the fin. Brief results on a delta aircraft model 
with a delta tailplane at the top of the fin are also included. 
—(1.8.1 1.6.1 x 1.3.3). 


Matrix methods for determining the longitudinal-stability 
derivatives of an airplane from transient flight data. J. J. 
Donegan. N.A.C.A. Report 1169. (1954). 
Three methods, illustrated by several examples using flight 
data, are presented for calculating the longitudinal stability 
derivatives from transient flight data. The results indicate 
the scatter which is typical of this type of analysis.—(1.8.2.1), 


WINGS AND AEROFOILS 


Some visual observations of the flow over a swept-back wing 
in a water tunnel, with particular reference to high incidences. 
N.C. Lambourne and P. §. Pusey. A.R.C., C.P. No. 192 (May 
1954, published 1955). 
A qualitative investigation has been made consisting of 
visual observations of the flow over a small sweptback wing. 
The tests were restricted to low Reynolds numbers (about 
10° based on mean chord). Visualisation was effected by 
(i) the addition of small aluminium particles to the water, 
and (ii) the introduction of air into the stream close to the 
model. In addition, the surface patterns produced by the 
flow after the wing had been coated with oil were obtained. 
Attention has been mainly concentrated on the type of flow 
that exists at high incidence—(1.10.1.2). 


The zero-lift wave drag of a particular family of unswept, 

tapered wings with linearly varying thickness ratio. A. Hender- 

son and J. M. Goodwin. N.A.C.A., T.N. 3418 (May 1955). 
On the basis of linear theory, the zero-lift wave drag of a 
particular family of unswept, tapered wings with linearly 
varying thickness ratio and symmetrical parabolic-arce 
sections has been calculated. By comparing the drag for 
these wings with that for a corresponding constant-thickness 
ratio wing with rhombic sections, it is found that the 
variable-thickness ratio wings can be used to advantage with 
no serious structural penalties if the wings are assumed to 
have the same given root thickness ratio or the same internal 
volume.—(1.10.1.3 x 1.2.3). 


Exploratory investigation of an airfoil with area suction applied 
to a porous, round trailing edge fitted with a lift-control vane. 
R. E. Dannenberg and J. A. Weiber. N.A.C.A., T.N. 3498 
(April 1955). 
A two-dimensional wind tunnel investigation has been made 
of an N.A.C.A. 651-012 aerofoil modified with a porous, 
round trailing edge fitted with a small vane. Area suction 
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was applied and the vane was used to fix the rear stagnation 
point. The effects of variations of the chordwise extent of 
suction and of the vane geometry were investigated with the 
aim of attaining high maximum lift with low suction 
quantity.—(1.10.2.1 x 1.1.5). 


Optimum two-dimensional multiplanes in supersonic flow. 

R. M. Licher. Douglas Report No. SM-18688 (January 1955). 
Biplanes and other multiplanes are studied by linear theory 
to find optimum configurations having minimum drag due 
to lift in supersonic flight. It is shown that suitable camber 
can often reduce the drag of a two-dimensional biplane. It 
is sometimes possible to further reduce the drag below that 
of the cambered biplanes by introduction of a thickness dis- 
tribution. A general criterion is developed to indicate when 
the optimum distribution of camber and thickness for mini- 
mum drag due to lift has been obtained. Distributions of 
camber and thickness having no net lift or drag, called zero 
loadings, are discussed. It is shown how zero loadings may 
be combined with optimum distributions of camber and 
thickness to provide more depth for structure without 
increasing the drag.—(1.10.1.1 x 1.2.3). 


Notes on the variation of drag with Mach number of a Buse- 
mann biplane. B. J. Beane. Douglas Report No. SM-18737 
(January 1955). 
It has been observed that the Busemann biplane, at zero lift, 
has zero wave drag at infinitely many Mach numbers. This 
result prompted further study of the wave drag of this con- 
figuration, with the object of finding the wave drag variation 
with Mach number.—(1.10.1.1 x 1.2.3). 


A study of interference and lift transfer for a staggered delta- 
wing biplane in supersonic flow. B.J. Beane. Douglas Report 
No. SM-18738 (January 1955). 
The problem of interference and lift transfer in supersonic 
flow is studied for a biplane having identical delta plan form 
wings. The wings are staggered so that there is interference 
at the upper wing only, while the lower wing is forward of 
the region of influence of the upper wing. Curves are 
included which show the downwash induced by the lower 
_—t along the centre line of the upper wing.—(1.10.1.1 x 


TESTING AND INSTRUMENTS 


The calculation of the wave drag of an arbitrary slender body 

by means of an electrical analogy tank. P. J. Pocock. N.A.E. 

Laboratory Report LR-127 (March 1955). 
In view of the analogy between the supersonic wave drag of 
a slender body and the subsonic induced drag of an 
“equivalent” lifting line the electric analogy tank for 
solving the Prandtl lifting line equations can be used to 
calculate the wave drag of arbitrary slender bodies. A 
Malavard type electric analogy tank at the N.A.E. is 
described and its use illustrated by examples.—(1.12.6 x 1.6.1). 


The N.R.C. icing wind tunnels and some of their problems. 

C. K. Rush. N.A.E. Laboratory Report LR-133 (April 1955). 
The 44 ft. and blower icing tunnels of the N.R.C. are 
described. Some of the difficulties encountered due to the 
icing function and the methods of dealing with them are 
discussed. The nature of tests undertaken in the tunnels is 
outlined and a few typical tests described.—(1.12.1.1). 


AIRCRAFT OPERATION 


Theoretical analysis of light plane landing and take-off accidents 

due to encountering the wakes of large airplanes. Z. O. Bleviss. 

Douglas Report No, SM-18647 (December 1954). 
It is shown in this analysis that many accidents to light 
aircraft during landing and take-off are caused by the rota- 
tional field of air associated with the trailing vortices from 
the wings of nearby large aircraft. Detailed numerical 
results are given for the case of a 40 ft. span light aeroplane 
with a DC-6B as the large aircraft —(5.3). 


AIRPORTS 


A T.V.O.R. automatic-rebroadcast system. T. §. Wonnell and 
B. H. Boyle. C.A.A Technical Development Report No. 232 
(May 1955). 
This report describes the development and testing of a 
T.V.O.R. automatic rebroadcast system developed at the 
Technical Development and Evaluation Centre. Indianapolis, 
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Indiana, The system provides a means whereby pilots can 
communicate with one another over airports which have no 
control tower functioning.—(6.5). 


Analytical and simulation studies of terminal-area air traffic 

control. S. M. Berkowitz and E. L. Fritz. C.A.A. Technical 

Development Report No. 251 (May 1955). 
This report summarises the programme of terminal-areg 
simulation work completed jointly at the Franklin Institute 
Laboratories for Research and Development and at the 
Technical Development and Evaluation Centre of the Civil 
Aeronautics Administration during a nine-month period, 
The specific tasks involved the set-up, description, and 
analysis of the mathematical and physical system analogs and 
a determination of the inter-action of the many variables 
involved. Advanced radar-vectoring techniques previously 
developed to a high degree by T.D.E.C. were employed. Ag 
a matter of convenience the Washington National Airport 
terminal area, which is 40 miles in diameter, was used in 
establishing certain generalised problems with many specific 
applications and solutions. Representative high-density 
traffic samples of either arrivals only or of arrivals and 
departures, with peak durations of one to four hours, were 
employed. The experiments included theoretical acceptance 
rates ranging from 30 to 41 landings per hour and used 
random arrival rates of 35 aircraft per hour and peak 
departure rates of 31 aircraft per hour. Three methods of 
analysing the problems were used; an F.I.L. ideal technique, 
an FIL. graphical space-time technique, and the T.D.E.C, 
dynamic simulator. Correlation among the three was found 
to be excellent for the wide variety of traffic samples and 
of conditions tested.—({6.5 x 26). 


INSTRUMENTS AND EQUIPMENT 


Moving coil vibration generators and pick-ups. H. M., Nelson. 

A.R.L./ 1.54 (June 1954). 
The equations for an ideal vibration generator or pick-up 
are derived and the modifications required to allow for a 
force/current phase lag or force/displacement phase lag in 
the spring are indicated. A method is indicated whereby 
the vibrator constants under working conditions can be 
obtained by measuring the electrical impedance of the 
vibrator when locked, free and with added mass.—(18.1). 


self-excited, alternating-current, constant-temperature hot- 

wire anemometer. C. E. Shepard. N.A.C.A., T.N. 3406 (April 

1955). 
The hot-wire anemometer described herein was used to 
study turbo-jet-engine compressor rotating stall and surge. 
The system was capable of measuring flow over a frequency 
range of zero to 500 cycles per second. A self-sustaining, 
8 kilocycle-per-second oscillation was used to heat the wire 
in a constant temperature system. The relatively rugged, 
large-diameter wire allowed the use of the anemometer 
during the full-scale performance testing of compressors and 
turbo-jet engines. The commercial audio-frequency amplifier 
used reduced the cost and the time required to build the 
anemometer.—(18.1). 


MATERIALS 


A method for assessing the corrosion of iron in soils. Part 1. 

Theory of the soil probe. K. F. Lorking. A.R.L./Met.2 

(December 1954). 
A detailed explanation is given of corrosion processes 
occurring on iron in aqueous media. A description is given 
of the theoretical basis of methods for the interpretation of 
potentials assumed by iron and by platinum electrodes. An 
account is then given of a soil probe designed to assess rates 
of corrosion of iron in soils.—(21.2.1 x 21.6.1 x 21.6.2). 


FATIGUE 
See also STRUCTURES 
Scatter in fatigue life of 24S-T alclad specimens with drilled 
holes. W. Weibull. §.A.A.B. T.N. 32 (May 1955). 


The fatigue lifetimes of 1,088 flat 24S-T Alclad specimens 
with two drilled holes have been determined at a tension 
pulsating between 14 kg. sq. mm. and zero.—(31.2.2.3.2.1.8). 
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